• <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ā)東, 張 鑫, 畢直磊, 張 強(qiáng)

    ?

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

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

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

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

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

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

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

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

    1 材料與方法

    1.1 研究區(qū)概況

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

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

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

    潘莊引黃灌區(qū)農(nóng)墾歷史悠久, 土地利用率較高。研究區(qū)種植作物主要是夏玉米(L.)和冬小麥(L.)。主要使用的化肥種類有碳銨、尿素和緩釋肥等[13]。長(zhǎng)期以來(lái)存在施用單元肥料和過(guò)量施肥等問(wèn)題, 導(dǎo)致土壤的養(yǎng)分結(jié)構(gòu)發(fā)生較大變化, 養(yǎng)分比例失調(diào), 致使土壤肥力狀況仍然滿足不了農(nóng)作物產(chǎn)量日益增長(zhǎng)的需要, 各種環(huán)境問(wèn)題如NO3-淋失等問(wèn)題凸顯。除了農(nóng)業(yè)區(qū)氮肥施用是水體NO3-的潛在污染來(lái)源, 潘莊引黃灌區(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è)計(jì)技術(shù)規(guī)定》(GB12997—97)中規(guī)定的水質(zhì)量控制、質(zhì)量表征、污染物鑒別采樣方案的原則, 設(shè)置了水樣點(diǎn)網(wǎng)格布設(shè)。于2013年9月采集40個(gè)地下水樣品(圖1), 采集的地下水包括36個(gè)淺層地下水和4個(gè)中層承壓微咸水。對(duì)各個(gè)取樣點(diǎn)進(jìn)行GPS定位, 并在現(xiàn)場(chǎng)測(cè)定pH、EC、溶解氧(DO)和氧化還原電位(ORP)。

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

    1.3 貝葉斯混合模型

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

    式中:X表示混合物同位素的值;P表示來(lái)源的比例;S表示第個(gè)來(lái)源的第種同位素的值, 服從均值為方差的正態(tài)分布;c表示第個(gè)來(lái)源的同位素的分餾系數(shù), 服從均值為方差為的正態(tài)分布;是殘余誤差, 表示其他各個(gè)混合物間無(wú)法量化的方差, 其均值和標(biāo)準(zhǔn)差在通常情況下均為0。本文采用貝葉斯混合模型定量計(jì)算了4種潛在來(lái)源(大氣降水、土壤氮、化肥、污水和糞便)對(duì)地下水中NO3-的貢獻(xiàn)比率, 這4種來(lái)源的同位素值參照文獻(xiàn)[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-次之; 陽(yáng)離子以Na+離子為主, Ca2+離子次之。陰離子中HCO3-離子含量294.5~ 1 259.4 mg?L-1; 陽(yáng)離子中Na+離子含量為48.7~627.2 mg?L-1。水組成反映了碳酸鹽巖溶解的影響。華北平原土壤鹽堿化嚴(yán)重, 鹽漬土中含有大量的NaHCO3和Na2HCO3, 隨灌溉和降水從土壤中溶解進(jìn)入地表水和地下水[27]。DO含量范圍為0.01~6.70 mg?L-1, 反映了補(bǔ)給條件和反硝化發(fā)生的可能[28]。NO3-含量介于0.1~197.0 mg?L-1, 平均值為34.2 mg?L-1。從變異系數(shù)來(lái)看, 離子中NO3-變異系數(shù)最大, 達(dá)139.4%, 大于100%, 表明NO3-離子空間變異大, 可能是因?yàn)榈叵滤鲃?dòng)過(guò)程中, 不同的水文地質(zhì)條件以及人類活動(dòng)如農(nóng)田施肥、污水排放等導(dǎo)致地下水水化學(xué)特征的空間變異增大[29]。

    由圖2可以看出, EC、HCO3-、SO42-、Cl-、Ca2+、Na+、Mg2+、K+隨著地下水深的增加而出現(xiàn)先增加后減少的趨勢(shì), 最大值出現(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-平均含量最高, 達(dá)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%的樣點(diǎn)超過(guò)《生活飲用水衛(wèi)生標(biāo)準(zhǔn)》(GB5749—2006)中NO3-的上限值90 mg?L-1。總體上看, 本研究區(qū)西北部地下水NO3-含量很低, 普遍小于45 mg?L-1; 而東北部NO3-含量較高, 其中, 寧津縣、武城縣、平原縣和陵縣有大片區(qū)域NO3-含量超過(guò)45 mg?L-1(圖3)。

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

    2.3 地下水NO3-氮氧同位素組成及各種來(lái)源貢獻(xiàn)率

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

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

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

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

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

    圖6 潘莊引黃灌區(qū)地下水NO3-潛在來(lái)源對(duì)地下水NO3-的貢獻(xiàn)率

    3 討論與結(jié)論

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

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

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

    本研究中, 潘莊灌區(qū)地下水NO3-濃度介于0.1~197.01 mg?L-1, 平均濃度為34.2 mg?L-1。有10%的地下水樣品超過(guò)了《生活飲用水衛(wèi)生標(biāo)準(zhǔn)》中規(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‰。貝葉斯混合模型計(jì)算結(jié)果表明, 地下水中約有56.2%的NO3-來(lái)自糞便和污水, 19.3%的NO3-來(lái)自于農(nóng)業(yè)化肥。為保護(hù)和改善研究區(qū)地下水水質(zhì), 控制地下水NO3-的進(jìn)一步增長(zhǎng), 需要從污染源入手采取如下措施: 1)有關(guān)部門(mén)應(yīng)以現(xiàn)在最嚴(yán)格的水資源管理制度和限制納污紅線考核制度為契機(jī), 嚴(yán)格監(jiān)管各排污口污水的不達(dá)標(biāo)排放。同時(shí), 加強(qiáng)污水管道建設(shè), 進(jìn)一步提高污水處理率。2)加強(qiáng)畜禽糞便的管理, 糞便堆放處必須有防滲措施, 并妥善處理。同時(shí)鼓勵(lì)合理地將其轉(zhuǎn)化為有機(jī)肥。3)提高化肥的利用效率, 加強(qiáng)引黃灌區(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] 劉貫群, 周書(shū)玉, 黃修東, 等. 多種方法識(shí)別青島大沽河平原區(qū)地下水硝酸鹽污染來(lái)源[J]. 環(huán)境科學(xué)學(xué)報(bào), 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)境影響研究——以山東省桓臺(tái)縣為例[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 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)查與評(píng)價(jià)[J]. 中國(guó)農(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)定同位素追蹤石家莊市地下水中的硝酸鹽來(lái)源[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] 魯垠濤, 劉芳, 姚宏, 等. 北京密云水庫(kù)小流域地下水硝酸鹽污染來(lái)源示蹤[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)濟(jì), 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]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 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] 劉桂儀. 魯北平原深層地下水開(kāi)發(fā)與環(huán)境問(wè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] 賀國(guó)平, 劉培斌, 慕星, 等. 永定河沖洪積扇地下水中硝酸鹽來(lái)源的同位素識(shí)別[J]. 水利學(xué)報(bào), 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] 張琳, 高翔. 我國(guó)廢水中主要污染物排放情況分析[J]. 統(tǒng)計(jì)學(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]. 中國(guó)給水排水, 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)場(chǎng)氮素流動(dòng)特征及利用效率——以河北津龍循環(huán)農(nóng)業(yè)園區(qū)為例[J]. 應(yīng)用生態(tài)學(xué)報(bào), 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] 王仕琴, 鄭文波, 孔曉樂(lè). 華北農(nóng)區(qū)淺層地下水硝酸鹽分布特征及其空間差異性[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 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)計(jì)局. 山東統(tǒng)計(jì)年鑒1993—2013[M]. 北京: 中國(guó)統(tǒng)計(jì)出版社, 1993—2013 Shandong Provincial Bureau of Statistics. Statistical Yearbook of Shandong 1993-2013[M]. Beijing: China Statistics Press, 1993–2013

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

    [38] 楊靜, 肖天昀, 李海波, 等. 江漢平原地下水中硝酸鹽的分布及影響因素[J]. 中國(guó)環(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ā)東, 張?chǎng)? 畢直磊, 張強(qiáng). 基于穩(wěn)定同位素和貝葉斯模型的引黃灌區(qū)地下水硝酸鹽污染源解析[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)(中英文), 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

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

    張秋英, 主要研究方向?yàn)榱饔蛩h(huán)境。E-mail: zhangqy@craes.org.cn

    張妍, 主要研究方向?yàn)檗r(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).

    猜你喜歡
    污染
    河流被污染了嗎?
    什么是污染?
    什么是污染?
    堅(jiān)決打好污染防治攻堅(jiān)戰(zhàn)
    堅(jiān)決打好污染防治攻堅(jiān)戰(zhàn)
    可以喝的塑料:污染解決之道?
    飲用水污染 誰(shuí)之過(guò)?
    食品界(2016年4期)2016-02-27 07:36:15
    對(duì)抗塵污染,遠(yuǎn)離“霾”伏
    都市麗人(2015年5期)2015-03-20 13:33:49
    P265GH低合金鋼在模擬污染大氣環(huán)境中的腐蝕行為
    污染防治
    江蘇年鑒(2014年0期)2014-03-11 17:09:46
    av视频在线观看入口| 联通29元200g的流量卡| 中文亚洲av片在线观看爽| 91av网一区二区| 中文字幕熟女人妻在线| 午夜爱爱视频在线播放| 日本一本二区三区精品| 日本-黄色视频高清免费观看| 国内揄拍国产精品人妻在线| 99国产精品一区二区蜜桃av| 亚洲精品乱码久久久v下载方式| 欧美日韩国产亚洲二区| 久久精品人妻少妇| 国产成年人精品一区二区| 欧美成人一区二区免费高清观看| 久久精品国产亚洲av天美| 成人高潮视频无遮挡免费网站| 国产亚洲5aaaaa淫片| 麻豆国产97在线/欧美| 男插女下体视频免费在线播放| 蜜桃久久精品国产亚洲av| 人妻制服诱惑在线中文字幕| 一区福利在线观看| 97热精品久久久久久| 国产91av在线免费观看| 免费av不卡在线播放| 在线观看一区二区三区| 亚洲精品影视一区二区三区av| 欧美成人免费av一区二区三区| kizo精华| 少妇人妻一区二区三区视频| 亚洲aⅴ乱码一区二区在线播放| 免费看日本二区| 中国美白少妇内射xxxbb| 日韩一本色道免费dvd| 亚洲av成人av| 99久国产av精品国产电影| 内地一区二区视频在线| 欧美日韩综合久久久久久| 草草在线视频免费看| 99国产精品一区二区蜜桃av| 精品人妻一区二区三区麻豆| 寂寞人妻少妇视频99o| 免费电影在线观看免费观看| 男人的好看免费观看在线视频| 又爽又黄a免费视频| 亚洲av二区三区四区| 午夜a级毛片| 岛国在线免费视频观看| 国内久久婷婷六月综合欲色啪| 少妇高潮的动态图| 国产精品,欧美在线| 最后的刺客免费高清国语| 别揉我奶头 嗯啊视频| 欧美不卡视频在线免费观看| 久久久久久久久久黄片| 美女cb高潮喷水在线观看| 乱码一卡2卡4卡精品| 精品不卡国产一区二区三区| 又粗又硬又长又爽又黄的视频 | 欧美变态另类bdsm刘玥| 91aial.com中文字幕在线观看| 18禁在线播放成人免费| 欧美成人免费av一区二区三区| 亚洲最大成人中文| 亚洲美女视频黄频| 一区二区三区免费毛片| 国产精品麻豆人妻色哟哟久久 | 少妇高潮的动态图| 亚洲av.av天堂| 国模一区二区三区四区视频| 黄色欧美视频在线观看| av天堂在线播放| 中国国产av一级| 青春草国产在线视频 | 国产成人福利小说| 秋霞在线观看毛片| 国产毛片a区久久久久| 亚洲成人中文字幕在线播放| 欧美xxxx黑人xx丫x性爽| 日韩欧美国产在线观看| 欧美精品一区二区大全| 一级二级三级毛片免费看| 美女高潮的动态| 亚洲在久久综合| 精品久久久久久久久亚洲| 欧美zozozo另类| 成年av动漫网址| .国产精品久久| 亚洲av男天堂| 久久久国产成人免费| 免费人成在线观看视频色| 国产女主播在线喷水免费视频网站 | av在线老鸭窝| 91午夜精品亚洲一区二区三区| 激情 狠狠 欧美| 夜夜夜夜夜久久久久| 麻豆成人av视频| 五月伊人婷婷丁香| 六月丁香七月| 啦啦啦韩国在线观看视频| 全区人妻精品视频| 亚洲欧洲日产国产| 久久精品国产99精品国产亚洲性色| 成人永久免费在线观看视频| 免费av毛片视频| 日本一二三区视频观看| 亚洲欧美成人精品一区二区| 久久这里只有精品中国| 亚洲国产色片| 欧美日韩国产亚洲二区| 久久精品国产清高在天天线| 久久午夜亚洲精品久久| 国产午夜精品论理片| www日本黄色视频网| 国产精品久久久久久av不卡| 大又大粗又爽又黄少妇毛片口| 久久久久久九九精品二区国产| 成人特级av手机在线观看| 国产中年淑女户外野战色| 麻豆成人av视频| 日韩精品有码人妻一区| 99在线视频只有这里精品首页| 免费一级毛片在线播放高清视频| av国产免费在线观看| 亚洲精品影视一区二区三区av| 久久人妻av系列| 亚洲成人久久性| a级毛片a级免费在线| 麻豆av噜噜一区二区三区| 黄色配什么色好看| 亚洲国产欧洲综合997久久,| 性欧美人与动物交配| 色综合色国产| 国产黄a三级三级三级人| 少妇人妻精品综合一区二区 | 国国产精品蜜臀av免费| 午夜精品国产一区二区电影 | 色播亚洲综合网| 老师上课跳d突然被开到最大视频| 国产精品一区二区性色av| 嫩草影院精品99| 精品日产1卡2卡| 在线天堂最新版资源| 床上黄色一级片| 99热这里只有是精品在线观看| 亚洲欧美日韩东京热| 亚洲综合色惰| 麻豆精品久久久久久蜜桃| 国产精品国产高清国产av| 免费av毛片视频| 成人性生交大片免费视频hd| 久久精品国产亚洲av涩爱 | 看片在线看免费视频| 国产精品无大码| 亚洲五月天丁香| 午夜精品一区二区三区免费看| 久久欧美精品欧美久久欧美| 日本av手机在线免费观看| 精品一区二区三区视频在线| 人妻久久中文字幕网| 嫩草影院入口| 国产91av在线免费观看| 女的被弄到高潮叫床怎么办| 伦理电影大哥的女人| 欧美日本亚洲视频在线播放| 免费看美女性在线毛片视频| 国产精品一区www在线观看| 日本三级黄在线观看| 黄色一级大片看看| 婷婷色综合大香蕉| 综合色av麻豆| 热99re8久久精品国产| 久久久精品大字幕| 国产毛片a区久久久久| 欧美日韩在线观看h| 中文字幕熟女人妻在线| 午夜激情欧美在线| 欧美潮喷喷水| 国产精品综合久久久久久久免费| 精品午夜福利在线看| 亚洲精品久久国产高清桃花| 精品99又大又爽又粗少妇毛片| 国产精品乱码一区二三区的特点| 人人妻人人澡欧美一区二区| 成年女人永久免费观看视频| 春色校园在线视频观看| 又爽又黄a免费视频| 国产精品国产高清国产av| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品人妻久久久久久| 看非洲黑人一级黄片| 国产亚洲精品久久久com| 尾随美女入室| 国产国拍精品亚洲av在线观看| 国产一级毛片在线| 男女啪啪激烈高潮av片| 日韩高清综合在线| 欧美激情国产日韩精品一区| 男女边吃奶边做爰视频| 人妻久久中文字幕网| 悠悠久久av| av在线播放精品| 精品久久久久久久久久免费视频| 欧美bdsm另类| 一级黄片播放器| 在线观看66精品国产| 丰满乱子伦码专区| 99久久成人亚洲精品观看| 国产乱人视频| 国内精品宾馆在线| 国产精品一区二区在线观看99 | 亚洲精品久久久久久婷婷小说 | 三级毛片av免费| 白带黄色成豆腐渣| 一边亲一边摸免费视频| 蜜桃久久精品国产亚洲av| 中文字幕av成人在线电影| 天天躁日日操中文字幕| 久久久久久久久久成人| 亚洲国产欧美人成| 久久久久久久亚洲中文字幕| 热99在线观看视频| 婷婷色av中文字幕| 成人特级黄色片久久久久久久| 乱人视频在线观看| 国产精品电影一区二区三区| 国产成人91sexporn| 少妇熟女欧美另类| 在线播放国产精品三级| 黄片无遮挡物在线观看| 国产真实乱freesex| 三级经典国产精品| 国产黄片视频在线免费观看| 免费观看精品视频网站| 我的女老师完整版在线观看| 校园春色视频在线观看| 久久人人精品亚洲av| 毛片女人毛片| 在线观看66精品国产| 欧美性感艳星| 日韩大尺度精品在线看网址| 色噜噜av男人的天堂激情| 51国产日韩欧美| 黄色日韩在线| 国产高清有码在线观看视频| 97超碰精品成人国产| 久久中文看片网| 国产69精品久久久久777片| 三级国产精品欧美在线观看| 插阴视频在线观看视频| 国产精品1区2区在线观看.| or卡值多少钱| 成人午夜精彩视频在线观看| 国产成年人精品一区二区| 亚洲av不卡在线观看| 日韩一本色道免费dvd| 一本久久精品| 久久久久久久久久久免费av| 狠狠狠狠99中文字幕| 日韩一区二区视频免费看| 婷婷六月久久综合丁香| 特大巨黑吊av在线直播| 亚洲欧美精品专区久久| 少妇人妻精品综合一区二区 | 免费搜索国产男女视频| 少妇熟女aⅴ在线视频| 亚洲无线在线观看| 国产精品野战在线观看| 哪里可以看免费的av片| 精品久久久噜噜| 成人国产麻豆网| 一级黄片播放器| 午夜福利在线观看免费完整高清在 | 欧美成人a在线观看| 婷婷色av中文字幕| 亚洲欧美日韩高清专用| 啦啦啦韩国在线观看视频| 只有这里有精品99| 丰满乱子伦码专区| 如何舔出高潮| 国产精品久久久久久久电影| 精品久久久久久久久av| 欧美成人精品欧美一级黄| 人妻夜夜爽99麻豆av| 日本av手机在线免费观看| 国产在线精品亚洲第一网站| 91午夜精品亚洲一区二区三区| 日韩av在线大香蕉| 久久草成人影院| 一区二区三区高清视频在线| av在线播放精品| 啦啦啦观看免费观看视频高清| 精品人妻偷拍中文字幕| 听说在线观看完整版免费高清| 中出人妻视频一区二区| 国产精品一区二区三区四区免费观看| 美女国产视频在线观看| 精品欧美国产一区二区三| 亚洲国产日韩欧美精品在线观看| 哪个播放器可以免费观看大片| 毛片一级片免费看久久久久| 国产综合懂色| 又粗又爽又猛毛片免费看| 免费观看人在逋| 级片在线观看| 1000部很黄的大片| 国产精品精品国产色婷婷| 成熟少妇高潮喷水视频| 亚洲中文字幕日韩| 一个人观看的视频www高清免费观看| 午夜免费男女啪啪视频观看| 91久久精品国产一区二区三区| 亚洲图色成人| 午夜福利高清视频| 3wmmmm亚洲av在线观看| 亚洲精品自拍成人| 丰满乱子伦码专区| 热99re8久久精品国产| 特级一级黄色大片| 亚洲激情五月婷婷啪啪| 青春草亚洲视频在线观看| 校园春色视频在线观看| 中文亚洲av片在线观看爽| 色综合站精品国产| 嘟嘟电影网在线观看| 99视频精品全部免费 在线| 久久久久久国产a免费观看| 干丝袜人妻中文字幕| 在线a可以看的网站| 有码 亚洲区| 国产精品一二三区在线看| h日本视频在线播放| 亚洲精品色激情综合| 欧美色欧美亚洲另类二区| 久久久久国产网址| 变态另类丝袜制服| 一级av片app| 欧美一区二区精品小视频在线| 亚洲最大成人手机在线| 日本-黄色视频高清免费观看| 性色avwww在线观看| 看十八女毛片水多多多| 日本av手机在线免费观看| 欧美日韩国产亚洲二区| 校园人妻丝袜中文字幕| 亚洲国产精品合色在线| 男女视频在线观看网站免费| 三级男女做爰猛烈吃奶摸视频| 男女视频在线观看网站免费| 在线a可以看的网站| 18禁在线无遮挡免费观看视频| 色综合色国产| 天天躁日日操中文字幕| 99视频精品全部免费 在线| 97热精品久久久久久| av天堂中文字幕网| 午夜福利视频1000在线观看| 三级国产精品欧美在线观看| 乱系列少妇在线播放| 成人毛片60女人毛片免费| 亚洲欧美日韩高清专用| 悠悠久久av| 国产精品嫩草影院av在线观看| 全区人妻精品视频| 国产片特级美女逼逼视频| 亚洲成人久久性| 亚洲国产色片| 男人舔奶头视频| 久久人人爽人人片av| 此物有八面人人有两片| 成年免费大片在线观看| 婷婷色综合大香蕉| 亚洲不卡免费看| www.色视频.com| 久久精品国产鲁丝片午夜精品| 男女做爰动态图高潮gif福利片| 国产久久久一区二区三区| 亚洲欧美成人综合另类久久久 | 亚洲国产精品国产精品| 国产私拍福利视频在线观看| 国产成人a区在线观看| 国产一区亚洲一区在线观看| 精品国内亚洲2022精品成人| 亚洲欧洲日产国产| 日日摸夜夜添夜夜添av毛片| 一个人免费在线观看电影| 欧美+日韩+精品| 一个人免费在线观看电影| 国产精品乱码一区二三区的特点| 毛片一级片免费看久久久久| 你懂的网址亚洲精品在线观看 | 亚洲在线观看片| АⅤ资源中文在线天堂| 国产精品免费一区二区三区在线| 国产亚洲av片在线观看秒播厂 | 一区二区三区高清视频在线| 一级黄片播放器| 只有这里有精品99| 亚洲av不卡在线观看| 一区福利在线观看| 亚洲精品乱码久久久v下载方式| 菩萨蛮人人尽说江南好唐韦庄 | 人人妻人人看人人澡| 青春草亚洲视频在线观看| 搞女人的毛片| 国产成人a∨麻豆精品| 日本成人三级电影网站| 精品一区二区三区视频在线| 国产极品天堂在线| 国产一级毛片在线| 搡女人真爽免费视频火全软件| av黄色大香蕉| 美女黄网站色视频| 大又大粗又爽又黄少妇毛片口| 久久久久久国产a免费观看| 在线播放无遮挡| 校园人妻丝袜中文字幕| 在线观看66精品国产| 禁无遮挡网站| 一进一出抽搐gif免费好疼| 狂野欧美白嫩少妇大欣赏| 一区二区三区高清视频在线| 国产亚洲欧美98| 美女大奶头视频| 久久草成人影院| 久久这里只有精品中国| 国内精品久久久久精免费| 成年女人永久免费观看视频| 熟妇人妻久久中文字幕3abv| videossex国产| 亚洲一区二区三区色噜噜| 狂野欧美白嫩少妇大欣赏| 国产精品福利在线免费观看| 免费黄网站久久成人精品| 免费电影在线观看免费观看| 麻豆一二三区av精品| 成人漫画全彩无遮挡| 久久鲁丝午夜福利片| 欧美成人a在线观看| 午夜福利在线观看免费完整高清在 | 国产毛片a区久久久久| 欧美成人一区二区免费高清观看| 日本黄大片高清| 人妻少妇偷人精品九色| 国产真实伦视频高清在线观看| 精品少妇黑人巨大在线播放 | 一级av片app| 亚洲精品自拍成人| 麻豆国产av国片精品| 日本-黄色视频高清免费观看| 18禁在线播放成人免费| 国产日韩欧美在线精品| 国产免费一级a男人的天堂| 99久久中文字幕三级久久日本| 美女高潮的动态| 日本黄大片高清| 久久韩国三级中文字幕| 欧美变态另类bdsm刘玥| 草草在线视频免费看| 热99re8久久精品国产| 成人特级黄色片久久久久久久| 激情 狠狠 欧美| 免费观看精品视频网站| 97人妻精品一区二区三区麻豆| 国产69精品久久久久777片| 国产精品福利在线免费观看| 中文字幕制服av| 九九久久精品国产亚洲av麻豆| 尤物成人国产欧美一区二区三区| 美女黄网站色视频| 在线播放国产精品三级| 在线免费观看的www视频| 国产在线精品亚洲第一网站| 99在线人妻在线中文字幕| 亚洲精华国产精华液的使用体验 | 亚洲无线观看免费| 69人妻影院| 色综合色国产| 最近手机中文字幕大全| av在线蜜桃| 欧美xxxx性猛交bbbb| 国产三级在线视频| 国产v大片淫在线免费观看| 97人妻精品一区二区三区麻豆| 亚洲人成网站在线播| 国产精品一及| 免费一级毛片在线播放高清视频| 久久久久久大精品| 国产伦一二天堂av在线观看| 成人午夜高清在线视频| 成年版毛片免费区| 国产成人影院久久av| 国产女主播在线喷水免费视频网站 | 亚洲人成网站在线播放欧美日韩| 日韩亚洲欧美综合| 日产精品乱码卡一卡2卡三| 欧美3d第一页| 久久国产乱子免费精品| 久久久久性生活片| 插逼视频在线观看| 亚洲中文字幕一区二区三区有码在线看| or卡值多少钱| 国产av一区在线观看免费| 欧美日韩精品成人综合77777| 国产日本99.免费观看| 欧美精品一区二区大全| 在现免费观看毛片| 色哟哟哟哟哟哟| 天美传媒精品一区二区| kizo精华| 99九九线精品视频在线观看视频| 小蜜桃在线观看免费完整版高清| 国产一区亚洲一区在线观看| 天堂网av新在线| 国产精品免费一区二区三区在线| 久久久久久久久久久丰满| 久久精品国产亚洲av香蕉五月| 中文字幕免费在线视频6| 久久国产乱子免费精品| 干丝袜人妻中文字幕| 91av网一区二区| 91精品国产九色| 又粗又爽又猛毛片免费看| 久久人人爽人人片av| 久久久午夜欧美精品| 午夜老司机福利剧场| 亚洲经典国产精华液单| 综合色丁香网| 亚洲精华国产精华液的使用体验 | 日韩欧美精品免费久久| 黄色配什么色好看| а√天堂www在线а√下载| 久久亚洲国产成人精品v| 一级二级三级毛片免费看| 国产精品人妻久久久久久| 中文字幕熟女人妻在线| 免费av观看视频| 日韩强制内射视频| 亚洲丝袜综合中文字幕| 久久久久久久久久黄片| 亚洲在线观看片| 国产精品免费一区二区三区在线| av黄色大香蕉| 日韩中字成人| 嘟嘟电影网在线观看| 国产 一区精品| 午夜精品国产一区二区电影 | 国产精品一区www在线观看| 亚洲一区二区三区色噜噜| 好男人在线观看高清免费视频| 91久久精品电影网| 人妻系列 视频| 久久人人精品亚洲av| 午夜精品国产一区二区电影 | 国内揄拍国产精品人妻在线| 中文字幕制服av| 久久中文看片网| 国产在线男女| 岛国在线免费视频观看| 天天一区二区日本电影三级| 好男人视频免费观看在线| 我的老师免费观看完整版| 99热6这里只有精品| 美女大奶头视频| 亚洲天堂国产精品一区在线| 亚洲av.av天堂| 亚洲三级黄色毛片| 国国产精品蜜臀av免费| 亚洲av中文字字幕乱码综合| 一个人免费在线观看电影| 日本成人三级电影网站| 亚洲国产精品成人久久小说 | 婷婷色av中文字幕| 菩萨蛮人人尽说江南好唐韦庄 | 青春草国产在线视频 | 免费观看人在逋| 伦理电影大哥的女人| 毛片一级片免费看久久久久| 老熟妇乱子伦视频在线观看| 亚洲高清免费不卡视频| 国产精品.久久久| 亚洲七黄色美女视频| 久久久久久久久中文| 最近视频中文字幕2019在线8| 搞女人的毛片| 99riav亚洲国产免费| 日韩人妻高清精品专区| 欧美成人精品欧美一级黄| 亚洲一级一片aⅴ在线观看| 亚洲欧美成人精品一区二区| 夜夜爽天天搞| 精品熟女少妇av免费看| 麻豆精品久久久久久蜜桃| av免费观看日本| 国产高清三级在线| 综合色av麻豆| 亚洲,欧美,日韩| 亚洲高清免费不卡视频| 我的老师免费观看完整版| 国产精品伦人一区二区| 97人妻精品一区二区三区麻豆| 天堂影院成人在线观看| 丰满乱子伦码专区| 亚洲色图av天堂| 国产精品一区www在线观看| 亚洲成人中文字幕在线播放| 国产成人a区在线观看| 国产黄色小视频在线观看| 国产精品1区2区在线观看.| 国产黄色视频一区二区在线观看 | 久久久久久九九精品二区国产| 青春草国产在线视频 | 一区二区三区免费毛片| 两个人的视频大全免费|