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

    基于穩(wěn)定同位素和貝葉斯模型的引黃灌區(qū)地下水硝酸鹽污染源解析*

    2019-03-08 03:03:58張秋英李發(fā)東畢直磊
    關(guān)鍵詞:污染

    張 妍, 張秋英, 李發(fā)東, 張 鑫, 畢直磊, 張 強

    ?

    基于穩(wěn)定同位素和貝葉斯模型的引黃灌區(qū)地下水硝酸鹽污染源解析*

    張 妍1, 張秋英2**, 李發(fā)東3,4, 張 鑫1, 畢直磊1, 張 強1

    (1. 西北大學(xué)城市與環(huán)境學(xué)院陜西省地表系統(tǒng)與環(huán)境承載力重點實驗室 西安 710127; 2. 中國環(huán)境科學(xué)研究院 北京 100012; 3. 中國科學(xué)院地理科學(xué)與資源研究所生態(tài)系統(tǒng)網(wǎng)絡(luò)觀測與模擬重點實驗室 北京 100101; 4. 中國科學(xué)院大學(xué)資源與環(huán)境學(xué)院 北京 100190)

    地下水硝酸鹽(NO3-)污染已經(jīng)成為全球嚴重的水環(huán)境問題之一, 由于飲用水中高含量NO3-會轉(zhuǎn)化成亞硝酸鹽而增加各種疾病和癌癥風(fēng)險, 其來源的確定對于NO3-污染的預(yù)防和控制非常重要。本文以黃河下游第二大灌區(qū)——潘莊灌區(qū)為例, 首次采用NO3-的氮氧穩(wěn)定同位素結(jié)合貝葉斯模型追溯地下水NO3-的來源并量化各種來源的貢獻比例。結(jié)果表明, 地下水NO3-含量分布在0.1~197.0 mg?L-1, 平均值為34.2 mg·L-1。與《生活飲用水衛(wèi)生標準》中規(guī)定的地下水NO3-最大含量[20 mg(N)?L-1, 相當(dāng)于NO3-含量90 mg·L-1]相比, 有10%的樣品NO3-含量超標。井深<30 m、30~60 m和>60 m的地下水NO3-平均含量分別為25.9 mg·L-1、39.7 mg·L-1和20.1 mg·L-1??臻g上, 寧津縣、武城縣、平原縣和禹城市有大片區(qū)域地下水NO3-含量較高。地下水NO3-的δ15N組成范圍為0.72‰~23.93‰, 平均值為11.62‰; δ18O組成范圍為0.49‰~22.50‰, 平均值為8.46‰。同位素結(jié)果表明糞便和污水、農(nóng)業(yè)化肥是地下水中NO3-的主要污染來源。這反映了人類活動是引起地下水NO3-污染的主要原因。貝葉斯模型結(jié)果顯示, 糞便和污水對潘莊灌區(qū)地下水中NO3-平均貢獻率高達56.2%, 化肥的平均貢獻率為19.3%, 大氣降水和土壤的平均貢獻率分別為6.2%和12.3%。由于污水、糞便和化肥是地下水中NO3-的主要來源, 為保護和改善研究區(qū)地下水水質(zhì), 建議加強污水管道建設(shè), 強化畜禽糞便的管理以及提高化肥利用效率。

    氮污染;地下水污染; 硝酸鹽; 同位素; 貝葉斯模型

    地下水是主要的供水資源, 過量的硝酸鹽(NO3-)含量對飲用水安全構(gòu)成巨大威脅, 直接影響地下水資源作為供水資源的可持續(xù)性[1]。當(dāng)飲用水中NO3-濃度達90~140 mg?L-1時, 在人體中經(jīng)硝酸鹽還原菌的作用, NO3-被還原生成亞硝酸鹽, 而亞硝酸鹽可與血液中的血紅蛋白結(jié)合形成高鐵血紅蛋白, 引起嬰兒患高鐵血紅蛋白癥而發(fā)生窒息, 也可導(dǎo)致孕婦自發(fā)性流產(chǎn)[2]。長期飲用高NO3-含量的水, 也可引發(fā)一些其他的健康問題, 如增加患膀胱癌和卵巢癌的風(fēng)險[3]。高濃度的NO3-不僅威脅著人類的身體健康, 也會加重水體的富營養(yǎng)化[4]。地下水一旦受到NO3-污染, 治理難度很大且費用很高, 最根本的辦法就是找到NO3-的污染源頭, 從源頭控制和減少NO3-向地下水輸送, 這對控制和治理地下水NO3-污染具有非常重要的意義。

    地下水中NO3-來源復(fù)雜多樣, 識別其來源的方法有很多, 最簡單和傳統(tǒng)的方法是通過調(diào)查污染區(qū)的土地利用類型并結(jié)合水化學(xué)特征來分析污染源[5-6]。然而, 由于氮循環(huán)中普遍存在復(fù)雜的物理和化學(xué)作用以及生物轉(zhuǎn)化過程, 傳統(tǒng)方法難以準確判斷NO3-的來源及其經(jīng)歷的地球化學(xué)循環(huán)過程。氮氧同位素技術(shù)的出現(xiàn)為地下水NO3-污染的來源追溯提供了一個重要的方法。由于不同來源的NO3-具有不同的同位素組成, δ15N和δ18O雙同位素在示蹤地下水NO3-不同來源[7-10]、評價硝化/反硝化過程[11-12]等方面得到了廣泛應(yīng)用。劉君等[10]采用NO3-的δ15N和δ18O相結(jié)合的方法對石家莊市地下水中的NO3-來源、反硝化作用的發(fā)生進行了識別, 判定其NO3-主要來源于當(dāng)?shù)氐幕屎蛣游锛S肥。Liu等[11]利用δ15N和δ18O技術(shù)研究了貴陽市地下水NO3-的來源, 發(fā)現(xiàn)夏季城郊地下水NO3-主要來源于化肥和含氮有機物的硝化作用。通過δ15N和δ18O可以確定地下水中NO3-的主要污染來源及其地球化學(xué)循環(huán), 但不能精確地知道每種污染源的貢獻, 因此需要借助模型來進一步估算各種來源對地下水中NO3-的貢獻比例[13]。貝葉斯模型是基于狄利克雷分布, 在貝葉斯框架下構(gòu)建的一個邏輯先驗分布模型, 綜合考慮了同位素的時空變異性、反硝化反應(yīng)中同位素分餾作用、NO3-的多個污染來源等因素[14]。雙同位素數(shù)據(jù)結(jié)合貝葉斯模型可以精確量化不同污染源對地下水NO3-的貢獻比例[15]。

    黃河下游引黃灌區(qū)位于華北平原腹地, 人類活動歷史悠久, 經(jīng)歷了長期的引黃灌溉和施肥耕作, 氮肥施用量逐年增加, 但是農(nóng)業(yè)活動中氮肥的利用率低[9], 導(dǎo)致地下水NO3-含量增加[16-17], 且工業(yè)及居民生活排污也對地下水造成了巨大的環(huán)境污染[13]。目前, 華北平原地下水NO3-污染的嚴峻現(xiàn)狀得到較多關(guān)注[9,18]。黃河下游地區(qū)地下水NO3-平均含量高達45.3 mg?L-1[19], 黃河三角洲地區(qū)地下水中NO3-平均含量甚至高達101 mg?L-1[20]。潛水和承壓水中NO3-含量迅速上升會造成深層地下水污染[9]。然而, 目前還少有該區(qū)地下水NO3-污染來源的定量研究。本文選擇黃河下游第二大引黃灌區(qū)——潘莊灌區(qū)為研究區(qū), 應(yīng)用穩(wěn)定同位素技術(shù)結(jié)合水化學(xué)識別地下水NO3-的各種來源, 并采用貝葉斯模型量化各種來源對地下水NO3-的貢獻比例, 旨在為地下水資源管理和氮污染來源控制提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    潘莊引黃灌區(qū)位于華北平原山東省德州市西部(36.4°~38.0° N、115.75°~117.6° E), 是黃河下游第二大灌區(qū), 屬于全國大型灌區(qū)??偼恋孛娣e5.85×105hm2,其中耕地面積3.33×105hm2, 居民用地9.10×104hm2, 河道和坑塘面積6.80×104hm2, 其他占地面積9.30×104hm2[21]。灌溉范圍涉及山東省齊河、禹城、夏津、武城、平原、陵縣、德城、寧津8個縣市(圖1)。潘莊引黃灌區(qū)地處黃河沖積平原, 地勢平坦開闊, 自西南向東北渤海灣傾斜。地貌以黃河沖積平原為主, 屬于暖溫帶大陸性季風(fēng)氣候區(qū), 四季分明, 多年平均氣溫11.7~14 ℃, 多年平均蒸發(fā)量1 000~1 200 mm, 多年平均降雨量582 mm。降雨時空分布不均, 多集中在6—9月, 其中8月份降雨量最大。太陽輻射總量5 225 MJ?m-2, 日照時數(shù)2 640 h。光熱資源豐富, 雨熱同期, 利于農(nóng)業(yè)生產(chǎn)[13]。潘莊引黃灌區(qū)引水緩減了德州地區(qū)工農(nóng)業(yè)生產(chǎn)和城鄉(xiāng)人畜飲水等水資源短缺, 改良了鹽堿地, 補充了地下水資源, 促進了當(dāng)?shù)貒窠?jīng)濟的快速持續(xù)發(fā)展, 創(chuàng)造了巨大的社會經(jīng)濟效益。

    圖1 潘莊引黃灌區(qū)地下水采樣點(G1-G40)分布圖

    地下水主要賦存并運動于第四系及新第三系松散沉積物之孔隙中, 依據(jù)地下水的埋藏條件、水力性質(zhì)和水質(zhì)特征可劃分為3個主要含水層組。第1含水層組為淺層地下水, 埋深0~60 m, 含水層厚度10~20 m。第2含水層組為中層微承壓咸水, 埋深60~200 m, 含水層厚度20 m左右。第3含水層組為深層地下水, 埋深200~500 m, 含水層厚度15~40 m[22]。含水層巖性以粉細砂、中細砂、中粗砂為主, 單井出水量20~80 m3?h-1。因受區(qū)域隔水層及咸水體的阻隔, 深層地下水與淺層地下水基本無水力聯(lián)系。在天然狀態(tài)下, 淺層地下水補給來源主要是降水入滲補給, 其次是山前側(cè)向徑流流入補給、河道滲漏補給。深層地下水補給來源以水平徑流為主, 主要接受西部上游同層補給和南部基巖裂隙水的補給, 總體徑流方向由南西向北東運動, 平均水力坡度0.11‰~0.12‰。自20世紀60年代以來, 隨著地下水開采量的不斷增加, 研究區(qū)深層地下水超采漏斗相繼形成, 使地下水天然狀態(tài)下的流場發(fā)生了改變, 呈現(xiàn)出向漏斗中心匯流的特征[23]。

    潘莊引黃灌區(qū)農(nóng)墾歷史悠久, 土地利用率較高。研究區(qū)種植作物主要是夏玉米(L.)和冬小麥(L.)。主要使用的化肥種類有碳銨、尿素和緩釋肥等[13]。長期以來存在施用單元肥料和過量施肥等問題, 導(dǎo)致土壤的養(yǎng)分結(jié)構(gòu)發(fā)生較大變化, 養(yǎng)分比例失調(diào), 致使土壤肥力狀況仍然滿足不了農(nóng)作物產(chǎn)量日益增長的需要, 各種環(huán)境問題如NO3-淋失等問題凸顯。除了農(nóng)業(yè)區(qū)氮肥施用是水體NO3-的潛在污染來源, 潘莊引黃灌區(qū)內(nèi)還有大量不同種類的工廠企業(yè), 大部分企業(yè)污水未經(jīng)處理直接排入河道, 這些廢水中硝態(tài)氮和銨態(tài)氮平均含量分別為16.81 mg?L-1和0.55 mg?L-1[9,24]。

    1.2 樣品的采集與分析

    根據(jù)當(dāng)?shù)氐乃牡刭|(zhì)條件和《采樣方案設(shè)計技術(shù)規(guī)定》(GB12997—97)中規(guī)定的水質(zhì)量控制、質(zhì)量表征、污染物鑒別采樣方案的原則, 設(shè)置了水樣點網(wǎng)格布設(shè)。于2013年9月采集40個地下水樣品(圖1), 采集的地下水包括36個淺層地下水和4個中層承壓微咸水。對各個取樣點進行GPS定位, 并在現(xiàn)場測定pH、EC、溶解氧(DO)和氧化還原電位(ORP)。

    將全部水樣過0.45mmol?L-1醋酸纖維素濾膜, 然后用戴安離子色譜儀(ICS-900, 美國)測定陰離子(Cl-、NO3-、SO42-)含量, 采用ICP-OES分析陽離子(K+、Ca2+、Na+、Mg2+)含量, 用硫酸滴定法測定HCO3-含量。δ15N與δ18O經(jīng)“離子交換法”[25]前處理之后, 采用氣體穩(wěn)定同位素質(zhì)譜儀(MAT-253)測定。所有離子含量和δ15N與δ18O均在中國科學(xué)院地理科學(xué)與資源研究所理化分析中心測定。

    1.3 貝葉斯混合模型

    使用貝葉斯穩(wěn)定同位素混合模型計算各種污染來源的貢獻比例[26], 表達式為:

    式中:X表示混合物同位素的值;P表示來源的比例;S表示第個來源的第種同位素的值, 服從均值為方差的正態(tài)分布;c表示第個來源的同位素的分餾系數(shù), 服從均值為方差為的正態(tài)分布;是殘余誤差, 表示其他各個混合物間無法量化的方差, 其均值和標準差在通常情況下均為0。本文采用貝葉斯混合模型定量計算了4種潛在來源(大氣降水、土壤氮、化肥、污水和糞便)對地下水中NO3-的貢獻比率, 這4種來源的同位素值參照文獻[13]。

    2 結(jié)果與分析

    2.1 地下水的水文地球化學(xué)特征

    由表1可知, 潘莊引黃灌區(qū)地下水的pH含量范圍為6.99~8.28, 均值為7.48, 屬偏堿性水質(zhì)。EC的含量范圍為770~7 090mS?cm-1。地下水中陰離子以HCO3-離子為主, SO42-和Cl-次之; 陽離子以Na+離子為主, Ca2+離子次之。陰離子中HCO3-離子含量294.5~ 1 259.4 mg?L-1; 陽離子中Na+離子含量為48.7~627.2 mg?L-1。水組成反映了碳酸鹽巖溶解的影響。華北平原土壤鹽堿化嚴重, 鹽漬土中含有大量的NaHCO3和Na2HCO3, 隨灌溉和降水從土壤中溶解進入地表水和地下水[27]。DO含量范圍為0.01~6.70 mg?L-1, 反映了補給條件和反硝化發(fā)生的可能[28]。NO3-含量介于0.1~197.0 mg?L-1, 平均值為34.2 mg?L-1。從變異系數(shù)來看, 離子中NO3-變異系數(shù)最大, 達139.4%, 大于100%, 表明NO3-離子空間變異大, 可能是因為地下水流動過程中, 不同的水文地質(zhì)條件以及人類活動如農(nóng)田施肥、污水排放等導(dǎo)致地下水水化學(xué)特征的空間變異增大[29]。

    由圖2可以看出, EC、HCO3-、SO42-、Cl-、Ca2+、Na+、Mg2+、K+隨著地下水深的增加而出現(xiàn)先增加后減少的趨勢, 最大值出現(xiàn)在井深40~60 m。淺層地下水(<60 m)的NO3-含量高于中層地下水(>60 m), 最大值為197.0 mg?L-1, 出現(xiàn)在井深40~60 m。為分析探討地下水NO3-含量隨井深的變化規(guī)律, 依據(jù)所取地下水樣的埋藏深度, 將水樣劃分為3組: 井深<30 m、30~60 m和>60 m的地下水。不同井深的地下水NO3-平均含量存在差異。井深<30 m的地下水NO3-平均含量為 25.9 mg?L-1; 井深在30~60 m, NO3-平均含量最高, 達39.7 mg?L-1; 井深>60 m的地下水, NO3-平均含量最小, 為20.1 mg?L-1。

    表1 潘莊引黃灌區(qū)地下水水化學(xué)特征

    圖2 潘莊引黃灌區(qū)地下水EC和主要離子含量與井深的關(guān)系

    2.2 地下水NO3-空間分布特征

    由表1可以看出, 潘莊灌區(qū)地下水中存在NO3-污染, 有10%的樣點超過《生活飲用水衛(wèi)生標準》(GB5749—2006)中NO3-的上限值90 mg?L-1??傮w上看, 本研究區(qū)西北部地下水NO3-含量很低, 普遍小于45 mg?L-1; 而東北部NO3-含量較高, 其中, 寧津縣、武城縣、平原縣和陵縣有大片區(qū)域NO3-含量超過45 mg?L-1(圖3)。

    圖3 潘莊引黃灌區(qū)地下水NO3-含量的空間分布

    2.3 地下水NO3-氮氧同位素組成及各種來源貢獻率

    地下水NO3-的δ15N和δ18O組成范圍分別為0.72‰~23.93‰和0.49‰~22.50‰(圖4)。地下水中δ15N和δ18O與NO3-濃度之間的相關(guān)關(guān)系比較復(fù)雜, 不存在δ15N和δ18O值隨NO3-濃度減少而增加的趨勢, 說明大部分地下水樣品中NO3-的δ15N和δ18O值基本反映了源的同位素特征。

    圖4 潘莊引黃灌區(qū)地下水中δ15N和δ18O與NO3-的關(guān)系圖

    40個地下水樣品中, 有28個樣品氮氧同位素值落在糞便和污水區(qū)域內(nèi)(圖5), 說明潘莊灌區(qū)地下水中NO3-的主要來源是糞便和污水。有8個樣品的氮氧同位素值落在化肥區(qū)域內(nèi), 說明農(nóng)業(yè)灌區(qū)化肥的使用是地下水NO3-污染的另一個重要來源。由于化肥與糞便和污水在雙同位素分布圖上有大部分重疊, 農(nóng)田化肥施入經(jīng)土壤中礦化作用轉(zhuǎn)化為硝態(tài)氮, 進而淋溶到地下水。另外4個樣品的氮氧同位素值未落在劃定的區(qū)域內(nèi)(圖5), 這4個點NO3-的δ15N/δ18O比值均不為2∶1, 且這4個點的DO均大于0.2 mg?L-1, 不利于發(fā)生反硝化[30], 因此這些地下水樣點NO3-污染是大氣降水和污水共同作用的結(jié)果。

    圖5 潘莊引黃灌區(qū)地下水中NO3-的δ15N和δ18O關(guān)系圖(不同來源的δ15N和δ18O范圍來自文獻[9,13])

    利用貝葉斯混合模型計算了大氣降水、土壤氮、化肥、糞便和污水對潘莊灌區(qū)地下水NO3-污染的貢獻率。由圖6可知, 糞便和污水的貢獻率最大, 貢獻率介于35.1%~80.5%, 平均貢獻率為56.2%; 化肥的貢獻率僅次于糞便和污水, 介于4.9%~31.7%, 平均貢獻率為19.3%; 大氣降水和土壤的貢獻率較小, 貢獻率分別為6.2%和12.3%。

    圖6 潘莊引黃灌區(qū)地下水NO3-潛在來源對地下水NO3-的貢獻率

    3 討論與結(jié)論

    本研究利用氮氧同位素初步判斷出地下水NO3-的污染來源, 利用貝葉斯混合模型進一步揭示了大氣降水、土壤氮、化肥、糞便和污水對潘莊灌區(qū)地下水NO3-污染的貢獻率。NO3-的氮氧同位素結(jié)合貝葉斯模型結(jié)果表明, 糞便和污水的貢獻率最大, 平均貢獻率達56.2%。潘莊引黃灌區(qū)內(nèi)有大量不同種類的工廠企業(yè), 如印刷廠、造紙廠、橡塑鋁塑制品廠、鋼鐵廠、石油加工廠、機械加工廠等高污染企業(yè), 大部分工業(yè)污水未經(jīng)處理直接排入河道, 且這些企業(yè)采取晝關(guān)夜排的策略來躲避政府和環(huán)保部門的監(jiān)測[13,31], 所排放的工業(yè)廢水中污染物主要以氨態(tài)氮和總氮為主[32]。此外, 經(jīng)調(diào)查, 研究區(qū)居住地成片分布, 尚無完善的排污管道, 除了隨意排放的生活污水和生活垃圾外, 居民所用旱廁所容納的糞尿會攜帶大量營養(yǎng)物質(zhì)滲入土壤甚至進入地下水中[33]。研究區(qū)內(nèi)還有較多分散的養(yǎng)殖場, 以前畜禽廢棄物多作為有機肥施用于農(nóng)田, 近年來由于當(dāng)?shù)貏趧恿ν饬骱突适褂昧康脑龃? 畜禽糞便直接堆積。我國全年禽畜糞便的氮素總量約為16億t, 還田利用氮素量占50%, 15%的氮素通過揮發(fā)作用進入大氣, 13%的氮素堆置廢棄, 流失總量高達7.8×105萬t。因畜禽排泄物已遠遠超過配套農(nóng)田的消納量, 致使畜禽糞尿隨著徑流直接進入溝渠以及土壤[34]。NO3-具有很高的遷移率, 很容易從土壤中淋溶滲出進入地下影響地下水水質(zhì)[24]。

    化肥的貢獻率僅次于糞便和污水, 平均貢獻率為19.3%。禹城市、齊河縣部分地區(qū)也存在地下水NO3-高的現(xiàn)象。華北平原是我國主要的糧食主產(chǎn)區(qū), 地下水超采嚴重, 農(nóng)業(yè)活動中氮肥消耗量已達450~600 kg?hm-2, 遠遠超過作物生長所需氮量, 且高于世界其他國家的氮肥施用量[35]。研究表明, 肥料當(dāng)季利用率較低, 氮肥為30%~35%、磷肥為10%~20%、鉀肥為35%~50%[9]。為增加小麥和玉米產(chǎn)量, 華北平原化肥施用量逐年增長。據(jù)統(tǒng)計資料[36-37], 山東省化肥施用量從20世紀90年代至今呈穩(wěn)定增長的趨勢(圖7a)。1993年, 山東省化肥平均施用量為2.8×106t, 而到2013年, 該省化肥平均施用量增加到4.8×106t。目前, 山東省農(nóng)田施用的肥料仍以尿素、硝銨、磷酸二銨和亞硝銨等氮肥和復(fù)合肥為主。氮肥施用量占34.21%, 復(fù)合肥占45.53%, 磷肥占10.49%, 鉀肥占9.76%。潘莊灌區(qū)是山東省的重要農(nóng)業(yè)生產(chǎn)區(qū), 近年來的化肥施用量也呈逐年增加的趨勢(圖7b)。1993年, 該區(qū)域化肥施用總量為1.6×105t, 到2013年, 潘莊灌區(qū)化肥施用總量增長到2.9×105t, 復(fù)合肥增長最為明顯。潘莊引黃灌區(qū)化肥施用主要以氮肥和復(fù)合肥為主, 在2013年這兩種肥料施用量分別占化肥總施用量的43.6%和33.7%。與山東省肥料構(gòu)成相比, 潘莊灌區(qū)氮肥施用比例過高, 而過高的氮素投入和低的作物吸收率導(dǎo)致農(nóng)田土壤氮的大量累積。土壤剖面累積的硝態(tài)氮如果不及時被作物吸收利用, 在較強降雨和大量灌溉條件下, 將向下移動逐漸脫離作物根區(qū), 使其生物有效性降低, 同時可能淋失直接進入淺層地下水[38]。因此, 農(nóng)田氮隨降水和農(nóng)作物灌溉等流失是農(nóng)業(yè)灌區(qū)水體NO3-含量增加的一個主要原因。黃河下游地區(qū)普遍存在過量施肥、肥料利用率偏低、NO3-流失嚴重等現(xiàn)象, 對地下水水質(zhì)安全構(gòu)成了極大的威脅。

    圖7 1993—2013年山東省(a)和潘莊灌區(qū)(b)化肥施用量年變化

    本研究中, 潘莊灌區(qū)地下水NO3-濃度介于0.1~197.01 mg?L-1, 平均濃度為34.2 mg?L-1。有10%的地下水樣品超過了《生活飲用水衛(wèi)生標準》中規(guī)定的NO3-上限值90 mg?L-1。從地理位置上分析, 寧津縣、武城縣、平原縣和禹城市地下水NO3-含量均較高。井深為30~60 m的地下水NO3-平均含量(39.7 mg?L-1)高于井深<30 m和井深>60 m。地下水NO3-的δ15N和δ18O均值分別為11.62‰和8.46‰。貝葉斯混合模型計算結(jié)果表明, 地下水中約有56.2%的NO3-來自糞便和污水, 19.3%的NO3-來自于農(nóng)業(yè)化肥。為保護和改善研究區(qū)地下水水質(zhì), 控制地下水NO3-的進一步增長, 需要從污染源入手采取如下措施: 1)有關(guān)部門應(yīng)以現(xiàn)在最嚴格的水資源管理制度和限制納污紅線考核制度為契機, 嚴格監(jiān)管各排污口污水的不達標排放。同時, 加強污水管道建設(shè), 進一步提高污水處理率。2)加強畜禽糞便的管理, 糞便堆放處必須有防滲措施, 并妥善處理。同時鼓勵合理地將其轉(zhuǎn)化為有機肥。3)提高化肥的利用效率, 加強引黃灌區(qū)的化肥施用管理。

    [1] SAJEDI-HOSSEINI F, MALEKIAN A, CHOUBIN B, et al. A novel machine learning-based approach for the risk assessment of nitrate groundwater contamination[J]. Science of the Total Environment, 2018, 644: 954–962

    [2] DAN-HASSAN M A, OLASEHINDE P I, AMADI A N, et al. Spatial and temporal distribution of nitrate pollution in groundwater of Abuja, Nigeria[J]. International Journal of Chemistry, 2012, 4(3): 104

    [3] LIU J, JIANG L H, ZHANG C J, et al. Nitrate-nitrogen contamination in groundwater: Spatiotemporal variation and driving factors under cropland in Shandong Province, China[J]. Earth and Environmental Science, 2017, 82(1): 012059

    [4] 劉貫群, 周書玉, 黃修東, 等. 多種方法識別青島大沽河平原區(qū)地下水硝酸鹽污染來源[J]. 環(huán)境科學(xué)學(xué)報, 2017, 37(1): 347–356 LIU G Q, ZHOU S Y, HUANG X D, et al. Multiple methods to recognize sources of underground water nitrate contamination in plain area of Dagu River, Qingdao, China[J]. Acta Scientiae Circumstantiae, 2017, 37(1): 347–356

    [5] 劉光棟, 吳文良, 劉仲蘭, 等. 華北農(nóng)業(yè)高產(chǎn)糧區(qū)地下水面源污染特征及環(huán)境影響研究——以山東省桓臺縣為例[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報, 2005, 13(2): 125–129 LIU G D, WU W L, LIU Z L, et al. Characteristics and environmental impact of non-point pollution of groundwater under the high-yield farmlands of North China — A case study from Huantai County, Shandong Province[J]. Chinese Journal of Eco-Agriculture, 2005, 13(2): 125–129

    [6] 高新昊, 江麗華, 劉兆輝, 等. 山東省農(nóng)村地區(qū)地下水硝酸鹽污染現(xiàn)狀調(diào)查與評價[J]. 中國農(nóng)業(yè)氣象, 2011, 32(1): 89–93 GAO X H, JIANG L H, LIU Z H, et al. Investigation and evaluation on nitrate contamination status of groundwater in rural areas of Shandong Province[J]. Chinese Journal of Agrometeorology, 2011, 32(1): 89–93

    [7] FERNáNDEZ E, GRILLI A, ALVAREZ D, et al. Evaluation of nitrate levels in groundwater under agricultural fields in two pilot areas in central Chile: A hydrogeological and geochemical approach[J]. Hydrological Processes, 2017, 31(6): 1206–1224

    [8] STOEWER M M, KNLLER K, STUMPP C. Tracing freshwater nitrate sources in pre-alpine groundwater catchments using environmental tracers[J]. Journal of Hydrology, 2015, 524: 753–767

    [9] ZHANG Y, LI F D, ZHANG Q Y, et al. Tracing nitrate pollution sources and transformation in surface- and ground-waters using environmental isotopes[J]. Science of the Total Environment, 2014, 490: 213–222

    [10] 劉君, 陳宗宇. 利用穩(wěn)定同位素追蹤石家莊市地下水中的硝酸鹽來源[J]. 環(huán)境科學(xué), 2009, 30(6): 1602–1607 LIU J, CHEN Z Y. Using stable isotope to trace the sources of nitrate in groundwater in Shijiazhuang[J]. Environmental Science, 2009, 30(6): 1602–1607

    [11] LIU C Q, LI S L, LANG Y C, et al. Using δ15N- and δ18O-values to identify nitrate sources in karst ground water, Guiyang, Southwest China[J]. Environmental Science & Technology, 2006, 40(22): 6928–6933

    [12] 魯垠濤, 劉芳, 姚宏, 等. 北京密云水庫小流域地下水硝酸鹽污染來源示蹤[J]. 環(huán)境化學(xué), 2016, 35(1): 180–188 LU Y T, LIU F, YAO H, et al. Source analysis of nitrate pollution source in groundwater in a small watershed of Miyun Reservoir in Beijing[J]. Environmental Chemistry, 2016, 35(1): 180–188

    [13] ZHANG Y, SHI P, LI F D, et al. Quantification of nitrate sources and fates in rivers in an irrigated agricultural area using environmental isotopes and a Bayesian isotope mixing model[J]. Chemosphere, 2018, 208: 493–501

    [14] PARNELL A C, INGER R, BEARHOP S, et al. Source partitioning using stable isotopes: Coping with too much variation[J]. PLoS One, 2010, 5(3): e9672

    [15] KIM K H, YUN S T, MAYER B, et al. Quantification of nitrate sources in groundwater using hydrochemical and dual isotopic data combined with a Bayesian mixing model[J]. Agriculture, Ecosystems & Environment, 2015, 199: 369–381

    [16] SERIO F, MIGLIETTA P P, LAMASTRA L, et al. Groundwater nitrate contamination and agricultural land use: A grey water footprint perspective in southern Apulia Region (Italy)[J]. Science of the Total Environment, 2018, 645: 1425–1431

    [17] ZHAO T Q, XU H H, He Y X, et al. Agricultural non-point nitrogen pollution control function of different vegetation types in riparian wetlands: A case study in the Yellow River wetland in China[J]. Journal of Environmental Sciences, 2009, 21(7): 933–939

    [18] JU X T, KOU C L, ZHANG F S, et al. Nitrogen balance and groundwater nitrate contamination: Comparison among three intensive cropping systems on the North China Plain[J]. Environmental Pollution, 2006, 143(1): 117–125

    [19] CHEN J Y. Holistic assessment of groundwater resources and regional environmental problems in the North China Plain[J]. Environmental Earth Sciences, 2010, 61(5): 1037–1047

    [20] CHEN J Y, TANIGUCHI M, LIU G Q, et al. Nitrate pollution of groundwater in the Yellow River delta, China[J]. Hydrogeology Journal, 2007, 15(8): 1605–1614

    [21] 王學(xué)軍, 劉慧霞, 李元波. 潘莊灌區(qū)續(xù)建配套和節(jié)水改造現(xiàn)狀與成效[J]. 水利科技與經(jīng)濟, 2012, 18(6): 63–64 WANG X J, LIU H X, LI Y B, et al. Present situation and effect of continued construction and water-saving transformation in the Panzhuang Irrigation Area[J]. Water Conservancy Science and Technology and Economy, 2012, 18(6): 63–64

    [22] 張芳, 李發(fā)東, 李靜, 等. 德州引黃灌區(qū)主要河系水化學(xué)空間特征分析[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報, 2013, 21(4): 487–493 ZHANG F, LI F D, LI J, et al. Hydrochemical characteristics of surface water in main rivers of the irrigation districts in the downstream of Yellow River[J]. Chinese Journal of Eco-Agriculture, 2013, 21(4): 487–493

    [23] 劉桂儀. 魯北平原深層地下水開發(fā)與環(huán)境問題[J]. 水文地質(zhì)工程地質(zhì), 2001, 28(3): 43–45 LIU G Y. The exploitation and environmental problems of deep groundwater in northern Shandong Plain[J]. Hydrogeology & Engineering Geology, 2001, 28(3): 43–45

    [24] JIA L, WANG W Y, LI Y H, et al. Heavy metals in soil and crops of an intensively farmed area: A case study in Yucheng City, Shandong Province, China[J]. International Journal of Environmental Research and Public Health, 2010, 7(2): 395–412

    [25] SILVA S R, GING P B, LEE R W, et al. Forensic applications of nitrogen and oxygen isotopes in tracing nitrate sources in urban environments[J]. Environmental Forensics, 2002, 3(2): 125–130

    [26] PARNELL A C, INGER R, BEARHOP S, et al. Source partitioning using stable isotopes: Coping with too much variation[J]. PLoS One, 2010, 5(3): e9672

    [27] GHAZARYAN K, CHEN Y N. Hydrochemical assessment of surface water for irrigation purposes and its influence on soil salinity in Tikanlik oasis, China[J]. Environmental Earth Sciences, 2016, 75(5): 383

    [28] PI?A-OCHOA E, áLVAREZ-COBELAS M. Denitrification in aquatic environments: A cross-system analysis[J]. Biogeochemistry, 2006, 81(1): 111–130

    [29] 蔣勇軍, 吳月霞, GROVES C, 等. 利用因子分析確定巖溶地下河系統(tǒng)水質(zhì)的影響因素[J]. 水文地質(zhì)工程地質(zhì), 2009, 36(4): 1–7 JIANG Y J, WU Y X, GROVES C, et al. Identification of influencing factors of groundwater quality variation using factor analysis[J]. Hydrogeology & Engineering Geology, 2009, 36(4): 1–7

    [30] 賀國平, 劉培斌, 慕星, 等. 永定河沖洪積扇地下水中硝酸鹽來源的同位素識別[J]. 水利學(xué)報, 2016, 47(4): 582–588 HE G P, LIU P B, MU X, et al. Identification of nitrate sources in groundwater in the Yongding River alluvial fan with isotope technology[J]. Journal of Hydraulic Engineering, 2016, 47(4): 582–588

    [31] ZHANG Y, LI F D, ZHAO G S, et al. An attempt to evaluate the recharge source and extent using hydrogeochemistry and stable isotopes in North Henan Plain, China[J]. Environmental Monitoring and Assessment, 2014, 186(8): 5185–5197

    [32] 張琳, 高翔. 我國廢水中主要污染物排放情況分析[J]. 統(tǒng)計學(xué)與應(yīng)用, 2018, 7(2): 134–140 ZHANG L, GAO X. Analysis on the discharge of main pollutants from wastewater in China[J]. Statistics and Application, 2018, 7(2): 134–140

    [33] 閆亞男, 張列宇, 席北斗, 等. 改良化糞池/地下土壤滲濾系統(tǒng)處理農(nóng)村生活污水[J]. 中國給水排水, 2011, 27(10): 69–72 YAN Y N, ZHANG L Y, XI B D, et al. Treatment of rural domestic sewage by modified septic tanks/subsurface soil infiltration system[J]. China Water & Wastewater, 2011, 27(10): 69–72

    [34] 石鵬飛, 鄭媛媛, 趙平, 等. 華北平原種養(yǎng)一體規(guī)?;r(nóng)場氮素流動特征及利用效率——以河北津龍循環(huán)農(nóng)業(yè)園區(qū)為例[J]. 應(yīng)用生態(tài)學(xué)報, 2017, 28(4): 1281–1288 SHI P F, ZHENG Y Y, ZHAO P, et al. Nitrogen flow characteristic and use efficiency in mixed crop-livestock system in North China Plain: Hebei Jinlong Circular Agriculture Farm as an example[J]. Chinese Journal of Applied Ecology, 2017, 28(4): 1281–1288

    [35] 王仕琴, 鄭文波, 孔曉樂. 華北農(nóng)區(qū)淺層地下水硝酸鹽分布特征及其空間差異性[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報, 2018, 26(10): 1476–1482 WANG S Q, ZHENG W B, KONG X L. Spatial distribution characteristics of nitrate in shallow groundwater of the agricultural area of the North China Plain[J]. Chinese Journal of Eco-Agriculture, 2018, 26(10): 1476–1482

    [36] 山東省統(tǒng)計局. 山東統(tǒng)計年鑒1993—2013[M]. 北京: 中國統(tǒng)計出版社, 1993—2013 Shandong Provincial Bureau of Statistics. Statistical Yearbook of Shandong 1993-2013[M]. Beijing: China Statistics Press, 1993–2013

    [37] 德州市統(tǒng)計局. 德州統(tǒng)計年鑒1993—2013[M]. 德州: 德州市統(tǒng)計局, 1993—2013 Dezhou Bureau of Statistics. Dezhou Statistical Yearbook 1993-2013[M]. Dezhou: Dezhou Bureau of Statistics, 1993-2013

    [38] 楊靜, 肖天昀, 李海波, 等. 江漢平原地下水中硝酸鹽的分布及影響因素[J]. 中國環(huán)境科學(xué), 2018, 38(2): 710–718 YANG J, XIAO T Y, LI H B, et al. Spatial distribution and influencing factors of the NO3-N concentration in groundwater in Jianghan Plain[J]. China Environmental Science, 2018, 38(2): 710–718

    Source identification of nitrate contamination of groundwater in Yellow River Irrigation Districts using stable isotopes and Bayesian model*

    ZHANG Yan1, ZHANG Qiuying2**, LI Fadong3,4, ZHANG Xin1, BI Zhilei1, ZHANG Qiang1

    (1. Shaanxi Key Laboratory of Earth Surface System and Environmental Carrying Capacity, College of Urban and Environmental Sciences, Northwest University, Xi’an 710127, China; 2. Chinese Research Academy of Environmental Sciences, Beijing 100012, China; 3. Key Laboratory of Ecosystem Network Observation and Modeling, Institute of Geographical Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China; 4. College of Resources and Environment, University of Chinese Academy of Sciences, Beijing 100190, China)

    Nitrate (NO3-) pollution in groundwater has become a serious environmental problem across the world. It is very important to determine the sources of nitrogen contamination in order to prevent and control NO3-pollution in groundwater. This is because the intake of polluted water can increase health risk of methemoglobinemia and cancer in both aquatic lives and humans. There has been an increasing trend in NO3-pollution in groundwater in the Lower Yellow River Irrigation Districts. Once groundwater is polluted by NO3-, recovery efforts can be very daunting. The effective control and management of NO3-pollution require accurate identification of the actual sources of pollution.In this paper, the sources of NO3-in groundwater in the Lower Yellow River Irrigation District (Panzhuang Irrigation District) were identified using stable isotopes (δ15N and δ18O) and the Bayesian model. The results showed that the range of NO3-concentrations in groundwater in the study area was 0.1–197.0 mg·L?1, with a mean of 34.2 mg·L?1. About 10% of the groundwater samples had NO3-concentration in excess of the maximal standard of nitrate level in drinking water in China (90 mg·L?1). Samples were divided into three depths, including 0-30 m (shallow layer), 30-60 m (middle layer) and >60 m (deep layer). The average NO3-concentrations in shallow groundwater layer, middle layer and deep layer were 25.9 mg·L?1, 39.7 mg·L?1and 20.1 mg·L?1, respectively. There were high NO3-concentrations in groundwater across Ningjin County, Wucheng County, Pingyuan County and Yucheng City. The composition of δ15N was in the range of 0.72‰–23.93‰, with an average of 11.62‰. That of δ18O was 0.49‰–22.50‰, with an average of 8.46‰. The values of δ15N and δ18O indicated that NO3-in groundwater in the study area mainly originated from chemical fertilizers, manure and sewage. The contributions of the four sources of NO3-(precipitation, chemical fertilizer, soil, manure and sewage) were quantified and estimated using the Bayesian model. The results showed that manure and sewage contributed the most to the overall NO3-level, with a mean NO3-contributionratio of 56.2%. Chemical fertilizer was the second contributor, with a mean NO3-contribution ratio of 19.3%. The mean NO3-contribution ratio of precipitation and soil was6.2% and 12.3%, respectively. After identification of NO3-pollution levels and sources, measures were required to reduce NO3-pollution in groundwater. Based on this study, the necessary measures included the construction of sewage pipeline and improving the utilization rate of chemical fertilizers in order to reduce NO3-pollution and improve water quality.

    Nitrogen pollution; Groundwater contamination; Nitrate; Stable isotopes; Bayesian model

    , E-mail: zhangqy@craes.org.cn

    Sep. 29, 2018;

    Nov. 27, 2018

    X523

    A

    2096-6237(2019)03-0484-10

    10.13930/j.cnki.cjea.180887

    張妍, 張秋英, 李發(fā)東, 張鑫, 畢直磊, 張強. 基于穩(wěn)定同位素和貝葉斯模型的引黃灌區(qū)地下水硝酸鹽污染源解析[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(中英文), 2019, 27(3): 484-493

    ZHANG Y, ZHANG Q Y, LI F D, ZHANG X, BI Z L, ZHANG Q. Source identification of nitrate contamination of groundwater in Yellow River Irrigation Districts using stable isotopes and Bayesian model[J]. Chinese Journal of Eco-Agriculture, 2019, 27(3): 484-493

    * 國家自然科學(xué)基金項目(41601017, 41771292, 41271047)、陜西省自然科學(xué)基礎(chǔ)研究項目(2017JQ4001)、中國博士后基金項目(2015M572591)和陜西省教育廳科研計劃項目(17JK0771)資助

    張秋英, 主要研究方向為流域水環(huán)境。E-mail: zhangqy@craes.org.cn

    張妍, 主要研究方向為農(nóng)業(yè)生態(tài)水循環(huán)與環(huán)境氮污染。E-mail: yanz@nwu.edu.cn

    2018-09-29

    2018-11-27

    * The study was supported by the National Natural Science Foundation of China (41601017, 41771292, 41271047), the Natural Science Foundation Research Project of Shaanxi Province (2017JQ4001), the Postdoctoral Science Foundation of China (2015M572591), and the Scientific Research Plan Projects of Shaanxi Education Department (17JK0771).

    猜你喜歡
    污染
    河流被污染了嗎?
    什么是污染?
    什么是污染?
    堅決打好污染防治攻堅戰(zhàn)
    堅決打好污染防治攻堅戰(zhàn)
    可以喝的塑料:污染解決之道?
    飲用水污染 誰之過?
    食品界(2016年4期)2016-02-27 07:36:15
    對抗塵污染,遠離“霾”伏
    都市麗人(2015年5期)2015-03-20 13:33:49
    P265GH低合金鋼在模擬污染大氣環(huán)境中的腐蝕行為
    污染防治
    江蘇年鑒(2014年0期)2014-03-11 17:09:46
    精品国内亚洲2022精品成人| 夜夜爽天天搞| 午夜a级毛片| 久久久久久九九精品二区国产| 亚洲色图av天堂| 在线国产一区二区在线| 欧洲精品卡2卡3卡4卡5卡区| 在线观看免费午夜福利视频| 午夜免费激情av| 午夜亚洲福利在线播放| 久久久精品欧美日韩精品| 宅男免费午夜| 悠悠久久av| 99视频精品全部免费 在线 | 亚洲最大成人中文| 啦啦啦观看免费观看视频高清| 黑人欧美特级aaaaaa片| www.精华液| 久久草成人影院| 国产麻豆成人av免费视频| 毛片女人毛片| 久久精品91无色码中文字幕| 最近最新中文字幕大全免费视频| 国产一区二区激情短视频| 国产私拍福利视频在线观看| 男女下面进入的视频免费午夜| 久久国产乱子伦精品免费另类| 午夜影院日韩av| 999久久久精品免费观看国产| 女同久久另类99精品国产91| 亚洲午夜精品一区,二区,三区| 欧美乱色亚洲激情| 国产激情欧美一区二区| 日韩大尺度精品在线看网址| 91老司机精品| 九色国产91popny在线| 最近在线观看免费完整版| 欧美一区二区国产精品久久精品| 亚洲熟妇中文字幕五十中出| 久久精品aⅴ一区二区三区四区| 亚洲国产看品久久| 午夜福利欧美成人| ponron亚洲| 久久欧美精品欧美久久欧美| 久久久久亚洲av毛片大全| 亚洲 国产 在线| 日韩欧美 国产精品| 1000部很黄的大片| 日韩欧美 国产精品| 午夜视频精品福利| 欧美一区二区精品小视频在线| 黄色 视频免费看| 亚洲欧美精品综合一区二区三区| h日本视频在线播放| 小说图片视频综合网站| 精品不卡国产一区二区三区| 亚洲午夜精品一区,二区,三区| 婷婷亚洲欧美| 欧美日韩瑟瑟在线播放| 日韩精品中文字幕看吧| 亚洲午夜精品一区,二区,三区| www日本在线高清视频| 国产毛片a区久久久久| h日本视频在线播放| 亚洲熟女毛片儿| 日本黄大片高清| 亚洲av成人不卡在线观看播放网| 18禁观看日本| 亚洲美女黄片视频| 亚洲欧美日韩无卡精品| 国产欧美日韩精品亚洲av| 亚洲成人中文字幕在线播放| 中文字幕av在线有码专区| 国产在线精品亚洲第一网站| 丰满人妻熟妇乱又伦精品不卡| 中文在线观看免费www的网站| 亚洲在线自拍视频| 成人精品一区二区免费| 麻豆一二三区av精品| 色综合婷婷激情| 婷婷亚洲欧美| www日本在线高清视频| 99久久99久久久精品蜜桃| 成人精品一区二区免费| 亚洲专区字幕在线| 亚洲一区二区三区不卡视频| 亚洲自偷自拍图片 自拍| 曰老女人黄片| 麻豆成人av在线观看| 成年版毛片免费区| 国产精品99久久99久久久不卡| 亚洲国产精品合色在线| 又黄又爽又免费观看的视频| 女同久久另类99精品国产91| 国产精品亚洲一级av第二区| 国产一区二区三区视频了| 午夜视频精品福利| 999久久久精品免费观看国产| 99久久久亚洲精品蜜臀av| 国产爱豆传媒在线观看| 久久香蕉精品热| 全区人妻精品视频| 每晚都被弄得嗷嗷叫到高潮| 男人舔女人的私密视频| 人人妻,人人澡人人爽秒播| 又粗又爽又猛毛片免费看| 午夜免费成人在线视频| 巨乳人妻的诱惑在线观看| 国内揄拍国产精品人妻在线| 国产私拍福利视频在线观看| 日韩三级视频一区二区三区| 19禁男女啪啪无遮挡网站| 久久久国产精品麻豆| 日韩欧美免费精品| 69av精品久久久久久| 夜夜爽天天搞| 久久午夜亚洲精品久久| xxxwww97欧美| 亚洲国产中文字幕在线视频| 日韩欧美国产在线观看| 可以在线观看毛片的网站| 搡老熟女国产l中国老女人| 99国产精品一区二区三区| 欧美日韩精品网址| 国产极品精品免费视频能看的| 欧美日韩黄片免| 欧美+亚洲+日韩+国产| 在线视频色国产色| 日韩成人在线观看一区二区三区| 午夜福利欧美成人| www.999成人在线观看| 综合色av麻豆| 国产精品自产拍在线观看55亚洲| 久久精品夜夜夜夜夜久久蜜豆| 在线十欧美十亚洲十日本专区| 亚洲精品一卡2卡三卡4卡5卡| 人人妻,人人澡人人爽秒播| 制服丝袜大香蕉在线| 18禁观看日本| 久久国产精品影院| 日韩欧美 国产精品| 成人三级黄色视频| 少妇丰满av| 国产精品免费一区二区三区在线| 亚洲国产欧洲综合997久久,| 91av网站免费观看| 免费看十八禁软件| 国产真实乱freesex| 国产亚洲欧美98| 色av中文字幕| 制服人妻中文乱码| 一级毛片高清免费大全| 不卡av一区二区三区| 亚洲熟女毛片儿| 精品人妻1区二区| 巨乳人妻的诱惑在线观看| 欧美黄色片欧美黄色片| 久久久久久国产a免费观看| 性色av乱码一区二区三区2| 日本一二三区视频观看| 久久久久性生活片| 91麻豆av在线| 亚洲在线观看片| 91在线观看av| 天天一区二区日本电影三级| 亚洲熟女毛片儿| 久久久久久人人人人人| 精华霜和精华液先用哪个| 宅男免费午夜| 最近最新免费中文字幕在线| 国产成人av教育| 欧美成人免费av一区二区三区| 啦啦啦观看免费观看视频高清| 好看av亚洲va欧美ⅴa在| 免费无遮挡裸体视频| 日韩免费av在线播放| 日韩成人在线观看一区二区三区| 国产私拍福利视频在线观看| 亚洲av成人av| 俄罗斯特黄特色一大片| 99国产精品一区二区蜜桃av| 男人和女人高潮做爰伦理| 免费在线观看亚洲国产| 欧美最黄视频在线播放免费| 精品福利观看| 国产99白浆流出| 久久久国产欧美日韩av| 国产伦精品一区二区三区视频9 | 成年版毛片免费区| 18禁观看日本| 精品国产美女av久久久久小说| 九九热线精品视视频播放| 天天躁日日操中文字幕| 国产美女午夜福利| 日本三级黄在线观看| 91在线观看av| 身体一侧抽搐| 亚洲国产看品久久| 久久久久久大精品| 日本一二三区视频观看| 国产单亲对白刺激| 日本在线视频免费播放| 怎么达到女性高潮| 国产精品国产高清国产av| 亚洲专区中文字幕在线| 精品国产美女av久久久久小说| 久久中文字幕人妻熟女| 最近在线观看免费完整版| 美女扒开内裤让男人捅视频| 人妻丰满熟妇av一区二区三区| 成人国产综合亚洲| 国内少妇人妻偷人精品xxx网站 | e午夜精品久久久久久久| 久久国产精品人妻蜜桃| 国产精品国产高清国产av| 欧美又色又爽又黄视频| 婷婷六月久久综合丁香| 综合色av麻豆| 久久午夜亚洲精品久久| 国产成人啪精品午夜网站| 一级毛片高清免费大全| 国产乱人伦免费视频| 午夜福利在线在线| 亚洲av成人av| 黄片大片在线免费观看| www国产在线视频色| 国产欧美日韩精品一区二区| 亚洲色图av天堂| 亚洲成人精品中文字幕电影| 欧美一区二区国产精品久久精品| 日本 欧美在线| 久久久久性生活片| 丁香六月欧美| 国语自产精品视频在线第100页| 国产男靠女视频免费网站| 国产欧美日韩精品亚洲av| 这个男人来自地球电影免费观看| 国产v大片淫在线免费观看| 欧美日本视频| 人妻夜夜爽99麻豆av| 亚洲,欧美精品.| av视频在线观看入口| 国产伦在线观看视频一区| 国产精品1区2区在线观看.| 欧美+亚洲+日韩+国产| 男人舔女人的私密视频| 日韩 欧美 亚洲 中文字幕| 热99在线观看视频| 久久天堂一区二区三区四区| 精品免费久久久久久久清纯| 免费在线观看成人毛片| 国产成人av激情在线播放| 精品福利观看| 中文字幕人成人乱码亚洲影| 淫秽高清视频在线观看| tocl精华| 久久人妻av系列| 黑人巨大精品欧美一区二区mp4| 搡老岳熟女国产| 九色成人免费人妻av| 香蕉国产在线看| www.熟女人妻精品国产| 日本熟妇午夜| 亚洲精品在线观看二区| 九色成人免费人妻av| 欧美日韩福利视频一区二区| 精品久久久久久久久久久久久| 欧美绝顶高潮抽搐喷水| 成人国产一区最新在线观看| 亚洲精品国产精品久久久不卡| 久久午夜亚洲精品久久| 国产又黄又爽又无遮挡在线| 激情在线观看视频在线高清| 国产高潮美女av| 十八禁网站免费在线| 日本一二三区视频观看| 色尼玛亚洲综合影院| 婷婷精品国产亚洲av| 美女大奶头视频| 综合色av麻豆| 女警被强在线播放| 怎么达到女性高潮| 国产亚洲av高清不卡| 男女做爰动态图高潮gif福利片| 两个人视频免费观看高清| 亚洲成a人片在线一区二区| 亚洲人成网站高清观看| 欧美日韩综合久久久久久 | 真人做人爱边吃奶动态| 成年女人看的毛片在线观看| 九九久久精品国产亚洲av麻豆 | 在线免费观看不下载黄p国产 | 久久人妻av系列| 激情在线观看视频在线高清| 久久久久久人人人人人| 国产私拍福利视频在线观看| 精品欧美国产一区二区三| 国产综合懂色| 国产精品99久久久久久久久| 国产精品久久久久久亚洲av鲁大| 国产欧美日韩精品亚洲av| 嫩草影院入口| 精品久久蜜臀av无| 午夜日韩欧美国产| 一进一出抽搐gif免费好疼| 国产精品久久久久久精品电影| 偷拍熟女少妇极品色| aaaaa片日本免费| 国产激情欧美一区二区| 午夜视频精品福利| 岛国视频午夜一区免费看| 18禁裸乳无遮挡免费网站照片| 国产成人啪精品午夜网站| 黄色成人免费大全| 搡老岳熟女国产| 欧美成人性av电影在线观看| 精品一区二区三区四区五区乱码| 麻豆一二三区av精品| 久久精品人妻少妇| 女生性感内裤真人,穿戴方法视频| 美女cb高潮喷水在线观看 | 18禁观看日本| 久久香蕉国产精品| 99在线视频只有这里精品首页| 午夜久久久久精精品| 国产乱人伦免费视频| 国产成人系列免费观看| 香蕉丝袜av| 成熟少妇高潮喷水视频| 在线观看舔阴道视频| 精品电影一区二区在线| 他把我摸到了高潮在线观看| 天天躁日日操中文字幕| 欧美日韩亚洲国产一区二区在线观看| 欧美乱码精品一区二区三区| 久久精品91蜜桃| 在线永久观看黄色视频| 国产精品久久视频播放| 日韩欧美国产在线观看| 精品国产亚洲在线| 亚洲熟妇中文字幕五十中出| 最近最新免费中文字幕在线| 最新美女视频免费是黄的| 久久精品91蜜桃| 久久中文看片网| 亚洲成av人片免费观看| 啦啦啦韩国在线观看视频| 国产高清视频在线播放一区| 一二三四社区在线视频社区8| 国产淫片久久久久久久久 | 999精品在线视频| 国产黄a三级三级三级人| 欧美黄色淫秽网站| 国产黄片美女视频| 久久中文看片网| 亚洲成人久久爱视频| 欧美乱码精品一区二区三区| 啦啦啦免费观看视频1| 天天躁日日操中文字幕| 中文在线观看免费www的网站| 亚洲国产高清在线一区二区三| 久久这里只有精品中国| 国产黄色小视频在线观看| 成熟少妇高潮喷水视频| 三级男女做爰猛烈吃奶摸视频| 99热只有精品国产| 18禁黄网站禁片午夜丰满| 国产亚洲精品一区二区www| 久久国产乱子伦精品免费另类| 欧美在线一区亚洲| 丰满人妻熟妇乱又伦精品不卡| 不卡一级毛片| 久久久久久久午夜电影| 亚洲人与动物交配视频| 国产午夜精品论理片| avwww免费| 母亲3免费完整高清在线观看| 亚洲国产中文字幕在线视频| 精品福利观看| 国产精品免费一区二区三区在线| 真人一进一出gif抽搐免费| 一级毛片精品| 村上凉子中文字幕在线| 两个人视频免费观看高清| 亚洲精品乱码久久久v下载方式 | 国产毛片a区久久久久| 1024香蕉在线观看| bbb黄色大片| 亚洲av成人一区二区三| 97人妻精品一区二区三区麻豆| 午夜福利在线观看吧| 人妻丰满熟妇av一区二区三区| 97碰自拍视频| 少妇裸体淫交视频免费看高清| 成人18禁在线播放| 精品国产超薄肉色丝袜足j| 欧美一区二区国产精品久久精品| 久久精品国产99精品国产亚洲性色| 免费av不卡在线播放| 变态另类丝袜制服| 国产单亲对白刺激| 日本成人三级电影网站| 亚洲中文av在线| 亚洲精品456在线播放app | 国产精品1区2区在线观看.| 亚洲熟妇中文字幕五十中出| 久久草成人影院| 精品一区二区三区视频在线观看免费| 亚洲av日韩精品久久久久久密| 性欧美人与动物交配| 国产精品久久久久久久电影 | 女生性感内裤真人,穿戴方法视频| 久久这里只有精品19| 久久久久国产精品人妻aⅴ院| 久久精品国产亚洲av香蕉五月| 国产欧美日韩一区二区精品| 好看av亚洲va欧美ⅴa在| 成人国产一区最新在线观看| 88av欧美| a在线观看视频网站| 亚洲av电影不卡..在线观看| 成人精品一区二区免费| 999久久久精品免费观看国产| 国产高清激情床上av| 一个人观看的视频www高清免费观看 | 成人鲁丝片一二三区免费| 精品国内亚洲2022精品成人| 国产一区二区三区视频了| 久久天堂一区二区三区四区| 一个人免费在线观看电影 | 国产av麻豆久久久久久久| 国产视频一区二区在线看| 国产精品香港三级国产av潘金莲| 久久久久久久午夜电影| 一级作爱视频免费观看| 国产探花在线观看一区二区| 动漫黄色视频在线观看| 国产精品av视频在线免费观看| 91在线精品国自产拍蜜月 | 久久久精品大字幕| 亚洲乱码一区二区免费版| 美女cb高潮喷水在线观看 | 国产三级黄色录像| 校园春色视频在线观看| 变态另类成人亚洲欧美熟女| 久久久久久大精品| 久久伊人香网站| 婷婷精品国产亚洲av在线| 禁无遮挡网站| 免费在线观看视频国产中文字幕亚洲| 99热只有精品国产| 女人被狂操c到高潮| 成人av一区二区三区在线看| 亚洲aⅴ乱码一区二区在线播放| 国产在线精品亚洲第一网站| 国产成人影院久久av| 亚洲一区高清亚洲精品| 九九在线视频观看精品| or卡值多少钱| xxx96com| 国产一区二区在线观看日韩 | 国产精品一区二区三区四区久久| 两人在一起打扑克的视频| 免费av不卡在线播放| 国产精品自产拍在线观看55亚洲| 国产综合懂色| 丁香欧美五月| 欧美日韩福利视频一区二区| 一二三四在线观看免费中文在| 观看美女的网站| 美女被艹到高潮喷水动态| 久久久成人免费电影| 桃色一区二区三区在线观看| 欧美三级亚洲精品| 在线观看美女被高潮喷水网站 | 欧美+亚洲+日韩+国产| 日韩欧美三级三区| 日本 欧美在线| 成人亚洲精品av一区二区| 欧美精品啪啪一区二区三区| 日本与韩国留学比较| 我要搜黄色片| 久久久国产成人精品二区| 一个人看视频在线观看www免费 | 亚洲av片天天在线观看| 国产精品影院久久| av在线天堂中文字幕| 欧美黄色淫秽网站| 久久久久国产精品人妻aⅴ院| 九九久久精品国产亚洲av麻豆 | 男女之事视频高清在线观看| 亚洲成a人片在线一区二区| 真人做人爱边吃奶动态| 日韩精品青青久久久久久| 九九热线精品视视频播放| 精品熟女少妇八av免费久了| 日日摸夜夜添夜夜添小说| 国产精品久久久久久亚洲av鲁大| 亚洲性夜色夜夜综合| 欧美3d第一页| 757午夜福利合集在线观看| 热99在线观看视频| 日韩高清综合在线| 日本一本二区三区精品| 免费看a级黄色片| 精品欧美国产一区二区三| 国产精品99久久99久久久不卡| 国内毛片毛片毛片毛片毛片| 国产激情欧美一区二区| 国产精品,欧美在线| 99热这里只有精品一区 | 久久精品影院6| 国产久久久一区二区三区| 99国产极品粉嫩在线观看| 后天国语完整版免费观看| 麻豆成人av在线观看| 亚洲黑人精品在线| 老司机深夜福利视频在线观看| 亚洲av电影在线进入| 国产精品电影一区二区三区| 国产精品一及| 色av中文字幕| 日韩欧美精品v在线| aaaaa片日本免费| 午夜福利高清视频| 亚洲成av人片免费观看| 一个人看的www免费观看视频| 精品国产乱码久久久久久男人| 国产av麻豆久久久久久久| 99热只有精品国产| 日本在线视频免费播放| 97人妻精品一区二区三区麻豆| 国产视频内射| 亚洲av成人不卡在线观看播放网| 日韩高清综合在线| 国产精品99久久99久久久不卡| 亚洲aⅴ乱码一区二区在线播放| 亚洲午夜理论影院| 长腿黑丝高跟| 久久热在线av| 国产精品久久久人人做人人爽| 久久久久国产一级毛片高清牌| 亚洲精品一区av在线观看| 成人av在线播放网站| 亚洲中文字幕一区二区三区有码在线看 | 在线观看午夜福利视频| 国语自产精品视频在线第100页| 搡老妇女老女人老熟妇| netflix在线观看网站| 国产精品一区二区三区四区久久| 天堂影院成人在线观看| 久久久成人免费电影| 色av中文字幕| 亚洲国产精品999在线| 女人被狂操c到高潮| 国产野战对白在线观看| 免费在线观看成人毛片| 国产高清视频在线播放一区| 国产精品 国内视频| 99国产精品99久久久久| 亚洲欧美日韩无卡精品| 听说在线观看完整版免费高清| 久久久国产成人精品二区| 一个人看的www免费观看视频| 国产私拍福利视频在线观看| 中文字幕精品亚洲无线码一区| 成年版毛片免费区| 日韩三级视频一区二区三区| 中出人妻视频一区二区| 精品国产超薄肉色丝袜足j| 欧美乱妇无乱码| 国产99白浆流出| 午夜福利在线观看吧| 国产精品亚洲一级av第二区| 国产成人影院久久av| 精品国产亚洲在线| 1024手机看黄色片| 12—13女人毛片做爰片一| 国产精品美女特级片免费视频播放器 | 午夜影院日韩av| 欧美丝袜亚洲另类 | 精品人妻1区二区| 亚洲成a人片在线一区二区| 久久人人精品亚洲av| 亚洲一区高清亚洲精品| 成人三级黄色视频| 特大巨黑吊av在线直播| 男女下面进入的视频免费午夜| 美女扒开内裤让男人捅视频| 草草在线视频免费看| 国产又黄又爽又无遮挡在线| 美女高潮喷水抽搐中文字幕| 高潮久久久久久久久久久不卡| 国产高清视频在线观看网站| 中文字幕久久专区| 日韩三级视频一区二区三区| 看黄色毛片网站| 制服丝袜大香蕉在线| 又黄又爽又免费观看的视频| 91久久精品国产一区二区成人 | 亚洲第一欧美日韩一区二区三区| www日本黄色视频网| 亚洲精品乱码久久久v下载方式 | 亚洲激情在线av| 国产精品一区二区免费欧美| 亚洲精品美女久久av网站| 亚洲中文日韩欧美视频| 亚洲九九香蕉| 熟妇人妻久久中文字幕3abv| 91久久精品国产一区二区成人 | 午夜精品一区二区三区免费看| 天堂√8在线中文| 免费在线观看日本一区| 亚洲av电影在线进入| 少妇熟女aⅴ在线视频|