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

    林芝市農(nóng)業(yè)面源污染負(fù)荷時(shí)空變化與分布特征

    2017-08-11 08:46:40郝守寧付意成
    關(guān)鍵詞:林芝地區(qū)面源污染源

    郝守寧,付意成

    (1.西藏農(nóng)牧學(xué)院,西藏林芝 860000;2.中國(guó)水利水電科學(xué)研究院,北京 100038)

    林芝市農(nóng)業(yè)面源污染負(fù)荷時(shí)空變化與分布特征

    郝守寧1,付意成2*

    (1.西藏農(nóng)牧學(xué)院,西藏林芝 860000;2.中國(guó)水利水電科學(xué)研究院,北京 100038)

    為分析西藏林芝地區(qū)的農(nóng)業(yè)面源污染形勢(shì)與特征,基于林芝市2000—2015年的農(nóng)業(yè)統(tǒng)計(jì)資料,應(yīng)用輸出系數(shù)模型估算了該地區(qū)農(nóng)業(yè)面源污染總氮(TN)、總磷(TP)年輸出負(fù)荷。并對(duì)農(nóng)村生活污水、農(nóng)田化肥、畜禽養(yǎng)殖三類污染源分別進(jìn)行了分析。結(jié)果表明.TN、TP均呈現(xiàn)逐年平穩(wěn)上升的趨勢(shì);農(nóng)業(yè)面源污染高負(fù)荷區(qū)主要集中在工布江達(dá)縣、波密縣和林芝縣,而農(nóng)業(yè)面源污染負(fù)荷強(qiáng)度最高的行政區(qū)為米林縣;農(nóng)業(yè)面源污染負(fù)荷的三大營(yíng)養(yǎng)源貢獻(xiàn)中畜禽養(yǎng)殖輸出負(fù)荷的貢獻(xiàn)率隨時(shí)間的推移有明顯增加的趨勢(shì),90%的農(nóng)業(yè)面源污染負(fù)荷來(lái)源于畜禽養(yǎng)殖和農(nóng)田化肥;林芝市農(nóng)牧業(yè)發(fā)達(dá)地區(qū)農(nóng)業(yè)面源污染形勢(shì)不容樂觀,須大力開展防控工作,以改善該地區(qū)農(nóng)業(yè)面源污染日趨嚴(yán)重的狀況。

    農(nóng)業(yè)面源;總氮;總磷;林芝地區(qū);分布特征

    農(nóng)業(yè)面源污染是由沉積物、廢料、農(nóng)藥、致病菌等不確定污染源引起的對(duì)水體、河流、湖泊、大氣等生態(tài)系統(tǒng)的污染[1]。農(nóng)業(yè)面源污染的發(fā)生沒有固定的排污口,其污染源主要借助降雨或排水過程將地表存留污染物帶走導(dǎo)致水體富營(yíng)養(yǎng)化[2],危及飲用水安全[3],使土壤板結(jié)[4],破壞生物多樣性和生態(tài)系統(tǒng)平衡[5]等。與點(diǎn)源污染相比,我國(guó)農(nóng)業(yè)面源污染具有排放分散、隱蔽、隨機(jī)、復(fù)雜、不易檢測(cè)等特點(diǎn),加之單位面積上的污染負(fù)荷較小,從而導(dǎo)致人們往往忽視農(nóng)業(yè)面源污染的宏觀效應(yīng),已有的研究成果顯示,農(nóng)業(yè)面源污染是導(dǎo)致國(guó)內(nèi)外區(qū)域水環(huán)境污染的重要原因之一[6]。近幾年隨著西藏經(jīng)濟(jì)快速發(fā)展,城鎮(zhèn)化步伐加快,人口及牲畜養(yǎng)殖量增加,農(nóng)業(yè)生產(chǎn)面臨著嚴(yán)峻的資源緊張和生態(tài)環(huán)境問題,農(nóng)業(yè)面源污染對(duì)西藏生態(tài)安全屏障建設(shè)的影響也日益突出。因此,研究農(nóng)業(yè)面源污染的來(lái)源、輸出負(fù)荷,并據(jù)此制定控制對(duì)策與防治措施具有重大意義。

    國(guó)內(nèi)外很多科研工作者對(duì)流域農(nóng)業(yè)面源污染進(jìn)行了一系列的研究和探討,面源污染負(fù)荷的估算方法主要包括監(jiān)測(cè)法[7]、輸出系數(shù)法[8-9]、排污系數(shù)法[10]以及模型法[11]等,估算模型的選定取決于研究區(qū)域所搜集資料的完整程度。很多國(guó)外開發(fā)的面源污染負(fù)荷計(jì)算模型的參數(shù)校準(zhǔn)和驗(yàn)證較為復(fù)雜,國(guó)內(nèi)多數(shù)流域和地區(qū)缺乏長(zhǎng)時(shí)間序列監(jiān)測(cè)資料,從而導(dǎo)致很多模型在國(guó)內(nèi)應(yīng)用受到局限[12]。輸出系數(shù)法所需參數(shù)較少、應(yīng)用性強(qiáng),特別適合在缺乏資料的條件下采用[13],被國(guó)內(nèi)學(xué)者應(yīng)用于很多研究區(qū)域,在三峽庫(kù)區(qū)大寧河流域[14]、大遼河流域[15]、長(zhǎng)江流域[16-17]、鄱陽(yáng)湖流域[18]、呼蘭河流域[19]等地區(qū)都得到了較好的應(yīng)用效果。本研究以林芝地區(qū)為研究區(qū)域,采用輸出系數(shù)模型,對(duì)林芝地區(qū)2000—2015年的農(nóng)村生活污水、農(nóng)田化肥、畜禽養(yǎng)殖三類污染源年總氮(TN)、總磷(TP)負(fù)荷進(jìn)行估算,并探討研究區(qū)域內(nèi)農(nóng)業(yè)面源污染負(fù)荷的時(shí)空分布特征與變化趨勢(shì),為林芝地區(qū)農(nóng)業(yè)面源污染的防治規(guī)劃和措施的制定提供科學(xué)依據(jù)。

    1 研究區(qū)概況

    林芝市處于西藏東南部,雅魯藏布江中下游流域,該地區(qū)氣候宜人,降雨量豐富,氣象資料顯示年均降雨量接近800 mm,素有“西藏江南”之美稱,市中心所在地八一街道距離自治區(qū)首府拉薩市400多公里。東與昌都地區(qū)和云南省交界,西與拉薩市和山南地區(qū)相連,北與那曲地區(qū)毗鄰,南與緬甸、印度兩國(guó)接壤,東西長(zhǎng)646.7 km,南北寬353.2 km,邊境線長(zhǎng)約1 006.5 km,幅員面積11.7萬(wàn)km2。平均海拔3000 m,而最低處卻只有152 m,是世界陸地垂直地貌落差最大的地帶,平均海拔低于西藏其他地區(qū)。地區(qū)轄林芝、工布江達(dá)、朗縣、米林、波密、察隅和墨脫7個(gè)縣,居民除藏族以外,還有門巴族、珞巴族、怒族、獨(dú)龍族、傈僳族、納西族、白族以及未確認(rèn)民族成份的僜人。畜牧業(yè)與農(nóng)業(yè)是林芝市農(nóng)牧民的主要經(jīng)濟(jì)來(lái)源,截至2015年末,林芝市的總?cè)丝跀?shù)為21.24萬(wàn),其中農(nóng)業(yè)人口13.89萬(wàn),占總?cè)丝诘?5%;耕地面積1.8萬(wàn)hm2,菜地面積1 465.41 hm2,草場(chǎng)面積1 032.76 hm2;養(yǎng)殖大牲畜41.7萬(wàn)頭、豬24.3萬(wàn)頭、羊6.6萬(wàn)只、家禽26.2萬(wàn)只。受地形地貌條件的影響,區(qū)內(nèi)居民傍河而居,依河而耕地,耕地灌溉施肥以氮肥和磷肥為主,區(qū)內(nèi)牲畜以行政村為單位實(shí)行規(guī)模化養(yǎng)殖。林芝市國(guó)土資源局?jǐn)?shù)據(jù)顯示.區(qū)內(nèi)未利用土地面積與常年積雪的高山林地面積之和占全區(qū)土地總面積的75.23%,耕地面積占全區(qū)土地面積的0.17%左右,而居民用地僅為4453 hm2,研究區(qū)內(nèi)人口及養(yǎng)殖密度雖低,但伴隨西部開發(fā)的加速、旅游業(yè)的發(fā)展,加之區(qū)內(nèi)人口、牲畜的不斷增加,林芝已成為以畜牧業(yè)為主導(dǎo)產(chǎn)業(yè)的生態(tài)脆弱區(qū),農(nóng)業(yè)面源污染是林芝市生態(tài)環(huán)境的一個(gè)極大挑戰(zhàn)。

    2 材料及方法

    2.1 輸出系數(shù)模型

    隨著面源污染對(duì)人類的影響越來(lái)越大,人們對(duì)面源污染的研究也在不斷深入。輸出系數(shù)模型在早期的應(yīng)用是假定所有土地利用類型的輸出系數(shù)都相等[20],這種建立在假設(shè)的基礎(chǔ)上計(jì)算出的結(jié)果往往與實(shí)際情況相差較大,經(jīng)過長(zhǎng)期的探索研究,1996年Johnes等[20]在早期輸出系數(shù)模型的基礎(chǔ)上進(jìn)行了改進(jìn),對(duì)農(nóng)業(yè)人口生活污水的排放及牲畜的輸出系數(shù)的確定方法進(jìn)行了改進(jìn),同時(shí)考慮植物的固氮、氮的空氣沉降等因素,提高了模型在流域面源污染負(fù)荷研究中的靈敏性。在國(guó)內(nèi),蔡明等[21]對(duì)該輸出系數(shù)模型做了進(jìn)一步的改進(jìn),應(yīng)用中充分考慮水文因素的影響,使其對(duì)區(qū)域不同年份的TN、TP負(fù)荷估算結(jié)果更加準(zhǔn)確,本次研究采用的改進(jìn)輸出系數(shù)模型考慮了水文因素的影響,輸出系數(shù)模型如下.

    式中:Lj為營(yíng)養(yǎng)物質(zhì)j的總流失量;α為降雨影響因子,即研究區(qū)域內(nèi)降雨量影響系數(shù);Eij為污染物j在第i種土地利用類型中的輸出系數(shù)或第i種畜禽的排放系數(shù)或人口的輸出系數(shù);Ai為第i類土地利用類型的面積或第i種牲畜的數(shù)量或人口數(shù)量;P為降雨輸入的營(yíng)養(yǎng)物質(zhì)總量,本文主要對(duì)不同農(nóng)業(yè)面源污染源的負(fù)荷進(jìn)行比較,未考慮此項(xiàng)的影響。

    2.2 輸出系數(shù)選取

    輸出系數(shù)選用是輸出系數(shù)模型的核心,也是最重要、最難以確定的參數(shù),輸出系數(shù)選取是否合理與準(zhǔn)確直接影響研究區(qū)域面源污染負(fù)荷的估算精確度[22]。林芝地區(qū)位于雅魯藏布江中下游流域,受藏族生活風(fēng)俗的影響,區(qū)內(nèi)不存在水產(chǎn)養(yǎng)殖,根據(jù)模型的選取及流域特點(diǎn)將流域農(nóng)業(yè)面源污染源分為農(nóng)村生活污水、農(nóng)田種植流失、畜禽養(yǎng)殖流失三大類。各污染源估算過程中不同的土地利用狀況、農(nóng)村人畜數(shù)量、化肥使用量以及作物種植方式,會(huì)導(dǎo)致選用不同的輸出系數(shù),對(duì)牲畜的輸出系數(shù)選用依據(jù)其種類、數(shù)量和分布情況進(jìn)行確定,人口輸出系數(shù)則主要依據(jù)生活污水的排放和處理狀況進(jìn)行選定。農(nóng)村廢水一般多采用化糞池進(jìn)行收集和處理,作為農(nóng)家肥使用,故此估算過程中的參數(shù)均采用化糞池處理后的排放系數(shù)。根據(jù)《第一次全國(guó)污染源普查城鎮(zhèn)生活源產(chǎn)排污系數(shù)手冊(cè)》中的分區(qū),以及城鎮(zhèn)生活污染源污染物排放系數(shù),結(jié)合林芝市農(nóng)村生活特點(diǎn)進(jìn)行相應(yīng)的系數(shù)調(diào)整。農(nóng)業(yè)種植習(xí)慣、降雨、地形等因素導(dǎo)致不同區(qū)域的農(nóng)田化肥流失存在較大的差異,本研究中通過查詢文獻(xiàn)數(shù)據(jù)庫(kù)、研究報(bào)告等資料,結(jié)合研究區(qū)域的相關(guān)農(nóng)田徑流監(jiān)測(cè)數(shù)據(jù),參考《第一次全國(guó)污染源普查農(nóng)業(yè)污染源肥料流失系數(shù)手冊(cè)》的研究成果,建立農(nóng)田化肥流失負(fù)荷估算方法,估算過程中,綜合考慮流域內(nèi)主要耕地類型、種植類型、耕作方式、土壤類型、地形地貌對(duì)參數(shù)選擇的影響。區(qū)內(nèi)畜禽養(yǎng)殖是主導(dǎo)產(chǎn)業(yè),且養(yǎng)殖方式是以行政村為單位的規(guī)模養(yǎng)殖,依據(jù)畜禽污染物排泄系數(shù)及畜禽糞便處理狀況估算畜禽養(yǎng)殖污染物流失量。依據(jù)上海環(huán)保局2000年對(duì)集約化畜禽養(yǎng)殖場(chǎng)污染情況的調(diào)查結(jié)果,及國(guó)家環(huán)??偩衷?000年對(duì)全國(guó)23個(gè)省、自治區(qū)、直轄市規(guī)模化畜禽養(yǎng)殖業(yè)污染情況的調(diào)查結(jié)果,確定研究區(qū)域內(nèi)輸出系數(shù),見表1和表2。

    2.3 數(shù)據(jù)獲取

    利用輸出系數(shù)模型對(duì)農(nóng)業(yè)面源污染負(fù)荷進(jìn)行估算時(shí),除了采用研究區(qū)域輸出系數(shù)數(shù)據(jù)外,還需要各污染源的統(tǒng)計(jì)數(shù)據(jù)和研究區(qū)的空間數(shù)據(jù),包含研究區(qū)域內(nèi)的不同土地利用類型面積、畜禽養(yǎng)殖年存欄量以及農(nóng)業(yè)人口數(shù)量??臻g數(shù)據(jù)來(lái)源于中國(guó)科學(xué)院地理科學(xué)與自然資源研究所的縣級(jí)行政區(qū)矢量數(shù)據(jù)。通過查閱2001—2016年林芝市統(tǒng)計(jì)年鑒,受當(dāng)?shù)刭Y料條件限制,整理匯集了林芝市8年的統(tǒng)計(jì)數(shù)據(jù),如表3所示。

    3 結(jié)果與分析

    3.1 各縣污染負(fù)荷強(qiáng)度估算

    林芝地區(qū)7個(gè)縣的人口密集程度、交通狀況以及土地利用結(jié)構(gòu)各不相同,加之各縣經(jīng)濟(jì)發(fā)展水平的差異,導(dǎo)致林芝地區(qū)7個(gè)縣內(nèi)的污染負(fù)荷強(qiáng)度也不盡相同?;谘芯繀^(qū)域的統(tǒng)計(jì)數(shù)據(jù),應(yīng)用輸出系數(shù)模型模擬計(jì)算,同時(shí)對(duì)比分析7個(gè)縣的氮、磷負(fù)荷強(qiáng)度,結(jié)果見圖1。

    表1 畜禽糞便、農(nóng)村居民生活污水流失率(%)Table 1 Loss rate of sewage from livestock breeding and rural residents(%)

    表2 不同土地利用類型流失系數(shù)Table 2 Loss coefficients of various land utilization modes

    表3 林芝地區(qū)2000—2015年農(nóng)業(yè)土地利用、畜禽養(yǎng)殖及人口情況Table 3 Agricultural land area,livestock breeding and population from 2000 to 2015 in Nyingchi

    2000—2015年期間林芝地區(qū)7個(gè)縣的農(nóng)業(yè)面源污染TN、TP年負(fù)荷存在顯著的空間分異特征,各縣年負(fù)荷的變化也存在差異,污染高負(fù)荷區(qū)主要集中在工布江達(dá)縣、波密縣、林芝縣等3個(gè)縣級(jí)行政區(qū),這與林芝地區(qū)的農(nóng)業(yè)人口及其畜禽養(yǎng)殖分布基本相符。該3個(gè)縣級(jí)行政區(qū)內(nèi)人口占研究區(qū)總農(nóng)業(yè)人口的54.6%,且農(nóng)業(yè)人口分布相對(duì)集中,同時(shí)該3個(gè)縣的草場(chǎng)生產(chǎn)較為發(fā)達(dá),畜禽養(yǎng)殖量較其他縣多,總體來(lái)說農(nóng)業(yè)面源污染TN、TP年輸出負(fù)荷與農(nóng)業(yè)人口數(shù)量及畜禽養(yǎng)殖量有著密切的關(guān)系。拉林鐵路及拉林高速公路的修建,更為3個(gè)縣的交通提供了便利,農(nóng)牧民的生活水平也隨之提高,大量農(nóng)村生活污水及糞便缺少基礎(chǔ)設(shè)施治理對(duì)策,加之散養(yǎng)畜禽隨時(shí)隨地的排泄,造成此區(qū)域TN、TP高負(fù)荷分布;其次3個(gè)縣占有研究區(qū)42.8%的耕地面積,研究區(qū)內(nèi)農(nóng)業(yè)生產(chǎn)活動(dòng)發(fā)達(dá),土地利用率高,化肥、農(nóng)藥施用量的逐年增加及水土流失等原因使耕地內(nèi)大量未被吸收的氮、磷嚴(yán)重流失也是造成該研究區(qū)內(nèi)TN、TP高負(fù)荷分布的因素。墨脫縣是研究區(qū)內(nèi)TN、TP負(fù)荷分布最低的行政區(qū),主要原因是墨脫縣地處世界第一的雅魯藏布大峽谷的深處,至今是全中國(guó)最后一個(gè)不通公路的縣,農(nóng)業(yè)人口數(shù)量、耕地面積及畜禽養(yǎng)殖量遠(yuǎn)遠(yuǎn)低于研究區(qū)內(nèi)的其他行政區(qū)。就三類污染源來(lái)講,林芝地區(qū)畜禽養(yǎng)殖營(yíng)養(yǎng)源產(chǎn)排比較大,相對(duì)于整個(gè)林芝地區(qū)農(nóng)業(yè)面源污染負(fù)荷而言,拉林高速路沿途的行政區(qū)是整個(gè)林芝地區(qū)農(nóng)業(yè)面源污染物年輸出較為嚴(yán)重的區(qū)域。

    圖1 林芝地區(qū)各縣TN、TP負(fù)荷對(duì)比Figure 1 Contrast on loads of TN and TP in different counties in Nyingchi

    3.2 TN、TP年輸出負(fù)荷空間分布變化

    3.2.1 不同污染源TN、TP的年輸出負(fù)荷變化

    本文研究區(qū)的農(nóng)業(yè)面源污染源由農(nóng)村生活污水、農(nóng)田化肥流失、畜禽養(yǎng)殖三類營(yíng)養(yǎng)源組成,將研究區(qū)域的統(tǒng)計(jì)數(shù)據(jù)及表1、表2的輸出系數(shù)帶入公式(1),并將單位換算成t,估算林芝地區(qū)2000—2015年各污染源的TN、TP年輸出負(fù)荷,見表4。結(jié)果表明,2000—2015年TN、TP年輸出負(fù)荷具有相似的規(guī)律.TN的年輸出總負(fù)荷明顯地呈逐年增加的趨勢(shì);TP的年輸出總負(fù)荷雖然從2011年之后有減少和增加的波動(dòng),但從長(zhǎng)時(shí)間段看TP的年輸出負(fù)荷總體上也呈上升的趨勢(shì)。三類不同營(yíng)養(yǎng)源中,農(nóng)村生活污水、畜禽養(yǎng)殖兩大營(yíng)養(yǎng)源對(duì)TN、TP的貢獻(xiàn)量雖差別較大但具有相同的規(guī)律,貢獻(xiàn)量均隨著時(shí)間的推移而增加;營(yíng)養(yǎng)源畜禽養(yǎng)殖流失對(duì)研究區(qū)TN的貢獻(xiàn)量與另兩類營(yíng)養(yǎng)源存在同樣的規(guī)律,而對(duì)研究區(qū)TP的貢獻(xiàn)量總體上雖也為上升趨勢(shì),但期間有下降的年份,這與研究區(qū)內(nèi)隨著市場(chǎng)需求而改變耕地種植類型有直接關(guān)系。2015年研究區(qū)內(nèi)營(yíng)養(yǎng)源畜禽養(yǎng)殖流失與農(nóng)田化肥流失對(duì)TN、TP年輸出負(fù)荷總和分別是農(nóng)村生活廢水貢獻(xiàn)量的12倍與25倍,是林芝地區(qū)的主要營(yíng)養(yǎng)源,也是日后研究區(qū)內(nèi)整治的主要對(duì)象。

    3.2.2 不同污染源TN、TP貢獻(xiàn)率的空間變化

    農(nóng)業(yè)面源污染負(fù)荷的分布及空間變化與研究區(qū)內(nèi)土地利用類型、農(nóng)村人口數(shù)量及畜禽養(yǎng)殖量等密切相關(guān),因此農(nóng)業(yè)面源污染負(fù)荷的產(chǎn)排具有很強(qiáng)的空間性。本研究借助GIS空間分析平臺(tái),將三類營(yíng)養(yǎng)源對(duì)TN、TP的貢獻(xiàn)量估算成果添加到相應(yīng)行政區(qū)劃圖的屬性列表中,然后利用GIS強(qiáng)大的空間分析功能得出不同污染源對(duì)研究區(qū)內(nèi)TN、TP貢獻(xiàn)率的空間變化專題圖,非常直觀地表達(dá)出研究區(qū)內(nèi)不同污染源的時(shí)空分布及其變化特征,不同營(yíng)養(yǎng)源對(duì)TN貢獻(xiàn)率見圖2,對(duì)TP貢獻(xiàn)率見圖3。

    表4 林芝地區(qū)2000—2015年各污染源TN、TP年輸出負(fù)荷(t)Table 4 Annual emission load of TN/TP from various pollution sources from 2000 to 2015 in Nyingchi(t)

    圖2及圖3結(jié)果表明,2000—2015年來(lái)自農(nóng)村生活污水、農(nóng)田化肥流失、畜禽養(yǎng)殖的TN、TP年輸出負(fù)荷與研究區(qū)TN、TP年輸出負(fù)荷空間變化規(guī)律十分相似,除了墨脫縣、察隅縣、朗縣,其他行政區(qū)內(nèi)三類營(yíng)養(yǎng)源的貢獻(xiàn)量與TN、TP的年輸出總負(fù)荷均呈逐年增加的趨勢(shì)。2000—2015年三類營(yíng)養(yǎng)源的TN輸出及空間分布組成也存在差異,林芝縣、工布江達(dá)縣及朗縣TN負(fù)荷主要來(lái)源為畜禽養(yǎng)殖;米林縣TN負(fù)荷最初主要來(lái)源為農(nóng)田化肥,而到2015年主要污染源變?yōu)樾笄蒺B(yǎng)殖;其他3個(gè)縣TN負(fù)荷主要來(lái)自農(nóng)田化肥;由圖2結(jié)果得出,畜禽養(yǎng)殖對(duì)米林縣的TN貢獻(xiàn)率變化最為明顯,其他行政區(qū)內(nèi)各污染源貢獻(xiàn)率雖有變化但不明顯。2000—2015年三類營(yíng)養(yǎng)源的TP輸出及其空間分布組成與研究區(qū)TN的規(guī)律相似。各行政區(qū)內(nèi)不同污染源的年貢獻(xiàn)量及TN、TP的年負(fù)荷雖不盡相同,但存在共同規(guī)律.各行政區(qū)TN、TP的主要來(lái)源為畜禽養(yǎng)殖和農(nóng)田化肥,該兩類污染源的TN、TP貢獻(xiàn)量之和在各行政區(qū)都超過總量的90%。研究區(qū)內(nèi)不同污染源的空間分布及貢獻(xiàn)量隨時(shí)間的推移而發(fā)生變化,圖2及圖3結(jié)果顯示研究區(qū)內(nèi)墨脫縣、朗縣及察隅縣TN、TP的年輸出負(fù)荷沒有隨時(shí)間推移而呈上升趨勢(shì),這與該3個(gè)行政區(qū)的交通情況較差、農(nóng)村人口基數(shù)少、土地利用率低等條件有直接關(guān)系;米林縣、工布江達(dá)縣及林芝縣TN、TP的年輸出負(fù)荷隨時(shí)間推移有明顯的上升趨勢(shì),該3個(gè)行政區(qū)地理位置優(yōu)越,交通便利、社會(huì)經(jīng)濟(jì)發(fā)展迅速,從而導(dǎo)致TN、TP的年輸出負(fù)荷增加。

    3.3 TN、TP年輸出負(fù)荷強(qiáng)度空間分布

    面源污染負(fù)荷的產(chǎn)排具有很強(qiáng)的空間性,根據(jù)研究區(qū)7個(gè)縣2015年TN、TP年輸出負(fù)荷,借助GIS空間分析平臺(tái)采用非監(jiān)督分類、顏色渲染等手段分析各行政區(qū)平均負(fù)荷強(qiáng)度分布專題圖(圖4),更直觀地表達(dá)TN、TP的空間產(chǎn)排特性。圖4結(jié)果顯示,米林縣的TN、TP負(fù)荷強(qiáng)度最高,遠(yuǎn)遠(yuǎn)高出全區(qū)平均水平,成為林芝地區(qū)的面源污染重點(diǎn)治理行政區(qū);察隅縣與墨脫縣TN、TP負(fù)荷強(qiáng)度遠(yuǎn)低于全區(qū)水平,且兩者TN、TP總負(fù)荷量一直保持著較低的增速,為區(qū)內(nèi)環(huán)境最好地域;林芝縣及工布江達(dá)縣TN、TP負(fù)荷強(qiáng)度雖低于米林縣,但高于全區(qū)平均水平,考慮其TN、TP年負(fù)荷變化趨勢(shì),兩者污染相對(duì)嚴(yán)重。為此,應(yīng)該對(duì)米林縣、林芝縣及工布江達(dá)縣內(nèi)的畜禽養(yǎng)殖、農(nóng)田化肥污染源采取必要的預(yù)防控制措施。

    4 討論

    圖2 林芝地區(qū)2000—2015年不同污染源TN負(fù)荷空間分布Figure 2 Loads of TN from various source pollution in 2000,2005,2010 and 2015 in Nyingchi

    圖3 林芝地區(qū)2000—2015年不同污染源TP負(fù)荷空間分布Figure 3 Loads of TP from various pollution source in 2000,2005,2010 and 2015 in Nyingchi

    由于西藏高原生態(tài)區(qū)沒有農(nóng)業(yè)面源污染的研究區(qū)域可供參考,本次農(nóng)業(yè)面源污染輸出負(fù)荷估算中,是根據(jù)《第一次全國(guó)污染源普查城鎮(zhèn)生活源產(chǎn)排污系數(shù)手冊(cè)》,參考《第一次全國(guó)污染源普查農(nóng)業(yè)污染源肥料流失系數(shù)手冊(cè)》的研究成果,以及農(nóng)業(yè)面源污染研究的文獻(xiàn)資料來(lái)選取該區(qū)域的輸出系數(shù),結(jié)果雖無(wú)法得到有效驗(yàn)證,但成果與實(shí)地踏勘結(jié)果基本相符。將本研究結(jié)果與三峽庫(kù)區(qū)及哈爾濱地區(qū)的農(nóng)業(yè)面源污染研究成果作以比較.農(nóng)業(yè)面源污染源研究對(duì)象同為生活污水、農(nóng)田種植、畜禽養(yǎng)殖;林芝地區(qū)污染物TN和TP的單位面積年輸出負(fù)荷雖低于三峽庫(kù)區(qū)和哈爾濱地區(qū),但年輸出負(fù)荷總量及空間分布的特征和影響因素與三峽庫(kù)區(qū)及哈爾濱地區(qū)的規(guī)律相同,研究區(qū)污染物TN和TP的年輸出負(fù)荷均隨著農(nóng)村人口、畜禽養(yǎng)殖量及農(nóng)田種植面積的增加呈逐年上升的趨勢(shì)。

    圖4 林芝地區(qū)2015年TN、TP年負(fù)荷空間強(qiáng)度分布圖Figure 4 Spatial distribution of TN and TP loads in Nyingchi in 2015

    本次研究污染指標(biāo)只考慮了TN和TP,污染源只涉及了農(nóng)村生活、農(nóng)田化肥、畜禽養(yǎng)殖三類污染源,未考慮農(nóng)藥、農(nóng)膜和農(nóng)村生活垃圾等污染源以及NH3-N、COD、BOD5等污染指標(biāo)。研究區(qū)內(nèi)水資源較為豐富,河網(wǎng)密集,但缺乏水質(zhì)監(jiān)測(cè)分析數(shù)據(jù),未能將區(qū)內(nèi)農(nóng)村水質(zhì)監(jiān)測(cè)評(píng)價(jià)與農(nóng)業(yè)面源污染排放進(jìn)行響應(yīng)分析。因此在今后的研究中,應(yīng)將區(qū)內(nèi)農(nóng)業(yè)面源污染時(shí)空特征研究與區(qū)內(nèi)社會(huì)經(jīng)濟(jì)、人類活動(dòng)特征以及水質(zhì)監(jiān)測(cè)分析結(jié)果等因素結(jié)合起來(lái),建立起基于時(shí)空尺度的林芝市農(nóng)業(yè)面源污染與農(nóng)村水質(zhì)監(jiān)測(cè)的動(dòng)態(tài)響應(yīng)系統(tǒng),從而深入分析區(qū)內(nèi)各類污染源時(shí)空特征,完善西藏生態(tài)區(qū)農(nóng)業(yè)面源污染問題的研究體系,為西藏生態(tài)安全屏障建設(shè)提供一定的參考。

    5 結(jié)論

    本文利用輸出系數(shù)模型,對(duì)林芝地區(qū)2000—2015年農(nóng)業(yè)面源污染TN、TP年輸出負(fù)荷進(jìn)行了估算和空間分析,結(jié)果表明.

    (1)林芝地區(qū)農(nóng)業(yè)面源污染TN、TP年輸出總負(fù)荷呈逐年上升的趨勢(shì),尤其是米林縣、林芝縣及工布江達(dá)縣TN、TP年輸出總負(fù)荷增速尤為明顯。

    (2)畜禽養(yǎng)殖和農(nóng)田化肥是林芝地區(qū)農(nóng)業(yè)面源污染最大的影響因素,是TN、TP最主要的來(lái)源。

    (3)林芝地區(qū)農(nóng)業(yè)面源污染較為嚴(yán)重的行政區(qū)為米林縣、林芝縣和工布江達(dá)縣,是今后的重點(diǎn)防治區(qū)域。

    在僅為2.38萬(wàn)hm2的耕地與居住面積內(nèi),林芝地區(qū)2015年污染物TN、TP的年輸出負(fù)荷分別為4 924.67、791.59 t,且區(qū)內(nèi)農(nóng)業(yè)面源污染年輸出負(fù)荷呈穩(wěn)步上升的趨勢(shì),農(nóng)業(yè)面源污染形勢(shì)日益嚴(yán)峻,故開展對(duì)該地區(qū)的農(nóng)業(yè)面源污染防控工作已刻不容緩。應(yīng)合理規(guī)劃畜禽養(yǎng)殖區(qū)域和規(guī)模,優(yōu)化農(nóng)業(yè)種植結(jié)構(gòu),科學(xué)施肥,減少化肥的流失,健全農(nóng)村污染處理措施,提高農(nóng)民環(huán)保意識(shí),從而有效地改善林芝地區(qū)農(nóng)業(yè)面源污染現(xiàn)狀,保護(hù)雅江流域的生態(tài)環(huán)境。

    [1]匡啟義.舒城縣農(nóng)業(yè)面源污染現(xiàn)狀及治理對(duì)策[J].現(xiàn)代農(nóng)業(yè)科技, 2015(9).226-227.

    KUANG Qi-yi.Status quo of agricultural non-point pollution in Shucheng County and its control strategies[J].Modern Agricultural Sciences and Technology,2015(9).226-227.

    [2]劉亞琦,劉加珍.聊城市農(nóng)業(yè)面源污染的問題與治理措施分析[J].河北農(nóng)業(yè)科學(xué),2015,19(2).82-86.

    LIU Ya-qi,LIU Jia-zhen.Analysis on the problems and countermeasures of agricultural non-point source pollutions in Liaocheng City[J].Journal of Hebei Agricultural Sciences,2015,19(2).82-86.

    [3]張維理,冀宏杰,Kolbe H,等.中國(guó)農(nóng)業(yè)面源污染形勢(shì)估計(jì)及控制對(duì)策Ⅱ.歐美國(guó)家農(nóng)業(yè)面源污染狀況及控制[J].中國(guó)農(nóng)業(yè)科學(xué),2004, 37(7).1018-1025.

    ZHANG Wei-li,JI Hong-jie,Kolbe H,et al.Estimation of agricultural non-point source pollution in China and the alleviating strategiesⅡ. Status of agricultural non-point source pollution and the alleviating strategies in European and American countries[J].Scientia Agricultura Sinica,2004,37(7).1018-1025.

    [4]安娜.遼寧省農(nóng)業(yè)面源污染現(xiàn)狀及管控措施[J].現(xiàn)代農(nóng)業(yè)科技, 2015(21).212-213.

    AN Na.Current situation and management measures for agricultural non-point source pollution in Liaoning Province[J].Modern Agricultural Sciences and Technology,2015(21).212-213.

    [5]宮曉敏.農(nóng)業(yè)面源污染對(duì)濱岸帶微生物群落的影響研究[D].山東農(nóng)業(yè)大學(xué),2013.

    GONG Xiao-min.Study on impacts of agricultural non-point source pollution on microbial community in riparian zone[D].Shandong Agricultural University,2013.

    [6]林雪原,荊延德,孫笑笑,等.濰坊市農(nóng)業(yè)面源污染評(píng)價(jià)及控制對(duì)策[J].曲阜師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,41(4).84-88.

    LIN Xue-yuan,JING Yan-de,SUN Xiao-xiao,et al.Evaluation and control strategies of agricultural non-point source pollution in Weifang City[J].Journal of Qufu Normal University,2015,41(4).84-88.

    [7]應(yīng)蘭蘭,侯西勇,路曉,等.我國(guó)非點(diǎn)源污染研究中輸出系數(shù)問題[J].水資源與水工程學(xué)報(bào),2010,21(6).90-95.

    YING Lan-lan,HOU Xi-yong,LU Xiao,et al.Discussion on the export coefficient method in non-point source pollution studies in China[J].Journal of Water Resources&Water Engineering,2010,21(6).90-95.

    [8]茆峰,蘇饋?zhàn)?康加廷,等.基于改進(jìn)輸出系數(shù)法的礦區(qū)重金屬面源污染負(fù)荷核算模型[J].環(huán)境科學(xué)研究,2012,25(2).207-211.

    MAO Feng,SU Kui-zu,KANG Jia-ting,et al.Research on heavy metal load accounting in mining areas using improved export coefficient method[J].Research of Environmental Sciences,2012,25(2).207-211.

    [9]姚錫良.農(nóng)村非點(diǎn)源污染負(fù)荷核算研究[D].華南理工大學(xué),2012.

    YAO Xi-liang.Study of the rural non-point source pollution load calculation[D].South China University of Technology,2012.

    [10]胡蕓蕓,王永東,李廷軒,等.沱江流域農(nóng)業(yè)面源污染排放特征解析[J].中國(guó)農(nóng)業(yè)科學(xué),2015,48(18).3654-3665.

    HU Yun-yun,WANG Yong-dong,LI Ting-xuan,et al.Characteristics analysis of agricultural non-point source pollution on Tuojiang River Basin[J].ScientiaAgriculturaSinica,2015,48(18).3654-3665.

    [11]朱立安,王繼增,胡耀國(guó),等.畜禽養(yǎng)殖非點(diǎn)源污染及其生態(tài)控制[J].水土保持通報(bào),2005,25(2).40-43.

    ZHU Li-an,WANG Ji-zeng,HU Yao-guo,et al.Non-point sources pollution from livestock and poultry and ecological control[J].Bulletin of Soil and Water Conservation,2005,25(2).40-43.

    [12]陳媛,郭秀銳,程水源,等.SWAT模型在三峽庫(kù)區(qū)流域非點(diǎn)源污染模擬的適用性研究[J].安全與環(huán)境學(xué)報(bào),2012,12(2).146-152.

    CHEN Yuan,GUO Xiu-rui,CHENG Shui-yuan,et al.On the applicability of SWAT model to the non-point source pollution in the watershed of the Three-Gorge Reservoir[J].Journal of Safety&Environment,2012,12(2).146-152.

    [13]王和意,劉敏,劉巧梅,等.城市降雨徑流非點(diǎn)源污染分析與研究進(jìn)展[J].城市環(huán)境與城市生態(tài),2003,16(6).283-285.

    WANG He-yi,LIU Min,LIU Qiao-mei,et al.Analysis and research progress on urban rainfall runoff non-point source pollution[J].Urban Environment&Urban Ecology,2003,16(6).283-285.

    [14]田甜,劉瑞民,王秀娟,等.三峽庫(kù)區(qū)大寧河流域非點(diǎn)源污染輸出風(fēng)險(xiǎn)分析[J].環(huán)境科學(xué)與技術(shù),2011,34(6).185-190.

    TIAN Tian,LIU Rui-min,WANG Xiu-juan,et al.Risk assessment of nitrogen and phosphorus export in Daning River of Three Gorges Reservoir Area[J].Environmental Science&Technology,2011,34(6). 185-190.

    [15]劉瑞民,何孟常,王秀娟.大遼河流域上游非點(diǎn)源污染輸出風(fēng)險(xiǎn)分析[J].環(huán)境科學(xué),2009,30(3).663-667.

    LIU Rui-min,HE Meng-chang,WANG Xiu-juan.Risk assessment of nitrogen and phosphorus export in upper reach of Daliao River Watershed[J].Environmental Science,2009,30(3).663-667.

    [16]馬亞麗,敖天其,張洪波,等.基于輸出系數(shù)模型瀨溪河流域?yàn)o縣段面源分析[J].四川農(nóng)業(yè)大學(xué)學(xué)報(bào),2013,31(1).53-59.

    MA Ya-li,AO Tian-qi,ZHANG Hong-bo,et al.Non-point source analysis of Laixi River Basin in Luxian County based on output coefficient model[J].Journal of Sichuan Agricultural University,2013,31(1). 53-59.

    [17]劉瑞民,沈珍瑤,丁曉雯,等.應(yīng)用輸出系數(shù)模型估算長(zhǎng)江上游非點(diǎn)源污染負(fù)荷[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2008,27(2).677-682.

    LIU Rui-min,SHEN Zhen-yao,DING Xiao-wen,et al.Application of export coefficient model in simulating pollution load of non-point source in upper reach of Yangtze River Basin[J].Journal of Agro-Environmental Science,2008,27(2).677-682.

    [18]陸建忠,陳曉玲,肖靖靖,等.改進(jìn)的輸出系數(shù)法在農(nóng)業(yè)污染源估算中的應(yīng)用[J].華中師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2012,46(3).373-378.

    LU Jian-zhong,CHEN Xiao-ling,XIAO Jing-jing,et al.Application of improved export coefficient model in estimation of agricultural source pollutant[J].Journal of Huazhong Normal University(Nat Sci),2012, 46(3).373-378.

    [19]張立坤,香寶,胡鈺,等.基于輸出系數(shù)模型的呼蘭河流域非點(diǎn)源污染輸出風(fēng)險(xiǎn)分析[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2014,33(1).148-154.

    ZHANG Li-kun,XIANG Bao,HU Yu,et al.Risk assessment of nonpoint source pollution in Hulan River Basin using an output coefficient model[J].Journal of Agro-Environment Science,2014,33(1).148-154.

    [20]Johnes P J,Heathwaite A L.Modeling the impact of land use change on water quality in agricultural catchments[J].Hydrological Processes, 1997,11(3).269-286.

    [21]蔡明,李懷恩,莊詠濤,等.改進(jìn)的輸出系數(shù)法在流域非點(diǎn)源污染負(fù)荷估算中的應(yīng)用[J].水利學(xué)報(bào),2004,35(7).40-45.

    CAI Ming,LI Huai-en,ZHUANG Yong-Tao,et al.Application of modified export coefficient method in polluting load estimation of nonpoint source pollution[J].Journal of Hydraulic Engineering,2004,35 (7).40-45.

    [22]李懷恩,莊詠濤.預(yù)測(cè)非點(diǎn)源營(yíng)養(yǎng)負(fù)荷的輸出系數(shù)法研究進(jìn)展與應(yīng)用[J].西安理工大學(xué)學(xué)報(bào),2003,19(4).307-312.

    LI Huai-en,ZHUANG Yong-tao.The export coefficient modeling approach for load prediction of nutrient from nonpoint source and its application[J].Journal of Xi′an University of Technology,2003,19(4). 307-312.

    Spatio-temporal changes and distribution characteristics of pollutant loads from agricultural non-point sources in Nyingchi,Tibet

    HAO Shou-ning1,FU Yi-cheng2*
    (1.XiZang Agriculture and Animal Husbandry College,Nyingchi 860000,China;2.China Institute of Water Resources and Hydropower Research,Beijing 100038,China)

    This paper is aimed at analyzing the trends and characteristics of agricultural non-point source pollution in Nyingchi,Tibet.The annual emission loads of total nitrogen(TN)and total phosphorus(TP)from the agricultural non-point source pollution were estimated by applying the export coefficient model to the statistical data of agriculture in Nyingchi from 2000—2015.The agricultural non-point source pollution can be caused by rural domestic sewage,application of chemical fertilizers,and livestock and poultry breeding.TN and TP showed a steady upward trend year by year.The loads of agricultural non-point pollution were higher in the Gongbujiangda,Bomi,and Nyingchi counties,whereas that in the Milin county was the highest.The pollution load of TN and TP owing to livestock and poultry breeding showed a significant increasing trend with time among the three nutrient sources.Overall,90%of the agricultural non-point source pollution load resulted from livestock and poultry breeding and the application of chemical fertilizers.Importantly,the situation of the agricultural non-point source pollution in Nyingchi is not optimistic,especially in the developed areas of agriculture and animal husbandry.It is necessary to carry out prevention and control measures for improving the situation of agricultural non-point source pollution in this district.

    agricultural non-point source;TN;TP;Nyingchi;distribution characteristics

    X592

    A

    1672-2043(2017)07-1308-08

    10.11654/jaes.2017-0499

    郝守寧,付意成.林芝市農(nóng)業(yè)面源污染負(fù)荷時(shí)空變化與分布特征[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2017,36(7).1308-1315.

    HAO Shou-ning,FU Yi-cheng.Spatio-temporal changes and distribution characteristics of pollutant loads from agricultural non-point sources in Nyingchi, Tibet[J].Journal of Agro-Environment Science,2017,36(7).1308-1315.

    2017-04-06

    郝守寧(1986—),男,河南滑縣人,講師,從事流域水資源管理與保護(hù)研究。E-mail:645427848@qq.com

    西藏農(nóng)牧學(xué)院高層次人才科研啟動(dòng)費(fèi)項(xiàng)目(RC201503);國(guó)家自然科學(xué)基金項(xiàng)目(51409269);西藏自治區(qū)科技廳重點(diǎn)科研項(xiàng)目(2015XZ01G34);西藏自治區(qū)科技廳地區(qū)基金項(xiàng)目(2016ZR-15-3)

    Project supported:The Special Foundation for Young Scientists of XiZang Agriculture and Animal Husbandry College(RC201503);The National Natural Science Foundation of China(51409269);The Key Program Scientific of the Science and Technology of Tibet,China(2015XZ01G34); The Scientific Research Fundation of the Science and Technology of Tibet,China(2016ZR-15-3)

    猜你喜歡
    林芝地區(qū)面源污染源
    農(nóng)業(yè)面源污染的危害與治理
    澄江市農(nóng)業(yè)面源污染成因及對(duì)策
    持續(xù)推進(jìn)固定污染源排污許可管理全覆蓋
    基于污染源解析的空氣污染治理對(duì)策研究
    十二五”期間佳木斯市污染源排放狀況分析
    西藏林芝地區(qū)零度層亮帶識(shí)別
    西藏一房屋建筑工程混凝土冬季施工案例分析
    基于SWAT模型的漳河流域面源污染模擬研究
    看不見的污染源——臭氧
    農(nóng)業(yè)面源污染對(duì)水質(zhì)的影響及防治對(duì)策
    岛国在线免费视频观看| 国产私拍福利视频在线观看| 亚洲欧美日韩东京热| 亚洲欧美精品综合久久99| 男插女下体视频免费在线播放| 亚洲最大成人av| 精品人妻偷拍中文字幕| 国产免费男女视频| 蜜桃久久精品国产亚洲av| 女人十人毛片免费观看3o分钟| 99久久中文字幕三级久久日本| 国产精品福利在线免费观看| 日本与韩国留学比较| 午夜影院日韩av| 日日啪夜夜撸| 亚洲自拍偷在线| 在线观看av片永久免费下载| 欧美三级亚洲精品| 亚洲中文字幕日韩| 少妇熟女aⅴ在线视频| 亚洲无线观看免费| 欧美潮喷喷水| 日日夜夜操网爽| 特级一级黄色大片| 欧美日本视频| 国内精品久久久久精免费| 亚洲在线自拍视频| 亚洲天堂国产精品一区在线| 干丝袜人妻中文字幕| 亚洲va在线va天堂va国产| 两个人视频免费观看高清| 偷拍熟女少妇极品色| 国产91精品成人一区二区三区| 国产av在哪里看| 国语自产精品视频在线第100页| 亚洲第一区二区三区不卡| 国产欧美日韩精品一区二区| 欧美高清性xxxxhd video| 1024手机看黄色片| 露出奶头的视频| 国产精品久久久久久精品电影| 两个人的视频大全免费| 99久久九九国产精品国产免费| 国产伦精品一区二区三区四那| 中文字幕免费在线视频6| 成人av在线播放网站| 久久久久久久精品吃奶| 久久人人爽人人爽人人片va| 97人妻精品一区二区三区麻豆| 一进一出抽搐gif免费好疼| 国产伦在线观看视频一区| 又爽又黄a免费视频| 极品教师在线视频| 熟女人妻精品中文字幕| 中亚洲国语对白在线视频| 国产成人福利小说| 97超视频在线观看视频| 精品久久久久久久久亚洲 | 国产精品国产三级国产av玫瑰| 亚洲精品成人久久久久久| 全区人妻精品视频| 国产亚洲精品av在线| 国产伦精品一区二区三区四那| 九九爱精品视频在线观看| 一级av片app| 熟女人妻精品中文字幕| 久久久久九九精品影院| 亚洲精华国产精华精| 日本免费a在线| 久久精品综合一区二区三区| 99久久精品热视频| 精品国产三级普通话版| 99久久精品热视频| 神马国产精品三级电影在线观看| 国产一级毛片七仙女欲春2| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品乱码久久久v下载方式| 成人性生交大片免费视频hd| 亚洲五月天丁香| 听说在线观看完整版免费高清| 日韩中字成人| 男人舔女人下体高潮全视频| 婷婷色综合大香蕉| 欧美又色又爽又黄视频| 欧美日韩综合久久久久久 | 又紧又爽又黄一区二区| 亚洲经典国产精华液单| 欧美日韩综合久久久久久 | 欧美绝顶高潮抽搐喷水| 如何舔出高潮| 国产精品美女特级片免费视频播放器| 丝袜美腿在线中文| 国国产精品蜜臀av免费| 色吧在线观看| 国产一区二区三区av在线 | 免费搜索国产男女视频| av在线天堂中文字幕| 天堂√8在线中文| 久久国产乱子免费精品| 永久网站在线| 一级黄色大片毛片| 特级一级黄色大片| 中国美白少妇内射xxxbb| 免费一级毛片在线播放高清视频| 最近最新中文字幕大全电影3| 老司机午夜福利在线观看视频| 别揉我奶头~嗯~啊~动态视频| 51国产日韩欧美| 能在线免费观看的黄片| 免费观看的影片在线观看| 亚洲精华国产精华精| 少妇猛男粗大的猛烈进出视频 | 亚洲熟妇熟女久久| 亚洲内射少妇av| 婷婷六月久久综合丁香| 97超视频在线观看视频| 久久99热这里只有精品18| 国内少妇人妻偷人精品xxx网站| av在线亚洲专区| 窝窝影院91人妻| 国产高清三级在线| 亚洲va日本ⅴa欧美va伊人久久| 精品99又大又爽又粗少妇毛片 | 九九在线视频观看精品| av天堂中文字幕网| 在线观看免费视频日本深夜| 欧美精品啪啪一区二区三区| 看免费成人av毛片| 亚洲成av人片在线播放无| 午夜视频国产福利| 天天一区二区日本电影三级| 在线播放无遮挡| 欧美日韩瑟瑟在线播放| 国产精品亚洲美女久久久| 久久久久久久久久成人| 人人妻人人看人人澡| 久久人妻av系列| 校园春色视频在线观看| 亚洲av不卡在线观看| 中文在线观看免费www的网站| 日韩欧美国产在线观看| 伦精品一区二区三区| 亚洲电影在线观看av| 日韩一本色道免费dvd| 91av网一区二区| 在线免费观看的www视频| 欧美性猛交黑人性爽| 成人午夜高清在线视频| 成人无遮挡网站| 国产真实乱freesex| 国产精品伦人一区二区| 黄色欧美视频在线观看| 在线国产一区二区在线| 变态另类成人亚洲欧美熟女| 亚洲第一区二区三区不卡| 欧美又色又爽又黄视频| 婷婷六月久久综合丁香| 亚洲av日韩精品久久久久久密| 精品乱码久久久久久99久播| 国产国拍精品亚洲av在线观看| 国产一区二区激情短视频| 欧美日韩中文字幕国产精品一区二区三区| 床上黄色一级片| 长腿黑丝高跟| 国国产精品蜜臀av免费| 亚洲av日韩精品久久久久久密| 国产av不卡久久| 亚洲av免费在线观看| 十八禁网站免费在线| 亚洲精品粉嫩美女一区| 天美传媒精品一区二区| 亚洲精品影视一区二区三区av| 精品久久久久久久末码| 69av精品久久久久久| 美女cb高潮喷水在线观看| 麻豆成人午夜福利视频| 在线看三级毛片| 精品久久久久久成人av| 老女人水多毛片| 精品一区二区三区av网在线观看| 国产精品永久免费网站| 精品人妻一区二区三区麻豆 | 国产精品av视频在线免费观看| eeuss影院久久| 一区二区三区高清视频在线| 特级一级黄色大片| 国产大屁股一区二区在线视频| 最新在线观看一区二区三区| 夜夜夜夜夜久久久久| 欧美国产日韩亚洲一区| 中文字幕高清在线视频| 免费无遮挡裸体视频| 动漫黄色视频在线观看| 村上凉子中文字幕在线| 免费av毛片视频| 最近视频中文字幕2019在线8| 在线观看午夜福利视频| 别揉我奶头~嗯~啊~动态视频| 乱系列少妇在线播放| 欧美丝袜亚洲另类 | 色综合色国产| 亚洲va日本ⅴa欧美va伊人久久| 波野结衣二区三区在线| 永久网站在线| 精品人妻1区二区| 在线免费观看的www视频| 久久精品影院6| 嫩草影院精品99| 欧美成人性av电影在线观看| 黄色视频,在线免费观看| 男女之事视频高清在线观看| 最近最新免费中文字幕在线| 亚洲va在线va天堂va国产| 可以在线观看的亚洲视频| 国内精品久久久久久久电影| 中出人妻视频一区二区| 国产又黄又爽又无遮挡在线| 啦啦啦韩国在线观看视频| 中文资源天堂在线| 黄色一级大片看看| 熟女电影av网| 此物有八面人人有两片| 精华霜和精华液先用哪个| 美女黄网站色视频| 亚洲18禁久久av| 女生性感内裤真人,穿戴方法视频| 69人妻影院| 能在线免费观看的黄片| 亚洲人成网站高清观看| 成年版毛片免费区| 国产三级中文精品| 久久久久久伊人网av| 99九九线精品视频在线观看视频| 两个人视频免费观看高清| 色精品久久人妻99蜜桃| 麻豆精品久久久久久蜜桃| 精品久久久噜噜| 国产精品乱码一区二三区的特点| 小说图片视频综合网站| 成人毛片a级毛片在线播放| eeuss影院久久| АⅤ资源中文在线天堂| 亚洲性夜色夜夜综合| av视频在线观看入口| 成人无遮挡网站| 欧美xxxx性猛交bbbb| 欧美另类亚洲清纯唯美| 久久精品国产亚洲av天美| 色综合色国产| 国产精品一区二区三区四区免费观看 | 亚洲精品成人久久久久久| АⅤ资源中文在线天堂| 男女视频在线观看网站免费| 女人被狂操c到高潮| 久久久久久久久中文| 动漫黄色视频在线观看| 午夜福利欧美成人| 偷拍熟女少妇极品色| 18禁黄网站禁片午夜丰满| 亚洲天堂国产精品一区在线| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利在线观看免费完整高清在 | 国产精品久久久久久亚洲av鲁大| 久久精品夜夜夜夜夜久久蜜豆| 床上黄色一级片| 欧美成人a在线观看| 欧美3d第一页| 日韩大尺度精品在线看网址| 人人妻人人看人人澡| 黄色一级大片看看| 国产极品精品免费视频能看的| 亚洲av中文av极速乱 | 国内精品久久久久久久电影| 99热精品在线国产| 国产精品亚洲一级av第二区| 国产不卡一卡二| 国产黄色小视频在线观看| 欧美色欧美亚洲另类二区| av专区在线播放| 国产免费一级a男人的天堂| 亚州av有码| 中亚洲国语对白在线视频| 欧美色欧美亚洲另类二区| 高清在线国产一区| 成人永久免费在线观看视频| 麻豆久久精品国产亚洲av| 少妇被粗大猛烈的视频| 中文字幕熟女人妻在线| 一本一本综合久久| 亚洲图色成人| 国产伦在线观看视频一区| 欧美一区二区精品小视频在线| 丰满的人妻完整版| 国产伦人伦偷精品视频| 简卡轻食公司| 中文字幕人妻熟人妻熟丝袜美| 赤兔流量卡办理| 欧美日韩黄片免| 久久这里只有精品中国| 天美传媒精品一区二区| 亚洲av免费高清在线观看| 内地一区二区视频在线| 尾随美女入室| 老熟妇仑乱视频hdxx| 亚洲美女视频黄频| 国产精品精品国产色婷婷| 国产欧美日韩一区二区精品| 在线播放国产精品三级| 免费看光身美女| 日本五十路高清| 亚洲av成人av| 好男人在线观看高清免费视频| 久久精品国产鲁丝片午夜精品 | 亚洲av.av天堂| 最近中文字幕高清免费大全6 | 国产aⅴ精品一区二区三区波| 内射极品少妇av片p| 亚洲精品在线观看二区| 国产国拍精品亚洲av在线观看| 97热精品久久久久久| 91精品国产九色| 国产免费男女视频| 亚洲自拍偷在线| 中文字幕熟女人妻在线| 国产成人aa在线观看| 亚洲三级黄色毛片| 麻豆久久精品国产亚洲av| 久久久久久久久久黄片| 高清毛片免费观看视频网站| 内地一区二区视频在线| 婷婷亚洲欧美| avwww免费| 亚洲一区高清亚洲精品| 国产探花在线观看一区二区| 全区人妻精品视频| 九九久久精品国产亚洲av麻豆| 国产高清有码在线观看视频| 久久久国产成人精品二区| 色综合婷婷激情| 22中文网久久字幕| 男女边吃奶边做爰视频| 熟女人妻精品中文字幕| 天堂av国产一区二区熟女人妻| 我的老师免费观看完整版| 亚洲av美国av| 午夜福利在线观看吧| 国产久久久一区二区三区| 一级黄色大片毛片| 一进一出好大好爽视频| 嫩草影院入口| 国产精品久久电影中文字幕| 国产精品98久久久久久宅男小说| 日韩欧美精品v在线| 亚洲性夜色夜夜综合| 欧美色欧美亚洲另类二区| 免费av观看视频| 人人妻,人人澡人人爽秒播| 琪琪午夜伦伦电影理论片6080| 能在线免费观看的黄片| 精品久久久久久久人妻蜜臀av| 精品人妻偷拍中文字幕| 丰满乱子伦码专区| 一级黄片播放器| 韩国av一区二区三区四区| 亚洲精品在线观看二区| 我要看日韩黄色一级片| 欧美中文日本在线观看视频| 尾随美女入室| 国产一区二区三区av在线 | av黄色大香蕉| 国产视频一区二区在线看| 免费电影在线观看免费观看| 久久国内精品自在自线图片| 乱码一卡2卡4卡精品| 动漫黄色视频在线观看| 日韩欧美三级三区| 久久九九热精品免费| 国产综合懂色| 国产极品精品免费视频能看的| 美女黄网站色视频| 国产欧美日韩一区二区精品| 日日摸夜夜添夜夜添小说| 成年女人看的毛片在线观看| 99热只有精品国产| 99久久精品国产国产毛片| 啦啦啦观看免费观看视频高清| 两个人的视频大全免费| 麻豆国产97在线/欧美| 久9热在线精品视频| 麻豆一二三区av精品| 在线观看av片永久免费下载| 亚洲精品国产成人久久av| x7x7x7水蜜桃| 中文字幕高清在线视频| av在线观看视频网站免费| 日本a在线网址| 精品99又大又爽又粗少妇毛片 | 欧美成人a在线观看| 亚洲图色成人| 欧美一区二区精品小视频在线| 日韩欧美免费精品| 亚洲avbb在线观看| 中文亚洲av片在线观看爽| 欧美黑人欧美精品刺激| 简卡轻食公司| 欧美绝顶高潮抽搐喷水| 一区福利在线观看| 亚洲国产色片| 亚洲av中文字字幕乱码综合| 亚洲成人中文字幕在线播放| 久久午夜福利片| 久久人妻av系列| 人人妻,人人澡人人爽秒播| 国产高清有码在线观看视频| 高清在线国产一区| 国产午夜福利久久久久久| 久久人人爽人人爽人人片va| 99热6这里只有精品| 一个人看视频在线观看www免费| 亚洲精品日韩av片在线观看| av视频在线观看入口| 免费看日本二区| 熟女电影av网| 国产一区二区三区在线臀色熟女| xxxwww97欧美| 欧美日韩国产亚洲二区| 少妇丰满av| 男人和女人高潮做爰伦理| 日日摸夜夜添夜夜添小说| 自拍偷自拍亚洲精品老妇| 女人十人毛片免费观看3o分钟| 又爽又黄a免费视频| 久久久久久九九精品二区国产| 俄罗斯特黄特色一大片| 欧美成人免费av一区二区三区| av.在线天堂| 我要看日韩黄色一级片| 精品久久久久久成人av| 麻豆久久精品国产亚洲av| 欧美日韩乱码在线| 老女人水多毛片| 亚洲自拍偷在线| 亚洲国产色片| 欧美区成人在线视频| 久久久久久久久久黄片| 美女 人体艺术 gogo| 在线观看66精品国产| АⅤ资源中文在线天堂| 婷婷丁香在线五月| 我的老师免费观看完整版| 午夜亚洲福利在线播放| 国产在线精品亚洲第一网站| 2021天堂中文幕一二区在线观| 亚洲五月天丁香| 中文字幕av在线有码专区| 久久亚洲精品不卡| 国内精品宾馆在线| 亚洲无线观看免费| av天堂中文字幕网| 欧美激情久久久久久爽电影| 久久中文看片网| 国产成人av教育| 国产视频一区二区在线看| 嫩草影院精品99| 老女人水多毛片| 亚洲精品一区av在线观看| 午夜老司机福利剧场| 老熟妇乱子伦视频在线观看| 99九九线精品视频在线观看视频| 亚洲不卡免费看| 国内少妇人妻偷人精品xxx网站| 欧美成人性av电影在线观看| 国产色爽女视频免费观看| 国产男靠女视频免费网站| 99riav亚洲国产免费| 中文资源天堂在线| 亚洲av美国av| 精品午夜福利视频在线观看一区| av在线蜜桃| 性插视频无遮挡在线免费观看| 免费一级毛片在线播放高清视频| 久久99热这里只有精品18| 国产视频内射| 免费观看在线日韩| 亚洲黑人精品在线| 看片在线看免费视频| 少妇丰满av| 久久久色成人| 欧美不卡视频在线免费观看| 舔av片在线| 国产成人a区在线观看| 国产免费av片在线观看野外av| 可以在线观看的亚洲视频| 亚洲精品一卡2卡三卡4卡5卡| 色综合站精品国产| 又黄又爽又刺激的免费视频.| 亚洲精品日韩av片在线观看| 在线免费十八禁| 黄色日韩在线| 亚洲狠狠婷婷综合久久图片| 国产探花在线观看一区二区| 国产不卡一卡二| 亚洲人成伊人成综合网2020| 成人特级黄色片久久久久久久| 午夜亚洲福利在线播放| a级一级毛片免费在线观看| 88av欧美| 真实男女啪啪啪动态图| 在线观看一区二区三区| 美女高潮的动态| 91午夜精品亚洲一区二区三区 | 看片在线看免费视频| 免费不卡的大黄色大毛片视频在线观看 | 日韩国内少妇激情av| 噜噜噜噜噜久久久久久91| 亚洲三级黄色毛片| 在线免费观看的www视频| 亚洲avbb在线观看| 人人妻,人人澡人人爽秒播| 别揉我奶头~嗯~啊~动态视频| 在线a可以看的网站| 午夜亚洲福利在线播放| 亚洲午夜理论影院| 床上黄色一级片| 成人美女网站在线观看视频| 色综合亚洲欧美另类图片| 中文字幕免费在线视频6| 亚洲国产精品sss在线观看| 精品无人区乱码1区二区| 国内少妇人妻偷人精品xxx网站| 久久午夜亚洲精品久久| 国产精品1区2区在线观看.| 欧美日韩乱码在线| 亚洲性久久影院| 男女边吃奶边做爰视频| 日韩欧美精品v在线| 真人一进一出gif抽搐免费| 波野结衣二区三区在线| 一夜夜www| 18+在线观看网站| 99在线人妻在线中文字幕| 午夜久久久久精精品| 一a级毛片在线观看| 国产亚洲欧美98| 日韩欧美 国产精品| 久久九九热精品免费| av天堂中文字幕网| 嫩草影院新地址| 尾随美女入室| 欧美精品国产亚洲| 黄色丝袜av网址大全| 午夜福利欧美成人| 三级毛片av免费| 老熟妇仑乱视频hdxx| 亚洲最大成人手机在线| 欧美日韩瑟瑟在线播放| 少妇裸体淫交视频免费看高清| 午夜免费男女啪啪视频观看 | 男女啪啪激烈高潮av片| 国产亚洲av嫩草精品影院| 三级男女做爰猛烈吃奶摸视频| 日韩av在线大香蕉| 老熟妇乱子伦视频在线观看| 一级毛片久久久久久久久女| 亚洲中文日韩欧美视频| 精品久久久久久久久av| av在线老鸭窝| 免费电影在线观看免费观看| 国产大屁股一区二区在线视频| 精品一区二区免费观看| 亚洲人成网站在线播| 精华霜和精华液先用哪个| 色视频www国产| av中文乱码字幕在线| 国产麻豆成人av免费视频| 成人毛片a级毛片在线播放| 乱系列少妇在线播放| 99国产极品粉嫩在线观看| 久久久久性生活片| 亚洲乱码一区二区免费版| 国产精品美女特级片免费视频播放器| 国产高清不卡午夜福利| 亚洲五月天丁香| 一级黄色大片毛片| 天天躁日日操中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 一区二区三区免费毛片| 欧美不卡视频在线免费观看| 嫩草影院入口| 亚洲人成网站在线播放欧美日韩| 夜夜爽天天搞| 成年版毛片免费区| 一个人观看的视频www高清免费观看| 网址你懂的国产日韩在线| 两个人视频免费观看高清| 亚洲无线观看免费| 我要看日韩黄色一级片| 特级一级黄色大片| 观看美女的网站| 看免费成人av毛片| 色播亚洲综合网| 国产精品1区2区在线观看.| 美女cb高潮喷水在线观看| 成人特级黄色片久久久久久久| 99久久精品国产国产毛片| 1000部很黄的大片| 亚洲成人久久爱视频| 婷婷精品国产亚洲av在线| 又爽又黄无遮挡网站| 国产一区二区在线观看日韩| 高清在线国产一区| 美女黄网站色视频| 亚洲三级黄色毛片| 国国产精品蜜臀av免费| 日日啪夜夜撸|