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

    渭河流域(陜西段)面源污染總磷時空特征評估

    2024-01-17 00:00:00馮愛萍王琦羅儀寧張建輝王雪蕾謝成玉黃莉王玉
    關鍵詞:面源污染渭河流域

    關鍵詞:面源污染;DPeRS模型;渭河流域;重點匯水區(qū)

    中圖分類號:X52 文獻標志碼:A 文章編號:1672-2043(2024)12-2776-10 doi:10.11654/jaes.2024-0932

    磷是限制湖泊和河流水體富營養(yǎng)化及藻類生長的關鍵因子,在來自工業(yè)和城鎮(zhèn)生活污水的點源磷污染得到控制后,面源磷污染逐步成為水體污染的主導因子[1]。面源污染是指污染物從非固定的排污口排放,在降雨或融雪過程的作用下,通過徑流的形式進入受納水體,并導致水體富營養(yǎng)化及其他形式的水污染[2]。農(nóng)業(yè)系統(tǒng)是面源污染的主要來源,第一次和第二次全國污染源普查結果顯示,來自農(nóng)業(yè)源水污染物排放量中,總磷分別為28.5萬t和21.20萬t,均占到排放總量的67%以上。因此,分析流域面源污染總磷流失的時空分布特征,有助于控制水體富營養(yǎng)化及水體污染。

    渭河是黃河第一大支流,被譽為陜西省的母親河,對陜西省的經(jīng)濟發(fā)展至關重要。隨著關中城市群經(jīng)濟和城鎮(zhèn)化的快速發(fā)展,從20世紀80年代開始渭河干流陜西段水質(zhì)不斷惡化。曹原等[3]以2014—2023年渭河干流陜西段9 個水質(zhì)監(jiān)測站逐月水質(zhì)監(jiān)測數(shù)據(jù)為依據(jù),選取渭河污染較嚴重的高錳酸鹽指數(shù)、氨氮和總磷共3種污染物指標進行水質(zhì)變化趨勢分析評價,結果表明,2014—2018年渭河干流水質(zhì)很差,沿程水質(zhì)基本為Ⅳ類、Ⅴ類,甚至劣Ⅴ類,且9個水質(zhì)監(jiān)測站中蔡家坡、咸陽公路橋和耿鎮(zhèn)3 個站點2014年總磷濃度均超過地表水Ⅲ類標準限值。隨著渭河流域水污染防治三年行動方案和綜合治理工程的實施,渭河流域點源污染基本得到控制,而面源污染因其來源廣泛、隨機性強、難以監(jiān)測等特點,已成為影響渭河流域水環(huán)境質(zhì)量的主要因素[4]。有研究表明,在各污染源中,對渭河干流陜西段污染貢獻率較多的是生活源和農(nóng)業(yè)面源,工業(yè)源貢獻率較少[5-6]。因此,分析面源污染的分布規(guī)律,識別典型斷面所在匯水區(qū)面源污染主導類型,對于改善渭河流域水環(huán)境質(zhì)量具有重要意義。目前關于渭河流域(陜西段)面源污染的相關研究多集中在水污染特征與水質(zhì)評價[7-8]、流域沿岸農(nóng)業(yè)土壤肥力綜合評價[9]、農(nóng)業(yè)面源污染防治對策[10]、河流生境質(zhì)量分布格局[11]、水生態(tài)安全及健康評價[12-13]等方面,缺少針對該流域面源污染負荷的空間量化及污染類型的深入解析,難以進一步有效支撐渭河流域水污染治理和水環(huán)境持續(xù)改善。

    模型模擬是流域面源污染量化評估的有效手段,相關的模型包括土水評估工具模型(Soil Water Assess?ment Tool,SWAT)[14-15]、農(nóng)業(yè)面源年化模型(Annual?ized Agricultural Nonpoint Source,AnnAGNPS)[16-17]、水文模擬模型(Hydrological Simulation Program-Fortran,HSPF)[18-19]和輸出系數(shù)模型(Export" Coefficient Model,ECM)[20]等,但不同模型的應用條件會受到流域自身特征和輸入數(shù)據(jù)的限制。這些模型中機理模型存在所需數(shù)據(jù)量較大、建模過程較繁雜、模擬所需時間較長等特點;輸出系數(shù)模型則主要依賴于統(tǒng)計數(shù)據(jù)和排污系數(shù)進行估算,這種方法產(chǎn)生的誤差通常較大。而面源污染遙感分布式估算模型(DPeRS模型)[21-22]可實現(xiàn)流域像元尺度多類型面源污染量的空間精細化評估,因此,本研究采用DPeRS模型,對2016—2020年渭河流域(陜西段)農(nóng)田種植源、畜禽養(yǎng)殖源、農(nóng)村生活源、城鎮(zhèn)生活源和水土流失源5個類型的總磷指標進行空間像元尺度評估,深入分析流域面源污染總磷排放負荷和入河負荷的時空分布特征,并篩選重點匯水區(qū)進行主導污染類型識別,研究結果可為管理部門開展渭河流域(陜西段)面源污染及水污染治理提供決策支撐。

    1 材料與方法

    1.1研究區(qū)概況

    本研究選取渭河流域(陜西段)作為研究區(qū)(圖1),對其面源污染時空分布特征進行評估分析。渭河是黃河的最大支流,發(fā)源于甘肅省定西市渭源縣鳥鼠山,主要流經(jīng)今甘肅天水及陜西關中平原的寶雞、咸陽、西安、渭南等地,至渭南市潼關縣匯入黃河,流域總面積為13.48萬km2,其中,陜西省境內(nèi)的流域面積為6.74 萬km2,占渭河流域總面積的50.0%,包括寶雞、咸陽、西安、渭南、銅川、延安、榆林、商洛等8市63個縣(市區(qū))的全部或部分。流域陜西省境內(nèi)土地利用類型以林地、耕地和草地為主,面積占比分別為54.42%、23.14%和13.94%,其中耕地主要分布在關中平原,草地主要分布在流域西北部的榆林和延安地區(qū)。渭河干流全長818 km,陜西境內(nèi)河長502 km,寶雞峽至咸陽為中游,咸陽至潼關為下游;渭河支流眾多,其中南岸的數(shù)量較多,但較大支流集中在北岸,水系呈扇狀分布;集水面積1000km2以上的支流有14條,其中涇河、洛河是渭河的2條最大支流。

    1.2 研究方法

    DPeRS模型是基于二元結構原理構建的遙感分布式面源污染模型[23-24],模型算法以遙感數(shù)據(jù)為驅(qū)動,以遙感像元為模擬單元,耦合定量遙感模型和生態(tài)水文過程模型,可實現(xiàn)流域面源污染負荷的月尺度空間估算和不同污染源在遙感像元上的空間量化表達,重點評估由降水引起的地表徑流型面源污染。DPeRS模型可評估包括農(nóng)田徑流型、農(nóng)村生活型、畜禽養(yǎng)殖型、城鎮(zhèn)生活型和水土流失型5個類型,溶解態(tài)和顆粒態(tài)兩種污染物形態(tài)。溶解態(tài)污染物即表明面源污染物具有水溶性,能隨地表徑流發(fā)生遷移,整個遷移過程受水循環(huán)控制,污染源主要來自農(nóng)田生產(chǎn)過程中肥料、農(nóng)村生活和城鎮(zhèn)生活中非集中處理和排放的生活垃圾和生活污水污染、散養(yǎng)和集中養(yǎng)殖過程未集中處理的糞便污染等;顆粒態(tài)污染物即通過附著于土壤顆粒體而實現(xiàn)遷移運動的污染元素,其遷移過程同水土流失密切相關。本研究采用DPeRS模型核算流域面源污染總磷排放負荷,并進一步采用HEQM模型[25-26]模擬總磷入河系數(shù),進而核算面源污染總磷入河負荷,其中核算總磷入河系數(shù)時采用了研究區(qū)主要水文站的月度觀測數(shù)據(jù)。DPeRS模型算法主要包括溶解態(tài)和顆粒態(tài)面源污染負荷評估方法,模型核心算法具體如下:

    (1)溶解態(tài)面源污染負荷評估方法

    溶解態(tài)面源污染共有4種污染類型,包括農(nóng)田徑流、城鎮(zhèn)生活、農(nóng)村生活和畜禽養(yǎng)殖污染。溶解態(tài)污染物負荷匡算模型統(tǒng)一表達式如下:

    1.3 模型數(shù)據(jù)庫構建

    DPeRS模型運行的數(shù)據(jù)庫主要包括地形數(shù)據(jù)、氣象數(shù)據(jù)、土壤數(shù)據(jù)、土地利用分布數(shù)據(jù)和月度植被覆蓋度數(shù)據(jù)及農(nóng)業(yè)統(tǒng)計數(shù)據(jù)等(表1),所有空間數(shù)據(jù)均統(tǒng)一轉換為30 m空間分辨率的柵格數(shù)據(jù)。其中,數(shù)字高程數(shù)據(jù)(DEM)用于研究區(qū)坡度和坡長數(shù)據(jù)的計算[27];利用薄板樣條滑動平均法進行降水量的空間插值[28];采用最大最小值定量反演算法進行流域植被覆蓋度反演[29];農(nóng)業(yè)統(tǒng)計數(shù)據(jù)為2016—2020年渭河流域(陜西段)縣域尺度數(shù)據(jù),通過輸入輸出法核算農(nóng)田氮磷表觀平衡量[30];基于全省土壤氮磷含量背景數(shù)據(jù),結合實地采樣檢測數(shù)據(jù),采用隨機森林(RF)算法對全省土壤全磷含量數(shù)據(jù)進行空間插值。此外,蘆村河斷面月度流量和水質(zhì)監(jiān)測數(shù)據(jù),為模型結果驗證提供數(shù)據(jù)支撐。

    2 結果與分析

    2.1 面源污染負荷時空分布特征

    基于渭河流域(陜西段)面源污染總磷排放負荷和入河負荷核算結果進行空間統(tǒng)計分析,結果表明:面源污染排放方面,2016—2020年渭河流域(陜西段)面源污染總磷排放負荷和排放量表現(xiàn)為波動上升的趨勢,2020年最高,其平均排放負荷為0.19 t·km-2、排放量為1.30萬t,2018年次之,2017年最低;面源污染入河方面,2016—2020年渭河流域(陜西段)面源污染總磷入河負荷和入河總量整體表現(xiàn)為先上升后下降的趨勢,2018年最為突出,其平均入河負荷為0.04 t·km-2,入河總量為0.26 萬t。2016—2020 年渭河流域(陜西段)面源污染總磷產(chǎn)排量統(tǒng)計信息詳見圖2。

    2016—2020年渭河流域(陜西段)面源污染總磷排放負荷空間分布如圖3所示。5年間排放負荷高值區(qū)主要集中分布在流域西北部的草地分布地帶和中部的耕地種植區(qū),具體在榆林市、延安市北部、渭南市西部及咸陽市中部等區(qū)域。從面源污染總磷排放負荷的時空分布變化上看,2016年關中平原的咸陽市中部和渭南市西部面源污染總磷排放負荷最為突出,這些區(qū)域旱地分布較為集中,排放負荷高與該區(qū)域2016年農(nóng)田磷平衡量相對較高有關;2018年和2020年流域西北部面源污染總磷排放負荷較為突出,其中2020年渭河流域(陜西段)面源污染總磷排放負荷值范圍為0~15.36 t·km-2,平均負荷為0.19 t·km-2,大部分區(qū)域的總磷排放負荷在0~1.00 t·km-2。

    2016—2020年渭河流域(陜西段)面源污染總磷入河負荷空間分布如圖4所示。5年間面源污染總磷入河負荷高值區(qū)主要分布在流域北部的榆林市和延安市西北部,其土地利用類型主要為草地。從時空變化趨勢上看,流域面源污染總磷入河負荷整體上呈現(xiàn)先加重后減輕的趨勢,2018年最為突出,2017年次之,這與流域面源污染總磷入河負荷高值區(qū)的降水量相對較高有關,建議這些區(qū)域進一步加強水土保持工作,以減少草地水土流失引發(fā)的面源污染物排放。

    2.2不同類型面源污染特征

    基于2016—2020年渭河流域(陜西段)面源污染分類核算結果,統(tǒng)計農(nóng)田徑流型、畜禽養(yǎng)殖型、農(nóng)村生活型、城鎮(zhèn)生活型和水土流失型5種類型的面源污染總磷排放量(圖5),結果表明:渭河流域(陜西段)面源污染總磷的首要污染類型為水土流失型,2016—2020年該類型平均排放量占比為67.92%;其次為農(nóng)田徑流型,該類型平均排放量占比為31.77%。從時間變化上看,2016—2020年流域面源污染水土流失型總磷排放量占比整體呈現(xiàn)先升后降的趨勢,其中2020年渭河流域(陜西段)總磷排放總量為14 065 t,水土流失型占比73.43%。因此,應進一步加強區(qū)域水土保持工作,尤其是流域西北部的延安市和榆林市的草地分布地區(qū),以減少水土流失引發(fā)的面源污染物排放,同時也應加強渭河流域中部平原地區(qū)農(nóng)田養(yǎng)分管理以減少養(yǎng)分流失。

    2.3 重點匯水區(qū)主導類型識別及面源污染貢獻率分析

    考慮到地表水監(jiān)測斷面水質(zhì)數(shù)據(jù)的一致性,本研究采用2017—2020年渭河流域(陜西段)范圍內(nèi)地表水監(jiān)測斷面總磷的月度監(jiān)測數(shù)據(jù)進行超標情況評價,以每年度12次監(jiān)測數(shù)據(jù)的算術平均值對地表水監(jiān)測斷面年度超標情況進行評價,并篩選近4年內(nèi)年度超標次數(shù)大于3次的斷面作為重點斷面,結果篩選出岔口、新河入渭和臨河入渭3個重點斷面,以重點斷面所在水文分區(qū)作為重點匯水區(qū)。

    基于2016—2020年渭河流域(陜西段)不同類型面源污染總磷評估結果,統(tǒng)計重點匯水區(qū)內(nèi)不同類型面源污染總磷排放量年均值,進一步對重點匯水區(qū)內(nèi)面源污染主導類型進行解析,結果表明:岔口和新河入渭斷面所在匯水區(qū)的面源污染總磷主導類型為水土流失型(排放量占比分別為61.99%和86.08%),建議此區(qū)域應進一步加強區(qū)域水土保持工作;臨河入渭斷面所在匯水區(qū)的面源污染總磷主導類型為農(nóng)田徑流型(排放量占比58.44%),建議從節(jié)水節(jié)肥等方面提高農(nóng)業(yè)生產(chǎn)效率,在保證土壤養(yǎng)分供給的前提下減少養(yǎng)分流失導致的面源污染排放。重點匯水區(qū)面源污染主導類型空間分布詳見圖6。

    以岔口斷面所在的重點匯水區(qū)(主要地類為林地和旱地)為例,基于耀縣水文站(臨近匯水區(qū)出口)流量月監(jiān)測數(shù)據(jù)、岔口斷面水質(zhì)月度監(jiān)測數(shù)據(jù),結合DPeRS 模型評估的入河量結果,對2017—2019 年岔口斷面所在匯水區(qū)的面源污染貢獻率進行估算,結果表明,該匯水區(qū)總磷面源污染對水質(zhì)總磷指標的貢獻率相對較高,2017—2019 年其各年貢獻率分別為86.53%、78.48% 和67.14%。建議該匯水區(qū)上中游地區(qū)主要關注林地水土保持、下游地區(qū)重點從農(nóng)田養(yǎng)分管理加強面源污染防控。

    2.4 結果驗證

    基于渭河流域內(nèi)現(xiàn)有地面水質(zhì)水文監(jiān)測站分布情況,結合子流域空間區(qū)劃和土地利用空間分布,在充分考慮監(jiān)測數(shù)據(jù)完整性的條件下,本研究選取以面源污染為主的蘆村河斷面作為模型結果驗證點位,該斷面接近三水河入涇河的入口,其所在的小流域土地利用分布以林地、草地和耕地為主。基于蘆村河斷面2020年月度水質(zhì)水量監(jiān)測結果核算斷面年度總磷污染量,與DPeRS 模型核算的該斷面所在小流域面源污染入河量結果進行驗證分析。2020年蘆村河斷面年度總磷污染量為2.89 t,DPeRS模型核算的該斷面所在小流域總磷面源污染入河量為2.61 t,驗證結果表明,假設蘆村河斷面所在小流域范圍內(nèi)污染源全部為面源的前提下,總磷面源污染的模型模擬精度達到90.31%。

    3 討論

    3.1 DPeRS模型的適用性

    對比國內(nèi)外其他面源污染模型,DPeRS模型在模型結構、模擬指標和運行條件等方面具有較大的管理應用優(yōu)勢。該模型以遙感像元為基本模擬單元,與SWAT模型提出的水文響應單元相比,在保證模擬精度的前提下大幅提高了面源污染模擬的空間分辨率;同時DPeRS 模型耦合了定量遙感模型,彌補了無資料地區(qū)或缺資料地區(qū)模型估算的不足,且該模型的參數(shù)設置為開放模式,可以根據(jù)參數(shù)豐富度進行重新構架,還可結合管理需求開展多層次的面源污染綜合評估分析。該模型模擬的指標為總氮、總磷、氨氮和化學需氧量,與管理部門關注的指標相一致。此外,DPeRS模型實現(xiàn)了像元尺度多類型多指標面源污染負荷的空間可視化,可直觀提供面源污染的關鍵源匯區(qū),與傳統(tǒng)總量減排核算方法相比,實現(xiàn)了由“點”到“面”的突破,可為科學制定面源污染防治措施提供技術支撐。Wang 等[21]采用DPeRS 模型對2010年新安江流域面源污染量進行模擬,將模擬結果與2010年黃山市污染源調(diào)查結果進行對比分析,結果表明總氮、總磷和氨氮指標的平均模擬精度達到0.75;并將DPeRS模型與SWAT模型的模擬結果進行對比分析,結果表明兩種模型總磷的模擬結果擬合度很好(R2為0.82)。為保證小流域的閉合性和水文數(shù)據(jù)的準確性,本研究結合可獲取的水文站監(jiān)測數(shù)據(jù),僅選用以面源污染為主的蘆村河斷面所在的小流域作為驗證流域,結果表明總磷面源污染的模型模擬精度達到90.31%,但不同驗證流域的監(jiān)測結果所評估的模擬精度會有所差異,因此,拓展驗證流域范圍、深入開展DPeRS模型的區(qū)域適用性研究是今后的重點研究方向。

    3.2 結果對比分析

    李建勛等[31]利用輸出系數(shù)法對陜西省農(nóng)業(yè)面源總氮、總磷污染的排放量和排放負荷進行模擬,結果表明咸陽市、渭南市總磷排放負荷較高,結果與本研究中的總磷排放高值區(qū)相吻合,這與其地處關中平原、農(nóng)業(yè)產(chǎn)業(yè)化程度較高、規(guī)?;笄蒺B(yǎng)殖發(fā)達等因素息息相關,且其數(shù)值與本研究結果在同一數(shù)量級。DPeRS模型模擬的2018年渭河流域典型斷面匯水區(qū)面源污染結果,總磷污染排放量和入河量分別為2 697t和512 t,排放強度在0.09~1.63 t·km-2間波動[32],與本研究的渭河流域排放強度波動區(qū)間具有高度一致性,二者的差值是可接受的,這與模型選取、輸入數(shù)據(jù)差異、斷面選取等因素有關。

    3.3面源污染時空變化特征的影響因素分析

    降雨沖刷和地表徑流是面源污染的主要驅(qū)動力,是面源污染發(fā)生的動力和輸移條件的載體,故研究區(qū)面源污染總磷排放負荷和入河負荷的時空變化特征與降水特征之間存在一定的相關性。在年尺度上,對二者進行相關性分析,結果發(fā)現(xiàn),研究區(qū)面源污染總磷排放負荷和入河負荷與年降水量的Pearson相關系數(shù)分別為0.57和0.47,在0.01水平上顯著相關。在月尺度上,在研究區(qū)降水量集中的6—8月,5年面源污染總磷平均排放量和入河量占比分別為52.15% 和53.38%,因此每年的6—8月是研究區(qū)面源污染總磷防治的關鍵期。

    2016—2020年渭河流域(陜西段)面源污染總磷的主要類型為水土流失型和農(nóng)田徑流型,除2016年兩種類型總磷排放量占比相當外,其余年份均以水土流失型為主,2018年水土流失型排放量占比較為突出(占比為82.48%),且草地面源污染較為突出的區(qū)域(西北部的榆林段和延安段)年降水量也相對較高(超過1 000 mm),2020年降水量超過600 mm的區(qū)域面積占流域總面積的40.18%,故2018年和2020年流域面源污染總磷排放負荷和排放量均較為突出。2016—2020年流域面源污染總磷入河負荷和入河量整體表現(xiàn)為先上升后下降的趨勢(2018年為折點),且入河負荷高值區(qū)主要集中分布在流域西北部的榆林段和延安段,這與榆林段和延安段的年降水量變化趨勢相一致,流域內(nèi)各區(qū)域的年降水量變化詳見圖7。此外,農(nóng)田磷平衡量是農(nóng)業(yè)面源污染的初始源強,其與農(nóng)田徑流型排放高度相關,2016年研究區(qū)農(nóng)田磷平衡量相對較高(流域均值達到34.38 t·km-2),故2016年研究區(qū)農(nóng)田徑流型總磷排放量占比最高。

    氣候變化影響區(qū)域氣候要素與水文要素的變化,全球變暖背景下,氣候氣象要素變化的不穩(wěn)定性加劇,氣候變化帶來的次生災害增加會導致水土流失加劇,同時區(qū)域降雨變化也會造成大量營養(yǎng)鹽被沖刷后隨地表徑流進入水體,加速受納水體的富營養(yǎng)化。渭河流域?qū)儆跉夂蛎舾袇^(qū)和生態(tài)環(huán)境脆弱區(qū),隨著氣候變暖、水循環(huán)加快,流域的降水特征也會發(fā)生變化,從而導致流域地表徑流和洪水災害的增加,流域6—8月降水量較為集中,且極端降水主要集中在7、8月[33],而總磷面源污染的產(chǎn)生和輸送受降雨沖刷和地表徑流的直接影響,因此,渭河流域關中平原地區(qū)重點應從源頭上防范豐水期面源污染物排放,尤其是渭南市西部和咸陽市中部的面源污染總磷排放負荷高值區(qū),建議從農(nóng)田和果園養(yǎng)分管理措施方面入手,從節(jié)水、節(jié)肥等方面提高農(nóng)業(yè)生產(chǎn)效率,保證土壤養(yǎng)分供給的前提下減少養(yǎng)分流失,進而降低農(nóng)田生產(chǎn)和果園種植過程導致的面源污染總磷流失;流域西北部的榆林市和延安市的草地分布地帶建議在7—8月加強水土保持工作,以減少草地水土流失引發(fā)的土壤中顆粒態(tài)磷的遷移運動,以及隨地表徑流進入水體的污染問題。

    4 結論

    (1)流域總磷面源污染負荷的時空變化特征與降水特征有較強的相關性,同時也與區(qū)域土地利用分布特征有關。2016—2020年渭河流域(陜西段)面源污染總磷入河負荷表現(xiàn)為先上升后下降的趨勢(2018年為折點),這與流域西北部負荷高值區(qū)因降水量變化引發(fā)的草地水土流失面源污染有關。6—8月流域面源污染總磷排放量和入河量較大,是研究區(qū)面源污染總磷防治的關鍵期。

    (2)對于渭河流域(陜西段)面源污染總磷負荷高值區(qū),流域西北部的榆林市和延安市以水土流失型為主要污染類型,建議進一步加強水土保持工作以減少草地水土流失引發(fā)的面源污染排放;關中平原的渭南市西部及咸陽市中部的種植區(qū)以農(nóng)田徑流型為主要污染類型,建議加強農(nóng)田養(yǎng)分管理以減少養(yǎng)分流失。

    (3)流域3個重點匯水區(qū)中,臨河入渭斷面所在匯水區(qū)的面源污染總磷主導類型為農(nóng)田徑流型,建議以節(jié)水、節(jié)肥等控源措施為主,岔口和新河入渭斷面所在匯水區(qū)主導類型為水土流失型,其中,岔口斷面所在匯水區(qū)2017—2019年總磷面源污染對水質(zhì)總磷指標的平均貢獻率達到77.38%,建議該匯水區(qū)上中游地區(qū)主要關注林地水土保持、下游地區(qū)重點從農(nóng)田養(yǎng)分管理加強面源污染防控。

    猜你喜歡
    面源污染渭河流域
    寧夏海原縣農(nóng)業(yè)面源污染防治工作調(diào)研
    農(nóng)村生活性面源污染治理的制約因素
    中國市場(2017年2期)2017-02-28 21:06:42
    涇河流域陜西段水污染源分析及污染防治對策
    基于RS/GIS 渭河流域植被覆蓋時空變化特征研究
    燕太子回國
    麻栗坡縣農(nóng)業(yè)面源污染現(xiàn)狀分析及控制措施
    環(huán)境流域治理中的法律對策
    法制博覽(2015年10期)2015-11-06 12:16:49
    渭河流域香菜夏秋無公害栽培技術
    面源污染、碳排放雙重視角下中國耕地利用效率的時空分異
    生態(tài)系統(tǒng)服務支付意愿及其影響因素分析
    軟科學(2015年6期)2015-07-10 02:38:40
    搡老岳熟女国产| 国产精品综合久久久久久久免费| 国产视频一区二区在线看| 丁香欧美五月| 日韩三级视频一区二区三区| 精品熟女少妇八av免费久了| 精品久久久久久久久久免费视频| 精品不卡国产一区二区三区| 亚洲激情在线av| 日韩av在线大香蕉| 免费看日本二区| 国产精品电影一区二区三区| 老汉色∧v一级毛片| 日韩欧美一区二区三区在线观看| 搡老妇女老女人老熟妇| 特大巨黑吊av在线直播 | 久久天堂一区二区三区四区| 亚洲av电影不卡..在线观看| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品中文字幕在线视频| 老汉色av国产亚洲站长工具| tocl精华| 国产成人av激情在线播放| 日本成人三级电影网站| 一区二区日韩欧美中文字幕| 18禁美女被吸乳视频| av片东京热男人的天堂| 真人做人爱边吃奶动态| 91字幕亚洲| 亚洲熟妇熟女久久| 国产极品粉嫩免费观看在线| 国产1区2区3区精品| 国产伦在线观看视频一区| 老司机靠b影院| 久久亚洲精品不卡| 人成视频在线观看免费观看| 丰满人妻熟妇乱又伦精品不卡| 国产亚洲精品av在线| 国产一级毛片七仙女欲春2 | 精品一区二区三区av网在线观看| 免费在线观看影片大全网站| 国产亚洲欧美在线一区二区| 久久这里只有精品19| 国产精品电影一区二区三区| 久久婷婷成人综合色麻豆| 久热这里只有精品99| 性色av乱码一区二区三区2| 在线看三级毛片| 真人做人爱边吃奶动态| 999久久久国产精品视频| 一二三四社区在线视频社区8| 侵犯人妻中文字幕一二三四区| 人成视频在线观看免费观看| 午夜福利欧美成人| 亚洲精品av麻豆狂野| 丰满人妻熟妇乱又伦精品不卡| 国产激情久久老熟女| 国产日本99.免费观看| 深夜精品福利| 黄色视频不卡| 成人免费观看视频高清| 欧美黄色淫秽网站| 亚洲中文字幕一区二区三区有码在线看 | 桃色一区二区三区在线观看| 激情在线观看视频在线高清| 国产精品久久久人人做人人爽| 人人妻人人澡人人看| 黑人巨大精品欧美一区二区mp4| 国产午夜精品久久久久久| 午夜久久久在线观看| 亚洲av成人不卡在线观看播放网| a级毛片在线看网站| 黑人操中国人逼视频| 久久久久久久久久黄片| 男人操女人黄网站| 亚洲色图av天堂| 国产精品自产拍在线观看55亚洲| 高清在线国产一区| 国产成人欧美在线观看| 中文字幕精品免费在线观看视频| 欧美午夜高清在线| 婷婷丁香在线五月| 精品人妻1区二区| 一本一本综合久久| 国产黄片美女视频| 精品免费久久久久久久清纯| 国产精品野战在线观看| 777久久人妻少妇嫩草av网站| 亚洲无线在线观看| 给我免费播放毛片高清在线观看| 可以免费在线观看a视频的电影网站| 亚洲三区欧美一区| 少妇裸体淫交视频免费看高清 | 精品少妇一区二区三区视频日本电影| 欧美乱妇无乱码| 村上凉子中文字幕在线| 国产高清videossex| 色综合站精品国产| 欧美不卡视频在线免费观看 | 国产三级黄色录像| 99国产精品一区二区三区| 久久午夜亚洲精品久久| 欧美日本亚洲视频在线播放| 亚洲第一电影网av| 久久香蕉国产精品| 亚洲成人免费电影在线观看| 精品一区二区三区四区五区乱码| 一区福利在线观看| 国产激情欧美一区二区| 国产成人精品久久二区二区91| 色哟哟哟哟哟哟| 久久草成人影院| 亚洲精品美女久久久久99蜜臀| 搡老熟女国产l中国老女人| 亚洲欧美一区二区三区黑人| 亚洲精品在线美女| 十分钟在线观看高清视频www| 久久精品国产亚洲av香蕉五月| 窝窝影院91人妻| 最近最新免费中文字幕在线| 嫩草影视91久久| 亚洲第一电影网av| 国内毛片毛片毛片毛片毛片| 亚洲三区欧美一区| 日本免费一区二区三区高清不卡| 一本大道久久a久久精品| 成人18禁高潮啪啪吃奶动态图| av福利片在线| 99国产精品一区二区蜜桃av| 亚洲一码二码三码区别大吗| 国产精品久久久久久亚洲av鲁大| 久久久久久久久免费视频了| 久久亚洲真实| 天天躁狠狠躁夜夜躁狠狠躁| 男女那种视频在线观看| 亚洲av熟女| 最好的美女福利视频网| 国产99久久九九免费精品| a级毛片在线看网站| 天天躁狠狠躁夜夜躁狠狠躁| 国产av在哪里看| 亚洲第一av免费看| 午夜视频精品福利| 亚洲av五月六月丁香网| 日日爽夜夜爽网站| 午夜久久久在线观看| 国产真实乱freesex| 亚洲午夜精品一区,二区,三区| 男人操女人黄网站| 亚洲熟妇中文字幕五十中出| 国产亚洲精品久久久久久毛片| 免费看日本二区| 欧美激情高清一区二区三区| 国产精品久久久久久亚洲av鲁大| 久久 成人 亚洲| av福利片在线| 亚洲中文av在线| 每晚都被弄得嗷嗷叫到高潮| 50天的宝宝边吃奶边哭怎么回事| 丝袜人妻中文字幕| 亚洲精品久久成人aⅴ小说| 成年版毛片免费区| 精品久久久久久成人av| 国产成人精品无人区| 国产一区在线观看成人免费| 色婷婷久久久亚洲欧美| 狂野欧美激情性xxxx| 欧美性猛交╳xxx乱大交人| 99国产极品粉嫩在线观看| 热re99久久国产66热| 女警被强在线播放| 美女高潮到喷水免费观看| 亚洲全国av大片| 午夜日韩欧美国产| 亚洲欧美激情综合另类| 嫩草影院精品99| 午夜影院日韩av| 熟妇人妻久久中文字幕3abv| 免费看日本二区| 亚洲aⅴ乱码一区二区在线播放 | 777久久人妻少妇嫩草av网站| 午夜成年电影在线免费观看| 亚洲aⅴ乱码一区二区在线播放 | 脱女人内裤的视频| 麻豆国产av国片精品| 久久久久久久午夜电影| av天堂在线播放| 欧美乱色亚洲激情| 色av中文字幕| 一级片免费观看大全| 亚洲成人国产一区在线观看| 国产成人欧美| 色播在线永久视频| 久久草成人影院| 欧美黑人巨大hd| 日韩欧美一区二区三区在线观看| 精品不卡国产一区二区三区| 成人18禁高潮啪啪吃奶动态图| 制服诱惑二区| 亚洲久久久国产精品| 亚洲国产欧美日韩在线播放| 国产精品野战在线观看| 日韩一卡2卡3卡4卡2021年| 制服丝袜大香蕉在线| 国产不卡一卡二| 久久国产精品影院| 国内精品久久久久精免费| 午夜老司机福利片| 91av网站免费观看| 日韩精品免费视频一区二区三区| 黑人巨大精品欧美一区二区mp4| 50天的宝宝边吃奶边哭怎么回事| 国产麻豆成人av免费视频| 欧美日韩乱码在线| 久久婷婷人人爽人人干人人爱| 久久久久久免费高清国产稀缺| 国产精品电影一区二区三区| 日韩一卡2卡3卡4卡2021年| 国产成人精品久久二区二区91| 1024香蕉在线观看| 欧美一区二区精品小视频在线| 中文资源天堂在线| 欧美成人免费av一区二区三区| 18禁在线无遮挡免费观看视频 | 男人狂女人下面高潮的视频| 亚洲高清免费不卡视频| 国产成人福利小说| 精品久久久久久久久亚洲| 你懂的网址亚洲精品在线观看 | 91午夜精品亚洲一区二区三区| 特级一级黄色大片| 真实男女啪啪啪动态图| 最近2019中文字幕mv第一页| 国产乱人偷精品视频| 男女那种视频在线观看| 美女被艹到高潮喷水动态| av在线亚洲专区| 欧美激情久久久久久爽电影| 我的老师免费观看完整版| 看十八女毛片水多多多| 波多野结衣巨乳人妻| 亚洲美女搞黄在线观看 | 久久鲁丝午夜福利片| 久久久精品大字幕| 欧美成人精品欧美一级黄| 精品一区二区三区视频在线| 亚洲av熟女| 国国产精品蜜臀av免费| 国产高清激情床上av| 97碰自拍视频| 日韩欧美 国产精品| 丰满的人妻完整版| 亚洲精品粉嫩美女一区| 校园人妻丝袜中文字幕| 午夜老司机福利剧场| 午夜爱爱视频在线播放| 男人狂女人下面高潮的视频| 久久久久久久亚洲中文字幕| 国产视频一区二区在线看| 国产午夜福利久久久久久| 老熟妇乱子伦视频在线观看| 久久久欧美国产精品| 国产精品一区二区免费欧美| 老司机福利观看| 99热全是精品| 美女cb高潮喷水在线观看| а√天堂www在线а√下载| 国产成人福利小说| 99热网站在线观看| 热99在线观看视频| 欧美又色又爽又黄视频| 成年女人看的毛片在线观看| 大香蕉久久网| 亚洲国产日韩欧美精品在线观看| 精品一区二区三区av网在线观看| 色av中文字幕| 啦啦啦啦在线视频资源| 亚洲人与动物交配视频| 亚洲精品国产av成人精品 | 国产色爽女视频免费观看| 男人狂女人下面高潮的视频| 最近2019中文字幕mv第一页| 给我免费播放毛片高清在线观看| 精品免费久久久久久久清纯| 99在线人妻在线中文字幕| 身体一侧抽搐| 成人特级黄色片久久久久久久| 欧美激情久久久久久爽电影| 国产午夜精品论理片| 麻豆一二三区av精品| 天堂网av新在线| 亚洲国产精品国产精品| 成人亚洲精品av一区二区| 欧美中文日本在线观看视频| 一级黄片播放器| 日韩 亚洲 欧美在线| 国产麻豆成人av免费视频| 国产伦精品一区二区三区四那| 日日啪夜夜撸| 超碰av人人做人人爽久久| 激情 狠狠 欧美| 国产精品日韩av在线免费观看| 国产精品国产三级国产av玫瑰| 日韩精品有码人妻一区| www日本黄色视频网| 亚洲精品日韩av片在线观看| 国产三级在线视频| 亚洲真实伦在线观看| 亚洲精品在线观看二区| 国产精品久久久久久亚洲av鲁大| ponron亚洲| 国产免费一级a男人的天堂| av.在线天堂| 国产午夜福利久久久久久| 亚洲第一区二区三区不卡| av.在线天堂| 日本-黄色视频高清免费观看| 熟女电影av网| 欧美三级亚洲精品| 免费不卡的大黄色大毛片视频在线观看 | 三级国产精品欧美在线观看| 亚洲第一区二区三区不卡| 精品熟女少妇av免费看| 亚洲欧美精品自产自拍| 人人妻人人澡人人爽人人夜夜 | 99久国产av精品| 小蜜桃在线观看免费完整版高清| 亚洲人成网站高清观看| 男女那种视频在线观看| 精品久久久久久久久久免费视频| 国产高清三级在线| 人妻丰满熟妇av一区二区三区| 欧美一级a爱片免费观看看| 99热精品在线国产| 欧美成人精品欧美一级黄| 国产真实乱freesex| 国产精品一区www在线观看| 91在线精品国自产拍蜜月| 精品久久久久久久久av| 午夜福利在线在线| 日日啪夜夜撸| 亚洲国产精品sss在线观看| 少妇被粗大猛烈的视频| 婷婷六月久久综合丁香| 国产蜜桃级精品一区二区三区| 哪里可以看免费的av片| 亚洲成人久久爱视频| 国产成人freesex在线 | 成年女人永久免费观看视频| 国产一区二区在线av高清观看| 久久久久久久午夜电影| 欧美性猛交╳xxx乱大交人| 欧美日韩乱码在线| 亚洲国产精品成人综合色| 国产成人a∨麻豆精品| 俺也久久电影网| 成人av在线播放网站| 久久99热6这里只有精品| 亚洲一区高清亚洲精品| 欧美极品一区二区三区四区| 亚洲久久久久久中文字幕| 国产aⅴ精品一区二区三区波| 亚洲人成网站高清观看| 午夜视频国产福利| 国产视频一区二区在线看| 国产精品一区二区性色av| 波多野结衣高清无吗| 精品久久久久久久久亚洲| 久久久久久伊人网av| 少妇的逼好多水| 亚洲人与动物交配视频| 一个人看视频在线观看www免费| 国产真实伦视频高清在线观看| 我要搜黄色片| 成年女人毛片免费观看观看9| 欧美人与善性xxx| 无遮挡黄片免费观看| 欧美日韩乱码在线| 成人av在线播放网站| 成人亚洲精品av一区二区| 亚洲欧美日韩高清在线视频| 午夜免费男女啪啪视频观看 | 国产精品人妻久久久影院| 中国国产av一级| 成人无遮挡网站| 一进一出抽搐动态| 精品午夜福利视频在线观看一区| 亚洲第一电影网av| 赤兔流量卡办理| 色综合色国产| ponron亚洲| 久久精品久久久久久噜噜老黄 | 波多野结衣巨乳人妻| 亚洲av中文av极速乱| 听说在线观看完整版免费高清| 99热全是精品| 亚洲,欧美,日韩| 搡老岳熟女国产| 男女视频在线观看网站免费| 大又大粗又爽又黄少妇毛片口| 精品不卡国产一区二区三区| 自拍偷自拍亚洲精品老妇| 我的女老师完整版在线观看| 我要看日韩黄色一级片| 18禁裸乳无遮挡免费网站照片| 欧美色欧美亚洲另类二区| 老师上课跳d突然被开到最大视频| av专区在线播放| 欧美绝顶高潮抽搐喷水| 日本黄色视频三级网站网址| 久久欧美精品欧美久久欧美| 美女免费视频网站| 亚州av有码| 我要搜黄色片| 亚洲欧美清纯卡通| 夜夜夜夜夜久久久久| 精品人妻一区二区三区麻豆 | 欧美一级a爱片免费观看看| 国产亚洲精品综合一区在线观看| 色综合色国产| 99久久久亚洲精品蜜臀av| 久久精品国产自在天天线| 99在线人妻在线中文字幕| 天天躁日日操中文字幕| 亚洲高清免费不卡视频| 欧美色视频一区免费| 亚洲国产精品国产精品| 日韩大尺度精品在线看网址| 欧洲精品卡2卡3卡4卡5卡区| 久久久久免费精品人妻一区二区| 成人特级黄色片久久久久久久| 精华霜和精华液先用哪个| 亚洲精品国产成人久久av| 99riav亚洲国产免费| 成人亚洲精品av一区二区| 男女那种视频在线观看| 久久久久久久久久成人| 国产高清视频在线播放一区| 国产精品一及| 亚洲国产精品成人综合色| 又爽又黄无遮挡网站| 亚洲中文日韩欧美视频| 寂寞人妻少妇视频99o| 日日撸夜夜添| 欧美zozozo另类| 成年免费大片在线观看| 国产色婷婷99| av.在线天堂| 免费av毛片视频| 免费人成视频x8x8入口观看| 又爽又黄无遮挡网站| 国产精品人妻久久久久久| 99久久久亚洲精品蜜臀av| 日韩欧美三级三区| 久久天躁狠狠躁夜夜2o2o| 老司机午夜福利在线观看视频| 在线播放无遮挡| 久久精品人妻少妇| 欧美精品国产亚洲| 给我免费播放毛片高清在线观看| 欧美国产日韩亚洲一区| 国产精品一区二区三区四区免费观看 | 欧美日韩乱码在线| 午夜爱爱视频在线播放| 国产综合懂色| 亚洲精品一卡2卡三卡4卡5卡| 看非洲黑人一级黄片| 最近中文字幕高清免费大全6| 成年免费大片在线观看| 亚洲国产日韩欧美精品在线观看| 日韩一本色道免费dvd| 成年女人永久免费观看视频| 亚洲国产色片| av专区在线播放| 在线观看av片永久免费下载| 十八禁网站免费在线| 欧美色视频一区免费| 人妻丰满熟妇av一区二区三区| 久久精品91蜜桃| 国产一区亚洲一区在线观看| 国产精华一区二区三区| 男女那种视频在线观看| 成年女人看的毛片在线观看| 国产成年人精品一区二区| 午夜精品在线福利| 亚洲三级黄色毛片| 欧美一区二区精品小视频在线| 床上黄色一级片| 亚洲人与动物交配视频| 国产精品人妻久久久久久| av在线观看视频网站免费| 免费看美女性在线毛片视频| 亚洲精品456在线播放app| 欧美激情久久久久久爽电影| 精品国产三级普通话版| 国产精品久久视频播放| 国产综合懂色| 我的老师免费观看完整版| 国产精品三级大全| aaaaa片日本免费| 搞女人的毛片| 成人国产麻豆网| 久久久久九九精品影院| 国内精品久久久久精免费| 22中文网久久字幕| 听说在线观看完整版免费高清| 男人狂女人下面高潮的视频| 色综合色国产| 美女黄网站色视频| 精品一区二区三区视频在线| 国产伦精品一区二区三区四那| 日产精品乱码卡一卡2卡三| 一夜夜www| 日本在线视频免费播放| 亚洲第一区二区三区不卡| 国产一区二区激情短视频| 夜夜夜夜夜久久久久| 此物有八面人人有两片| 美女 人体艺术 gogo| 精品一区二区三区视频在线观看免费| 两个人视频免费观看高清| 久久久久久久久中文| 国产私拍福利视频在线观看| 久久草成人影院| 国产三级中文精品| 春色校园在线视频观看| 婷婷精品国产亚洲av| 国产一区二区激情短视频| 蜜桃亚洲精品一区二区三区| 插逼视频在线观看| 国国产精品蜜臀av免费| 人妻久久中文字幕网| 成年女人看的毛片在线观看| 亚洲最大成人av| 直男gayav资源| 日韩av不卡免费在线播放| 日日摸夜夜添夜夜爱| 午夜精品一区二区三区免费看| 男女啪啪激烈高潮av片| 狠狠狠狠99中文字幕| 丝袜喷水一区| 2021天堂中文幕一二区在线观| 国产精品久久久久久av不卡| 成人一区二区视频在线观看| 搡老岳熟女国产| 久久午夜亚洲精品久久| 欧美激情国产日韩精品一区| 免费电影在线观看免费观看| 国产aⅴ精品一区二区三区波| 久久精品久久久久久噜噜老黄 | 色噜噜av男人的天堂激情| 免费在线观看成人毛片| 国产精品一区二区性色av| 插阴视频在线观看视频| 国语自产精品视频在线第100页| 一级毛片久久久久久久久女| 午夜视频国产福利| 精品久久国产蜜桃| 国产伦精品一区二区三区四那| www.色视频.com| 日本黄色片子视频| 亚洲中文字幕一区二区三区有码在线看| 深爱激情五月婷婷| 亚洲在线观看片| 3wmmmm亚洲av在线观看| 熟女人妻精品中文字幕| 波野结衣二区三区在线| 亚洲性夜色夜夜综合| 日本一本二区三区精品| 亚洲自拍偷在线| 国产片特级美女逼逼视频| 看十八女毛片水多多多| 日韩强制内射视频| 99热这里只有是精品50| 精华霜和精华液先用哪个| 美女免费视频网站| 国产精品一区www在线观看| .国产精品久久| 免费大片18禁| 国产探花在线观看一区二区| 欧美bdsm另类| 久久久久性生活片| 国产免费一级a男人的天堂| 免费搜索国产男女视频| 国产亚洲精品av在线| 十八禁网站免费在线| 亚洲不卡免费看| 国产精品无大码| 热99re8久久精品国产| 亚洲第一区二区三区不卡| 69av精品久久久久久| 日韩成人伦理影院| 九色成人免费人妻av| 国产aⅴ精品一区二区三区波| av国产免费在线观看| 少妇熟女欧美另类| 悠悠久久av| 久久久久九九精品影院| 少妇高潮的动态图| 在线观看美女被高潮喷水网站| a级一级毛片免费在线观看| 亚洲av.av天堂| 成熟少妇高潮喷水视频| 亚洲精品在线观看二区| 成人性生交大片免费视频hd| 你懂的网址亚洲精品在线观看 | 九九爱精品视频在线观看| 精品人妻视频免费看| 51国产日韩欧美| 国内精品美女久久久久久| 美女内射精品一级片tv| 国产欧美日韩一区二区精品| 搞女人的毛片| 美女xxoo啪啪120秒动态图| 亚洲中文日韩欧美视频|