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

    運(yùn)用氮、氧雙同位素技術(shù)研究永安江硝酸鹽來源

    2016-11-04 05:37:11吳文歡何小娟吳海露王欣澤
    關(guān)鍵詞:永安洱海牲畜

    吳文歡,何小娟,吳海露,王欣澤,沈 劍

    (上海交通大學(xué)環(huán)境科學(xué)與工程學(xué)院,上海 200240)

    ?

    運(yùn)用氮、氧雙同位素技術(shù)研究永安江硝酸鹽來源

    吳文歡,何小娟①,吳海露,王欣澤,沈 劍

    (上海交通大學(xué)環(huán)境科學(xué)與工程學(xué)院,上海 200240)

    以洱海入湖河流永安江為研究對(duì)象,利用硝酸鹽δ15N和δ18O雙同位素技術(shù)對(duì)永安江水體的硝酸鹽氮來源進(jìn)行識(shí)別。在永安江沿程共布置9個(gè)監(jiān)測(cè)點(diǎn),分析硝酸鹽的污染特征,并利用離子交換樹脂法對(duì)水樣進(jìn)行預(yù)處理后測(cè)試硝酸鹽δ15N和δ18O。結(jié)果表明,永安江硝酸鹽氮源負(fù)荷占永安江總氮源污染的50%左右,各采樣點(diǎn)ρ(硝酸鹽)為0.07~5.22 mg·L-1,均值為1.00~2.39 mg·L-1。經(jīng)同位素測(cè)試,各采樣點(diǎn)δ15N-NO3-均值為6.12‰~13.88‰,δ18O-NO3-均值為8.24‰~11.72‰;永安江河水中硝酸鹽主要來自于流域內(nèi)化學(xué)肥料、牲畜糞便、生活污水和土壤有機(jī)氮硝化;利用IsoSource混合模型對(duì)4種形態(tài)的硝酸鹽來源進(jìn)行定量分析,發(fā)現(xiàn)化學(xué)肥料占37.3%,牲畜糞便占34.6%,村落污水占18.2%,土壤有機(jī)氮占9.9%。利用IsoSource混合模型可為河流硝酸鹽來源定量研究提供新的研究思路,硝酸鹽貢獻(xiàn)比例與河流流經(jīng)村落位置及土地利用類型有關(guān)。

    永安江;硝酸鹽污染;氮氧同位素;來源示蹤;IsoSource模型;污染貢獻(xiàn)率

    洱海是云南省第2大湖泊,是大理飲用水水源地和重要的旅游資源,隨著洱海流域近些年來大蒜等經(jīng)濟(jì)作物種植量的不斷擴(kuò)大和牲畜養(yǎng)殖業(yè)的興起,大量化肥和牲畜糞便排放對(duì)洱海造成的污染越來越大,尤其是硝酸鹽氮的污染,常年監(jiān)測(cè)數(shù)據(jù)顯示部分入湖河流中ρ(硝酸鹽)已達(dá)11.61 mg·L-1,超出GB 3838—2002《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》中集中式生活飲用水地表水源地補(bǔ)充項(xiàng)目標(biāo)準(zhǔn)限值(10 mg·L-1)。識(shí)別入湖河流中硝酸鹽氮的主要來源是控制氮污染的前提。通常采用統(tǒng)計(jì)學(xué)手段估算氮從陸地遷移至水體的量,即較多地采用水質(zhì)理化指標(biāo)、土地利用類型、人口數(shù)量、化肥使用量和牲畜數(shù)量等來估算氮排放量的方法[1-2]。但這種方法可能存在資料不完善、統(tǒng)計(jì)不準(zhǔn)確的問題,且流域內(nèi)硝酸鹽氮有多種來源,不同污染源會(huì)受到遷移轉(zhuǎn)化過程中多種因素(地形、轉(zhuǎn)移途徑、處理工藝)的影響,使得氮在流域及水體中的遷移轉(zhuǎn)化過程非常復(fù)雜,不同區(qū)域的非點(diǎn)源污染隨著人口密度和土地利用方式不同也存在較大差異[3]。因此,筆者采用硝酸鹽氮氧雙穩(wěn)定同位素溯源技術(shù)來識(shí)別硝酸鹽的污染狀況,可以彌補(bǔ)統(tǒng)計(jì)學(xué)手段調(diào)查的局限性,為硝酸鹽的污染控制提供重要技術(shù)手段。該方法在國(guó)內(nèi)外多被應(yīng)用于地下水、海洋等水體[4-5],應(yīng)用于地表水研究還較少[6-11];且大多數(shù)研究是基于硝酸鹽氮氧穩(wěn)定同位素進(jìn)行的定性分析,定量研究較少[12]。以洱海流域永安江為研究對(duì)象,利用硝酸鹽δ15N和δ18O值對(duì)永安江水體中硝酸鹽的主要潛在來源進(jìn)行識(shí)別,并結(jié)合Iso-Source混合模型對(duì)污染源進(jìn)行定量分析。該研究將流域內(nèi)的村落污水、牲畜糞便、化學(xué)肥料和土壤有機(jī)氮等作為1個(gè)系統(tǒng)來分析,可為永安江流域氮污染的控制、治理提供科學(xué)依據(jù)。自然界中氮元素的穩(wěn)定同位素有15N和14N 2種,其氮原子中15N和14N的相對(duì)豐度分別為0.366 3%和99.633 7%。而氧原子的穩(wěn)定同位素有16O、17O和18O 3種,其相對(duì)豐度分別為99.759%、0.037%和0.204%[13]。

    樣品中的的同位素比值通常用它與某一標(biāo)準(zhǔn)物質(zhì)的同位素比值的千分差δ值來表示,其定義式如下:

    δ=[(Ra-Rs)/Rs]×1 000。

    (1)

    式(1)中,δ為同位素比率;Rs為樣品中15N/14N或18O/17O的比值;Ra為大氣氮和維也納海水中的15N/14N或18O/17O的比值。

    氮在自然界循環(huán)過程中會(huì)發(fā)生礦化、同化、硝化和反硝化作用等一系列生化反應(yīng),均會(huì)導(dǎo)致氮同位素發(fā)生分餾作用,因此不同來源的硝酸鹽δ15N值會(huì)有所差異[14]。同樣,硝酸鹽中氧的同位素值δ18O也受到不同來源的影響,具有不同的同位素值,這是氮源來源識(shí)別的理論基礎(chǔ)。此外,不同硝酸鹽氮源具有不同的穩(wěn)定同位素值,且具有一定的規(guī)律性,進(jìn)入水體后保持相對(duì)的穩(wěn)定性,因而可以利用δ15N和δ18O同位素的特征值,溯源水體硝態(tài)氮。目前多數(shù)學(xué)者利用不同氮源中硝酸鹽的氮氧穩(wěn)定同位素δ15N和δ18O值的特征值來識(shí)別水體中的硝酸鹽來源,進(jìn)行定性分析。大量的研究和文獻(xiàn)表明,不同氮源來源的硝酸鹽δ15N和δ18O值范圍如下:由無機(jī)化學(xué)肥料產(chǎn)生的δ15N值為-6‰~6‰,δ18O值為17‰~25‰[15];由牲畜糞便產(chǎn)生的δ15N值為5‰~25‰,δ18O值為-5‰~7‰[16];由生活污水產(chǎn)生的δ15N值為4‰~19‰,δ18O值為-5‰~10‰[17];土壤中有機(jī)氮礦化產(chǎn)生的δ15N值為-3‰~8‰,δ18O值為-10‰~10‰[15]。

    1 研究區(qū)概況

    永安江流域位于洱海北部,是洱海流域的一部分,地處北緯25°25′~26°16′,東經(jīng)99°32′~100°27′,位于瀾滄江、金沙江和元江3大水系分水嶺地帶,屬瀾滄江—湄公河水系,河流全長(zhǎng)18.35 km,流域總面積約110.25 km2,多年平均徑流量為0.38億m3,占洱海多年平均徑流量的5.2%。流域內(nèi)村落密集,人口眾多,農(nóng)田數(shù)量大,流域總?cè)丝?4 975人,占洱海流域總?cè)丝诘?.2%。流域內(nèi)工業(yè)較少,牲畜養(yǎng)殖和種植業(yè)較為發(fā)達(dá)。永安江是洱海流域典型的河流,流經(jīng)農(nóng)田、村落等多種土地利用類型,因而掌握永安江流域的污染狀況并加以控制可以削減洱海氮污染負(fù)荷。

    2 樣品采集與測(cè)試方法

    2.1 樣品采集

    如圖1所示,根據(jù)研究區(qū)域河流流經(jīng)區(qū)域土地利用類型、村落等情況,在永安江共布置9個(gè)點(diǎn),分別為下山口(Y1)、三枚村(Y2)、大花園(Y3)、大樹營(yíng)濕地入口(Y4)、簸箕村(Y5)、文筆湖(Y6)、青索(Y7)、江尾(Y8)和洱海入湖口(Y9)。用采樣器在永安江河流中間取水面下0.5 m深水樣,用高密度聚乙烯塑料采樣瓶取樣密封,測(cè)定總氮(TN)、硝態(tài)氮(NO3-)、氨氮(以NH4+計(jì))、氯離子(Cl-)濃度等化學(xué)指標(biāo),其中總氮濃度采用堿性過硫酸鉀消解紫外線分光光度法(HJ 636—2012)、硝態(tài)氮濃度采用紫外分光光度法(HJ/T 346—2007)、氨氮濃度采用水楊酸分光光度法(HJ 536—2009)、氯離子濃度采用硝酸銀滴定法(GB 11896—89)測(cè)定?,F(xiàn)場(chǎng)使用便攜式溶氧儀(HACH-HQ30d)進(jìn)行溶解氧(DO)、溫度等測(cè)試。另外,使用10 L高密度聚乙烯塑料桶取水樣,密封、冷藏運(yùn)回實(shí)驗(yàn)室進(jìn)行同位素樣品預(yù)處理。采樣時(shí)間為2015年2—5月,采樣頻率為每月2次,共采樣8次。此外,對(duì)流域內(nèi)所使用的化肥、牲畜糞便、土壤、村落生活污水采樣,其中所采化肥樣品為當(dāng)?shù)氐屎蛷?fù)合肥,按當(dāng)?shù)厥褂昧颗浔然旌稀瞱(氮肥)=60%,w(復(fù)合肥)=30%)〕,用去離子水浸泡24 h,預(yù)處理方式與水樣一致,共采樣6次;牲畜糞便和土壤按照每月2次的頻率共采樣8次,同樣使用去離子水浸泡24 h,放置搖床6 h,過濾后與水樣處理方式一致。預(yù)處理后的樣品供硝酸鹽δ15N和δ18O同位素值測(cè)試用。

    Y1—下山口;Y2—三枚村;Y3—大花園;Y4—大樹營(yíng)濕地入口;Y5—簸箕村;Y6—文筆湖;Y7—青索;Y8—江尾;Y9—洱海入湖口。

    2.2 樣品的測(cè)試方法

    同位素測(cè)試的預(yù)處理采用離子交換樹脂法[18-20]。

    δ15N測(cè)試預(yù)處理:帶回實(shí)驗(yàn)室的水樣用0.45 μm孔徑的聚碳酸酯膜過濾,取國(guó)產(chǎn)強(qiáng)堿201×7(701型)苯乙烯氯型陰離子交換樹脂6 g置于樹脂柱中,過濾后的水樣通過蠕動(dòng)泵控制流速8.33~16.67 mL·min-1,通過陰離子交換樹脂柱對(duì)硝酸根離子進(jìn)行富集,并用15 mL 3 mol·L-1的鹽酸分5次洗脫離子交換樹脂,在洗脫液中分多次加入總量為6.5 g的過量Ag2O(在水浴中進(jìn)行),充分?jǐn)嚢?中和混合液pH值為5.5~6.0,過濾去除中和液中沉淀的AgCl,過濾后的約為40 mL的硝酸鹽溶液收集于塑料杯中,將其分成2份,一份放在冷凍干燥機(jī)(SCIENTZ-12N)中凍干成AgNO3粉末,避光保存,用來測(cè)試δ15N;另一份進(jìn)行同位素氧測(cè)試的預(yù)處理。

    δ18O測(cè)試預(yù)處理:為了去除其他含氧離子(PO43-、CO32-和SO42-等)對(duì)于測(cè)試的影響,在AgNO3中和液中加入3 mol·L-1BaCl2溶液,置冷藏過夜后,將其形成的沉淀過濾去除,然后再通過陽離子交換樹脂(001×7,H+型樹脂)去除過多的Ba2+,再加入Ag2O中和后去除沉淀,確保pH值為6左右,然后在濾后液中加入一定量的活性炭并不斷攪動(dòng)除去溶解有機(jī)質(zhì),最后在凍干機(jī)中凍干成AgNO3粉末。

    預(yù)處理后的樣品在國(guó)家海洋局第三海洋研究所進(jìn)行δ15N和δ18O值測(cè)試。測(cè)試儀器是Thermo Fisher公司的元素分析儀-穩(wěn)定同位素質(zhì)譜儀聯(lián)機(jī)(Flash EA 1112 HT-Delta V Advantages)。δ15N值的測(cè)試是利用氦氣(He)作為載氣,流速為90 mL·min-1,在反應(yīng)管溫度為960 ℃、色譜柱溫度為50 ℃條件下,將樣品中的氮轉(zhuǎn)化為氮?dú)鉁y(cè)試δ15N;δ18O的測(cè)試同樣使用氦氣作為載氣,流速為90 mL·min-1,在高溫裂解管溫度1 400 ℃、色譜柱溫度85 ℃ 條件下,將樣品中的氧轉(zhuǎn)化為CO2測(cè)試δ18O[21]。δ15N和δ18O值分別以大氣氮(AIR)和大洋海水(SMOW)國(guó)際標(biāo)準(zhǔn)作為參考標(biāo)準(zhǔn),δ15N和δ18O 值的分析精度均為±0.3‰,能滿足該研究的精度要求。

    2.3 IsoSource模型

    IsoSource模型是以質(zhì)量平衡混合模型為基礎(chǔ),開發(fā)用于計(jì)算源頭值超過n+1(n為多種同位素值)的潛在貢獻(xiàn)比例,最初主要應(yīng)用于食物網(wǎng)、植物水分供給的研究[22-23]。在該模型中,在一定增量范圍內(nèi),使用標(biāo)準(zhǔn)線型混合模型來模擬每一種可能污染比例(和為1),模型在混合計(jì)算中滿足同位素質(zhì)量守恒,通過測(cè)試的同位素值信息來確定混合物各部分的比例范圍,因而同樣可以用來估算各種硝酸鹽來源對(duì)河流硝酸鹽污染的貢獻(xiàn)率,計(jì)算出的每個(gè)解代表了一個(gè)資源百分比的組合。在軟件中設(shè)置好模型的增量參數(shù)(increment)和容差參數(shù)(tolerance),模型利用迭代方法計(jì)算出水樣中不同污染來源所占貢獻(xiàn)率的概率分布圖,并給出所有來源計(jì)算結(jié)果的平均值。不同來源所有可能的百分比組合則按下式計(jì)算[12]:

    式(2)中,i為增量參數(shù);s為污染源數(shù)量。該模型可以將各種潛在污染源貢獻(xiàn)的可能組合都顯示出來。將各個(gè)組合中模擬出來的各比例結(jié)果相加,如果相加之和的值在設(shè)定的容差范圍之內(nèi)(如±0.01),即可認(rèn)為該模擬的組合值是合適的解[24]。

    3 結(jié)果與分析

    3.1 永安江流域硝酸鹽氮污染狀況

    從硝酸鹽監(jiān)測(cè)結(jié)果看,永安江流域硝酸鹽氮約占總氮的50%,是主要的氮源成分,在監(jiān)測(cè)期間ρ(硝酸鹽)為0.01~5.01mg·L-1(圖2),均值為1.07~2.97mg·L-1,總體來說硝酸鹽濃度從上游到下游呈下降趨勢(shì),這可能與近年來當(dāng)?shù)卣畬?duì)永安江河道進(jìn)行了一系列的整治有關(guān),尤其是在河水流經(jīng)流域中部的大樹營(yíng)濕地后,河水中硝酸鹽濃度減少近一半。Y5樣點(diǎn)硝酸鹽濃度較高,可能原因是該采樣點(diǎn)與村落居民生活區(qū)(距離2m)和農(nóng)田尾水排出口非常接近。

    從硝酸鹽氮氧同位素值的測(cè)試結(jié)果來看(表1),硝酸鹽δ15N值為8.03‰~13.88‰,各采樣點(diǎn)均值分別為8.24‰~11.72‰,其中Y3的δ15N均值最低,Y2的δ15N均值最高,這表明Y2牲畜糞便養(yǎng)殖污染較大。δ18O值為8.17‰~17.65‰,各采樣點(diǎn)均值分別為10.60‰~13.76‰,其中Y3的δ18O均值最低,Y7最高,可能是因?yàn)閅7人口較密集,受生活污水的污染增大,而且能看到Y(jié)7采樣點(diǎn)的ρ(Cl-)(均值27.92mg·L-1)較其他采樣點(diǎn)較高,這也從側(cè)面反映出該點(diǎn)受村落污水污染較其他采樣點(diǎn)嚴(yán)重。

    Y1—下山口;Y2—三枚村;Y3—大花園;Y4—大樹營(yíng)濕地入口;Y5—簸箕村;Y6—文筆湖;Y7—青索;Y8—江尾;Y9—洱海入湖口。

    表1 永安江采樣點(diǎn)水化學(xué)指標(biāo)及硝酸鹽氮、氧同位素

    Table 1 Chemical indices,15N-NO3-and18O-NO3-in the water samples collected from the Yongan River at the sampling sites

    采樣點(diǎn)ρ(DO)/(mg·L-1)ρ(NH4+)/(mg·L-1)ρ(TN)/(mg·L-1)ρ(NO3-)/(mg·L-1)ρ(Cl-)/(mg·L-1)δ15N-NO3-/‰δ18O-NO3-/‰Y17.55±0.510.22±0.205.57±3.972.22±1.0712.00±4.428.65±0.7412.07±0.62Y25.77±1.530.43±0.264.42±3.032.39±1.2810.85±3.129.72±1.1113.42±1.02Y38.04±0.990.26±0.134.88±2.492.09±0.8620.42±9.538.24±1.8510.60±2.51Y47.35±0.650.33±0.214.32±1.772.01±0.9321.66±11.225.84±1.5813.88±3.41Y59.54±4.070.26±0.244.69±1.812.91±1.2036.11±20.889.25±1.1116.00±0.34Y69.52±1.930.22±0.222.86±1.321.07±0.6720.91±17.048.67±3.1611.75±3.11Y75.19±1.500.17±0.132.67±1.411.00±0.5227.92±15.477.77±1.7513.76±1.47Y86.40±1.540.15±0.143.10±1.781.11±0.7319.66±16.969.23±2.9615.77±0.57Y96.35±1.730.09±0.062.65±1.861.34±0.6522.91±18.288.13±1.7115.52±0.65

    數(shù)據(jù)為平均值±標(biāo)準(zhǔn)差。Y1為下山口;Y2為三枚村;Y3為大花園;Y4為大樹營(yíng)濕地入口;Y5為簸箕村;Y6為文筆湖;Y7為青索;Y8為江尾;Y9為洱海入湖口。

    3.2 定性分析

    利用硝酸鹽氮氧雙同位素識(shí)別河流中的硝酸鹽來源,是基于硝酸鹽氮進(jìn)入水體保持一定的穩(wěn)定性,但實(shí)際研究中氮是一種易分餾的元素,且反硝化是影響氮穩(wěn)定同位素的主要影響因素。反硝化改變了硝酸鹽來源同位素的組成,會(huì)影響最終的溯源結(jié)果,因此識(shí)別河流中硝酸鹽是否發(fā)生了反硝化作用是識(shí)別硝酸鹽來源的前提。有研究表明,反硝化將導(dǎo)致δ15N和δ18O均在硝酸鹽中富集,如果δ15N與δ18O之間存在一種結(jié)余為1.3∶1~ 2.1∶1的線性富集比例關(guān)系,表明可能存在著一定的反硝化作用[25]。此外,當(dāng)水體中ρ(DO)超過3.1 mg·L-1,則不利于反硝化作用的發(fā)生[26]。整個(gè)采樣期間各采樣點(diǎn)ρ(DO)為3.42~13.33 mg·L-1(平均值為7.49 mg·L-1),可以確定污染源進(jìn)入水體不發(fā)生反硝化作用,從表1可以看出δ15N與δ18O的比值不存在1.3∶1~2.1∶1區(qū)間內(nèi)的線性比例關(guān)系,因而可以確定河流中幾乎不存在反硝化作用。

    將采樣測(cè)試的同位素?cái)?shù)據(jù)布置在同位素值典型值分布圖(圖3)中,可以看出,各采樣點(diǎn)δ15N-NO3-和δ18O-NO3-的值主要集中在化肥或者牲畜糞便及污水區(qū)域或者其附近,結(jié)合污染源調(diào)查情況,確定永安江河流的硝態(tài)氮污染主要來自于化學(xué)肥料、牲畜糞便、村落污水和土壤有機(jī)氮。對(duì)上述這4種來源進(jìn)行同位素值測(cè)試,結(jié)果如表2所示。

    3.3 定量分析

    經(jīng)過IsoSource模型計(jì)算,該研究所得的2—5月各采樣點(diǎn)硝酸鹽各污染源比例如表3所示。通過IsoSource混合模型計(jì)算的結(jié)果顯示,3—4月變化不大,主要原因可能是洱海流域在該時(shí)段內(nèi)為旱季,但5月各污染來源變化較為明顯,其中化肥和牲畜糞便對(duì)于硝酸鹽貢獻(xiàn)比例明顯上升。通過實(shí)地采樣也發(fā)現(xiàn),5月初為洱海流域的水稻種植季節(jié),農(nóng)民會(huì)對(duì)農(nóng)田利用牲畜糞便和化肥進(jìn)行泡田,因而該時(shí)間段內(nèi)化肥和牲畜糞便的占比有較為明顯的上升,這與試驗(yàn)數(shù)據(jù)相吻合。

    圖3 不同硝酸鹽來源的δ15N和δ18O典型值與樣品值分布

    表2 永安江流域潛在硝酸鹽氮氧同位素組成比例

    Table 2 Isotopic composition of potential nitrate sources in the Yongan River basin

    來源δ15N-NO3-δ18N-NO3-化學(xué)肥料2.01±0.1220.10±3.21生活污水8.01±1.218.74±1.98牲畜糞便19.98±2.109.92±1.82土壤有機(jī)氮2.87±0.212.03±0.21

    表3 永安江河流硝酸鹽污染來源比例定量計(jì)算結(jié)果

    Table 3 Ratio of the contributions of nitrate pollutant in the Yongan River from various sources

    %

    “—”表示該同位素值沒測(cè)出。Y1為下山口;Y2為三枚村;Y3為大花園;Y4為大樹營(yíng)濕地入口;Y5為簸箕村;Y6為文筆湖;Y7為青索;Y8為江尾;Y9為洱海入湖口。

    4 結(jié)論

    1)永安江河水中硝酸鹽氮污染占總氮污染的50%,在永安江上游(大樹營(yíng)濕地以上)ρ(硝酸鹽)超過2 mg·L-1,經(jīng)過大樹營(yíng)濕地處理凈化后,硝酸鹽濃度下降50%左右。

    (2)根據(jù)IsoSource混合模型計(jì)算結(jié)果,就所有采樣點(diǎn)4個(gè)月份的均值來說,永安江硝酸鹽來自于化學(xué)肥料、牲畜糞便、村落污水和土壤有機(jī)氮的比例分別為37.3%、34.6%、18.2%和9.9%。從各個(gè)采樣點(diǎn)來看,硝酸鹽來源貢獻(xiàn)比有所差異,其中化肥占比為27.1%~56.0%,其次是牲畜糞便,占13.1%~56.0%,村落污水和土壤有機(jī)氮分別占5.0%~35.2%和2.0%~16.4%,河流中硝酸鹽來源與土地利用類型及村落位置有關(guān)。

    (3)永安江硝酸鹽污染源及貢獻(xiàn)比例各季節(jié)有所不同。5月化肥和牲畜糞便污染貢獻(xiàn)比例較其他3個(gè)月有所增加。

    [1] 夏曉樹,張松.貴州省點(diǎn)污染源調(diào)查分析及治理對(duì)策[J].貴州水力發(fā)電,2006(2):8-10.

    [2] 馮素珍,李暢游.內(nèi)蒙古烏梁素海污染源調(diào)查研究[J].中國(guó)農(nóng)村水利水電,2010(4):125-128.

    [3] 秦伯強(qiáng),楊柳燕,陳非洲,等.湖泊富營(yíng)養(yǎng)化發(fā)生機(jī)制與控制技術(shù)及其應(yīng)用[J].科學(xué)通報(bào),2006,51(16):1857-1866.

    [4] 吳登定,姜月華,賈軍遠(yuǎn),等.運(yùn)用氮、氧同位素技術(shù)判別常州地區(qū)地下水氮污染源[J].水文地質(zhì)工程地質(zhì),2006,33(3):11-15.

    [5] 張翠云,張勝,李政紅,等.利用氮同位素技術(shù)識(shí)別石家莊市地下水硝酸鹽污染源[J].地球科學(xué)進(jìn)展,2004,19(2):183-191.

    [6] LEE K S,BONG Y S,LEE D,etal.Tracing the Sources of Nitrate in the Han River Watershed in Korea,Usingδ15N-NO3-andδ18O-NO3-Values[J].Science of the Total Environment,2008,395:117-124.

    [7] PARDO L H,KENDALL C,PETT-RIDGE J,etal.Evaluating the Source of Streamwater Nitrate Usingδ15N andδ18O in Nitrate in Two Watersheds in New Hampshire,USA[J].Hydrological Processes,2004,18(14):2699-2712.

    [8] DEUTSCH B,MEWES M,LISKOWI I,etal.Quantification of Diffuse Nitrate Inputs Into a Small River System Using Stable Isotopes of Oxygen and Nitrogen in Nitrate[J].Organic Geochemistry,2006,37(10):1333-1342.

    [9] 張東,楊偉,趙建立.氮同位素控制下黃河及其主要支流硝酸鹽來源分析[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2012,28(6):622-627.

    [10]李清光,王仕祿.滇池流域硝酸鹽污染的氮氧同位素示蹤[J].地球與環(huán)境,2012,40(3):321-327[11]任玉芬,張心昱,王效科,等.北京城市地表河流硝酸鹽氮來源的氮氧同位素示蹤研究[J].環(huán)境工程學(xué)報(bào),2013,7(5):1636-1640.

    [12]盧麗,李文莉,裴建國(guó),等.基于IsoSource的桂林寨底地下河硝酸鹽來源定量研究[J].地球?qū)W報(bào),2014,35(2):248-254.

    [13]COOK G A,LAUER C M.Oxygen[M]∥CLIDFFORD,HAMPLE A.The Encyclopedia of the Chemical Elements.New York,USA:Wiley and Sons Publishers,1968:499-512.

    [14]PERTERSON B J,FRY B.Stable Isotopes in Ecosystem Studies[J].Annual Review of Ecology & Systernatics,2003,18(1):293-320.

    [15]DING J T,XI B D,GAO R T,etal.Idengtifing Diffused Nitrate Sources in a Stream in a Argricultural Field Using a Dural Isotopic Approach[J].Science of the Total Environment,2014,484:10-18.

    [16]FINALY J C,STERNER R W,KUMAR S.Isotopic Evidence for In-Lake Production of Accumulating Nitrate in Lake Superior[J].Ecological Application,2007,17(8):2323-2332.

    [17]FREYER H D.Seasonal Variation of 15N/14N Ratios in Atmospheric Nitrate Species[J].Tellus Series B:Chemical & Physical Meteorology,2002,43(1):30-44.

    [18]肖化云,劉叢強(qiáng).水樣硝酸鹽氮同位素分析預(yù)處理方法探討[J].巖礦測(cè)試,2002,21(2):105-108.

    [19]周愛國(guó),陳銀琢.水環(huán)境硝酸鹽氮污染研究新方法:δ15N和δ18O相關(guān)法[J].地球科學(xué),2003,28(2):219-224.

    [20]李利文,甘義群,周愛國(guó),等.離子交換樹脂法高效提取水中硝酸鹽研究[J].環(huán)境科學(xué)與技術(shù),2012(增刊2):54-57.

    [21]毛巍,梁志偉,李偉,等.利用氮、氧穩(wěn)定同位素識(shí)別水體硝酸鹽污染源研究進(jìn)展[J].應(yīng)用生態(tài)學(xué)報(bào),2013,24(4):1146-1152.

    [22]BENSTEAD J P,PRINGLE C M.Testing Isosource:Stable Isotope Analysis of a Tropical Fishery With Diverse Organic Matter Sources[J].Ecology,2008,87(2):326-330.

    [23]姜啟吳,歐志吉,左平.鹽沼植被對(duì)江蘇鹽城濕地生態(tài)系統(tǒng)有機(jī)質(zhì)貢獻(xiàn)的初步研究[J].海洋通報(bào),2012,31(5):547-551.

    [24]HOPKINS J B,FERGUSON J M.Estimating the Diets of Animals Using Stable Isotopes and a Comprehensive Bayesian Mixing Model[J].Plos One,2012,7(1):e28478.

    [25]盧義.基于同位素分析法的氮來源及遷移過程研究:以太滆運(yùn)河—竺山灣流域?yàn)槔齕D].重慶:重慶交通大學(xué),2013.

    [26]陳法錦,李學(xué)輝,賈國(guó)東.氮氧同位素在河流硝酸鹽研究中的應(yīng)用[J].地球科學(xué)進(jìn)展,2007,22(12):1251-1257.

    (責(zé)任編輯: 陳 昕)

    Identification of Sources of Nitrate in the Yongan River With Isotopic Technology.

    WU Wen-huan, HE Xiao-juan, WU Hai-lu, WANG Xin-ze, SHEN Jian

    (School of Environmental Science and Engineering, Shanghai Jiao Tong University, Shanghai 200240, China)

    A case study of the Yongan River, a major inflow of the Erhai Lake, was conducted to identify sources of nitrate-N in the river water with the15N-NO3-and18O-NO3-double isotope technique. Along the river, 9 monitoring sites were set up, responsible for characterizing nitrate pollution of the water, and measuring nitrate15N-NO3-and18O-NO3-in the water samples after treatment with ion exchange resin. Results show that the nitrate sources along the Yongan River are responsible for about 50% of the total nitrogen pollution load in the water. Nitrate concentration in the water samples at the monitoring sites varied in the range of 0.07-5.22 mg·L-1, with mean being in the range of 1.00-2.39 mg·L-1. The isotope test reveals that the mean15N-NO3-of the water samples at the monitoring sites ranged from 6.12‰ to 13.88‰, and the mean18O-NO3-did from 8.24‰ to 11.72‰. The findings demonstrate that the nitrate pollutant in the river comes mainly from chemical fertilization, livestock manure, domestic sewage, and nitrification of organic nitrogen in the soil in the valley. Quantification of contributions of the four sources with the IsoSource mixing model tells that 37.3% was attributed to chemical fertilization, 34.6% to disposal of livestock manure, 18.2% to discharge of rural domestic sewage and 9.9% to nitrification of organic nitrogen in the soil. The experiment shows that the use of the IsoSource mixed model to quantify sources of nitrate pollutant in the river may provide future studies with a new train of thought and that ratio of the contributions of the sources to the total nitrate pollutant in the river is related to location of the villages the river runs through and type of land use therein.

    Yongan River;nitrate pollution;nitrogen and oxygen isotope;source tracing;IsoSource model;pollution contribution rate

    2015-08-25

    國(guó)家水體污染控制與治理科技重大專項(xiàng)(2012ZX07105-003);上海交通大學(xué)青年教師啟動(dòng)計(jì)劃(13X100040042)

    X522

    A

    1673-4831(2016)05-0802-06

    10.11934/j.issn.1673-4831.2016.05.017

    吳文歡(1988—),男,安徽樅陽人,碩士生,主要從事湖泊河流生態(tài)修復(fù)研究。E-mail: wenhuanwu@foxmail.com

    ① 通信作者E-mail: hexiaojuan@sjtu.edu.cn

    猜你喜歡
    永安洱海牲畜
    內(nèi)鏡止血在急性非靜脈曲張性上消化道出血治療中的應(yīng)用
    長(zhǎng)安永安
    交換牲畜
    洱海月下
    黃河之聲(2020年11期)2020-07-16 13:27:12
    洱海,好美
    洱海
    永安橋 百年永安
    愛上洱海,只需要這十個(gè)瞬間
    海峽姐妹(2017年5期)2017-06-05 08:53:17
    竹鄉(xiāng)永安
    農(nóng)藥毒死牲畜 為何諸多說法
    亚洲精品成人av观看孕妇| 天天躁夜夜躁狠狠躁躁| 亚洲图色成人| 91成人精品电影| 久久久a久久爽久久v久久| 久热这里只有精品99| 考比视频在线观看| 十八禁网站网址无遮挡| 成人亚洲欧美一区二区av| 纵有疾风起免费观看全集完整版| 亚洲一区二区三区欧美精品| 99香蕉大伊视频| 日韩中文字幕视频在线看片| √禁漫天堂资源中文www| 欧美性感艳星| 丝袜美足系列| 下体分泌物呈黄色| 天堂俺去俺来也www色官网| 亚洲成人一二三区av| 亚洲伊人色综图| 美女中出高潮动态图| 欧美激情国产日韩精品一区| 香蕉丝袜av| av网站免费在线观看视频| 成人手机av| 桃花免费在线播放| 国产男女内射视频| xxxhd国产人妻xxx| 午夜91福利影院| 久久久久精品性色| 性色av一级| 久久久精品区二区三区| tube8黄色片| 日韩视频在线欧美| 成年动漫av网址| 一级爰片在线观看| 一区二区日韩欧美中文字幕 | 岛国毛片在线播放| 在线精品无人区一区二区三| 日韩成人伦理影院| 黄色一级大片看看| 免费大片黄手机在线观看| 黑丝袜美女国产一区| 你懂的网址亚洲精品在线观看| 全区人妻精品视频| 免费女性裸体啪啪无遮挡网站| 亚洲精品aⅴ在线观看| 草草在线视频免费看| 一本—道久久a久久精品蜜桃钙片| 国国产精品蜜臀av免费| 伦理电影免费视频| 精品少妇久久久久久888优播| 黄色毛片三级朝国网站| 美国免费a级毛片| 亚洲精品成人av观看孕妇| 香蕉精品网在线| 成人综合一区亚洲| 亚洲国产精品999| 一区二区日韩欧美中文字幕 | 亚洲成av片中文字幕在线观看 | 久久亚洲国产成人精品v| 久久久久久久久久成人| 国产在线视频一区二区| 久久久久久久久久久免费av| 精品国产露脸久久av麻豆| 天天躁夜夜躁狠狠久久av| 欧美日韩成人在线一区二区| 国产一区二区在线观看日韩| 最新中文字幕久久久久| 成人国产麻豆网| 国产精品秋霞免费鲁丝片| 日本vs欧美在线观看视频| 日韩精品有码人妻一区| 爱豆传媒免费全集在线观看| 成人漫画全彩无遮挡| 精品国产国语对白av| av电影中文网址| 亚洲精品第二区| 女的被弄到高潮叫床怎么办| 精品亚洲成a人片在线观看| 秋霞伦理黄片| 欧美激情 高清一区二区三区| 免费人成在线观看视频色| 国产伦理片在线播放av一区| kizo精华| 在线精品无人区一区二区三| 看十八女毛片水多多多| 国产精品久久久久久久电影| 久久久久久久久久久久大奶| 波多野结衣一区麻豆| 18禁裸乳无遮挡动漫免费视频| 国产亚洲午夜精品一区二区久久| 制服人妻中文乱码| 高清不卡的av网站| 久久99热6这里只有精品| 亚洲av在线观看美女高潮| av视频免费观看在线观看| 亚洲精品美女久久久久99蜜臀 | 亚洲美女黄色视频免费看| 日本vs欧美在线观看视频| 国产精品秋霞免费鲁丝片| 亚洲国产欧美在线一区| 免费久久久久久久精品成人欧美视频 | 午夜av观看不卡| 国产成人免费观看mmmm| 亚洲婷婷狠狠爱综合网| 国产老妇伦熟女老妇高清| 下体分泌物呈黄色| 99热这里只有是精品在线观看| 国产精品麻豆人妻色哟哟久久| 满18在线观看网站| 久久热在线av| 国产在线一区二区三区精| 人体艺术视频欧美日本| 中文精品一卡2卡3卡4更新| 丝袜脚勾引网站| 国产精品人妻久久久久久| 国产黄频视频在线观看| 亚洲,一卡二卡三卡| 亚洲图色成人| 91午夜精品亚洲一区二区三区| 久久久久人妻精品一区果冻| 亚洲国产精品国产精品| 精品福利永久在线观看| 免费观看性生交大片5| 中文字幕制服av| 色网站视频免费| 精品久久蜜臀av无| 在线观看www视频免费| www.av在线官网国产| 精品久久久久久电影网| 日韩 亚洲 欧美在线| 免费看不卡的av| 91国产中文字幕| 久久 成人 亚洲| 精品亚洲成国产av| 欧美精品一区二区大全| 一区二区三区四区激情视频| 午夜影院在线不卡| 哪个播放器可以免费观看大片| 黄色毛片三级朝国网站| 天天影视国产精品| 这个男人来自地球电影免费观看 | 国产男女内射视频| 在线亚洲精品国产二区图片欧美| 日本黄大片高清| 免费av中文字幕在线| 免费看av在线观看网站| 成年美女黄网站色视频大全免费| 在线观看美女被高潮喷水网站| 18禁观看日本| 美女大奶头黄色视频| 男女国产视频网站| a级毛片黄视频| 婷婷成人精品国产| 美女中出高潮动态图| 国产一区二区三区综合在线观看 | 久久国产亚洲av麻豆专区| 国产精品熟女久久久久浪| 9热在线视频观看99| 免费黄网站久久成人精品| 亚洲欧美成人精品一区二区| 精品福利永久在线观看| 在线观看免费高清a一片| 97在线人人人人妻| 丝袜人妻中文字幕| 中文字幕精品免费在线观看视频 | 亚洲,一卡二卡三卡| 丝袜喷水一区| 九九在线视频观看精品| 看十八女毛片水多多多| 国产精品久久久av美女十八| 熟女电影av网| 超色免费av| 久久久精品94久久精品| 日韩在线高清观看一区二区三区| 成人毛片60女人毛片免费| 男女高潮啪啪啪动态图| 丰满少妇做爰视频| 啦啦啦啦在线视频资源| 蜜桃在线观看..| 成人亚洲精品一区在线观看| freevideosex欧美| 啦啦啦在线观看免费高清www| 亚洲国产最新在线播放| 自线自在国产av| 看免费成人av毛片| 国产精品一国产av| 日本欧美国产在线视频| 国产精品.久久久| 亚洲av福利一区| 在线观看www视频免费| 国产av精品麻豆| 久久久精品94久久精品| 女人被躁到高潮嗷嗷叫费观| 成人毛片60女人毛片免费| xxxhd国产人妻xxx| 午夜免费观看性视频| 最黄视频免费看| 国产片特级美女逼逼视频| 国产在线视频一区二区| 精品一区在线观看国产| 久久久a久久爽久久v久久| 只有这里有精品99| av电影中文网址| 中国美白少妇内射xxxbb| 亚洲伊人久久精品综合| 亚洲成人一二三区av| 免费观看在线日韩| 久久韩国三级中文字幕| 国产精品不卡视频一区二区| 欧美精品国产亚洲| av黄色大香蕉| 国产精品国产三级专区第一集| 18禁在线无遮挡免费观看视频| 午夜视频国产福利| 中国三级夫妇交换| 免费看av在线观看网站| 1024视频免费在线观看| 一本久久精品| 国产欧美亚洲国产| 欧美性感艳星| 成人午夜精彩视频在线观看| 日本wwww免费看| 最近的中文字幕免费完整| 亚洲成色77777| 伦理电影免费视频| 精品国产一区二区三区四区第35| 美女国产视频在线观看| 纯流量卡能插随身wifi吗| 久久99一区二区三区| 啦啦啦中文免费视频观看日本| 精品一区二区三卡| 色婷婷av一区二区三区视频| 999精品在线视频| 欧美xxⅹ黑人| 久久精品国产综合久久久 | 少妇的逼好多水| 国产欧美日韩一区二区三区在线| 成人毛片60女人毛片免费| 国产成人a∨麻豆精品| 国产又色又爽无遮挡免| 国产男女内射视频| 欧美精品高潮呻吟av久久| 如日韩欧美国产精品一区二区三区| 国产亚洲精品久久久com| 日本黄色日本黄色录像| 狂野欧美激情性xxxx在线观看| 伊人久久国产一区二区| 五月伊人婷婷丁香| 久久精品久久久久久久性| 欧美xxⅹ黑人| 美国免费a级毛片| 久久久久精品人妻al黑| 国产高清国产精品国产三级| 久久精品熟女亚洲av麻豆精品| 26uuu在线亚洲综合色| 97人妻天天添夜夜摸| 90打野战视频偷拍视频| 日日爽夜夜爽网站| 欧美+日韩+精品| 亚洲内射少妇av| 十八禁高潮呻吟视频| 一区二区三区精品91| 免费日韩欧美在线观看| 99视频精品全部免费 在线| 免费大片18禁| 亚洲美女搞黄在线观看| 国产精品人妻久久久久久| 中文字幕免费在线视频6| 午夜免费男女啪啪视频观看| 国产一区二区三区综合在线观看 | 国产精品一国产av| 欧美性感艳星| 韩国高清视频一区二区三区| 免费黄网站久久成人精品| 亚洲国产日韩一区二区| 精品人妻偷拍中文字幕| 日韩制服骚丝袜av| 春色校园在线视频观看| 考比视频在线观看| 欧美精品国产亚洲| 又黄又爽又刺激的免费视频.| 精品国产乱码久久久久久小说| 久久人人爽人人片av| 亚洲欧美精品自产自拍| 美国免费a级毛片| 看免费成人av毛片| 精品少妇久久久久久888优播| 国产av一区二区精品久久| 亚洲av男天堂| 高清毛片免费看| 黄片无遮挡物在线观看| 99热这里只有是精品在线观看| 午夜福利网站1000一区二区三区| 国产亚洲精品久久久com| 国产白丝娇喘喷水9色精品| 国产日韩一区二区三区精品不卡| 午夜视频国产福利| √禁漫天堂资源中文www| 午夜福利视频在线观看免费| 日韩制服骚丝袜av| 老熟女久久久| 草草在线视频免费看| 婷婷成人精品国产| 午夜av观看不卡| 秋霞在线观看毛片| 国产精品 国内视频| 日韩欧美精品免费久久| 精品人妻在线不人妻| 欧美日韩精品成人综合77777| 日日爽夜夜爽网站| 久久精品国产鲁丝片午夜精品| 全区人妻精品视频| 9色porny在线观看| 丝袜喷水一区| 免费人妻精品一区二区三区视频| 一二三四中文在线观看免费高清| 自线自在国产av| av网站免费在线观看视频| 久久ye,这里只有精品| 国产日韩欧美视频二区| 欧美人与善性xxx| 婷婷色综合大香蕉| 王馨瑶露胸无遮挡在线观看| 国产免费又黄又爽又色| 男的添女的下面高潮视频| 观看美女的网站| 精品亚洲乱码少妇综合久久| 日本与韩国留学比较| 亚洲欧美一区二区三区国产| 亚洲av欧美aⅴ国产| 一级毛片 在线播放| 熟妇人妻不卡中文字幕| 视频区图区小说| 欧美成人精品欧美一级黄| 免费人妻精品一区二区三区视频| 51国产日韩欧美| 永久免费av网站大全| 2021少妇久久久久久久久久久| 久久久久久久久久人人人人人人| 9191精品国产免费久久| 高清视频免费观看一区二区| av免费在线看不卡| a级毛色黄片| 国产高清国产精品国产三级| 丝袜在线中文字幕| 国产精品一区二区在线观看99| 五月伊人婷婷丁香| 青青草视频在线视频观看| 又黄又爽又刺激的免费视频.| 久久99蜜桃精品久久| 毛片一级片免费看久久久久| 国产xxxxx性猛交| 哪个播放器可以免费观看大片| 中文字幕av电影在线播放| 狠狠精品人妻久久久久久综合| 精品国产一区二区三区久久久樱花| 少妇的丰满在线观看| 嫩草影院入口| 亚洲丝袜综合中文字幕| 亚洲色图 男人天堂 中文字幕 | 男人添女人高潮全过程视频| 九色成人免费人妻av| a级毛色黄片| kizo精华| 免费大片黄手机在线观看| 少妇熟女欧美另类| 一二三四中文在线观看免费高清| 久久久久久人妻| 女人久久www免费人成看片| 亚洲精品视频女| 亚洲内射少妇av| 午夜福利影视在线免费观看| 少妇人妻精品综合一区二区| av线在线观看网站| 黄色一级大片看看| 国产极品粉嫩免费观看在线| 人人澡人人妻人| 热re99久久精品国产66热6| 亚洲精华国产精华液的使用体验| 中文字幕另类日韩欧美亚洲嫩草| 中文欧美无线码| 在线观看美女被高潮喷水网站| 色婷婷久久久亚洲欧美| 老熟女久久久| 亚洲国产精品999| 欧美精品一区二区大全| 极品人妻少妇av视频| 亚洲国产色片| 久久女婷五月综合色啪小说| 搡老乐熟女国产| 三上悠亚av全集在线观看| 18在线观看网站| 男男h啪啪无遮挡| 乱码一卡2卡4卡精品| 一区在线观看完整版| 亚洲国产精品成人久久小说| 不卡视频在线观看欧美| 久久久久久伊人网av| 黄色 视频免费看| 国产成人精品无人区| 免费观看性生交大片5| 在线观看免费视频网站a站| 亚洲精品,欧美精品| 日韩 亚洲 欧美在线| 国产成人精品久久久久久| 亚洲综合精品二区| 伦理电影大哥的女人| 久久久久久久精品精品| 精品人妻一区二区三区麻豆| 国产极品粉嫩免费观看在线| 人妻系列 视频| 一级,二级,三级黄色视频| 国内精品宾馆在线| 精品国产露脸久久av麻豆| 王馨瑶露胸无遮挡在线观看| 精品酒店卫生间| 看免费av毛片| 亚洲欧美成人精品一区二区| 视频区图区小说| 国产成人欧美| 街头女战士在线观看网站| 侵犯人妻中文字幕一二三四区| 亚洲欧美一区二区三区黑人 | 桃花免费在线播放| 日韩制服丝袜自拍偷拍| 亚洲美女黄色视频免费看| 亚洲欧美一区二区三区国产| 精品熟女少妇av免费看| 国语对白做爰xxxⅹ性视频网站| 免费高清在线观看日韩| 男女午夜视频在线观看 | 十八禁高潮呻吟视频| 99热国产这里只有精品6| 在线观看www视频免费| 精品福利永久在线观看| 久久精品国产综合久久久 | 国产日韩一区二区三区精品不卡| 在线观看人妻少妇| 欧美人与性动交α欧美精品济南到 | 肉色欧美久久久久久久蜜桃| 日韩成人av中文字幕在线观看| 成人综合一区亚洲| 视频区图区小说| 一级黄片播放器| 男女边摸边吃奶| 一本久久精品| 插逼视频在线观看| 国产成人a∨麻豆精品| 岛国毛片在线播放| 成人影院久久| 亚洲丝袜综合中文字幕| 一级片'在线观看视频| 国产亚洲精品第一综合不卡 | 麻豆精品久久久久久蜜桃| videosex国产| 国产综合精华液| 亚洲精品一二三| 中文字幕最新亚洲高清| 蜜桃在线观看..| 丰满乱子伦码专区| 考比视频在线观看| 亚洲精品国产色婷婷电影| 久久精品熟女亚洲av麻豆精品| 亚洲国产欧美在线一区| 国产xxxxx性猛交| 国产欧美另类精品又又久久亚洲欧美| 亚洲精品一区蜜桃| 视频区图区小说| 韩国av在线不卡| 夫妻性生交免费视频一级片| 国产精品三级大全| 三上悠亚av全集在线观看| 男女免费视频国产| 青青草视频在线视频观看| 99久久人妻综合| 人成视频在线观看免费观看| 人人澡人人妻人| 女人被躁到高潮嗷嗷叫费观| 国产色婷婷99| 亚洲美女黄色视频免费看| 久久人人爽人人片av| 香蕉丝袜av| 春色校园在线视频观看| 欧美日韩综合久久久久久| 精品少妇久久久久久888优播| 中国美白少妇内射xxxbb| 国产精品一区二区在线不卡| 一区二区日韩欧美中文字幕 | 视频中文字幕在线观看| kizo精华| 精品久久国产蜜桃| 国产亚洲精品第一综合不卡 | 乱码一卡2卡4卡精品| 精品熟女少妇av免费看| 婷婷色综合www| 免费女性裸体啪啪无遮挡网站| 成年人午夜在线观看视频| 成人毛片60女人毛片免费| 亚洲国产色片| 天天躁夜夜躁狠狠久久av| 欧美xxxx性猛交bbbb| 男女边摸边吃奶| 少妇 在线观看| 最新的欧美精品一区二区| 亚洲人与动物交配视频| 永久网站在线| 99久久精品国产国产毛片| 寂寞人妻少妇视频99o| 欧美国产精品va在线观看不卡| 寂寞人妻少妇视频99o| 午夜激情av网站| 汤姆久久久久久久影院中文字幕| 51国产日韩欧美| 男女啪啪激烈高潮av片| 亚洲人与动物交配视频| 婷婷色麻豆天堂久久| 国产无遮挡羞羞视频在线观看| 亚洲色图 男人天堂 中文字幕 | 久久亚洲国产成人精品v| 久久综合国产亚洲精品| 亚洲久久久国产精品| 97在线视频观看| 亚洲精品久久午夜乱码| 久久亚洲国产成人精品v| 久久综合国产亚洲精品| 欧美激情国产日韩精品一区| 亚洲av电影在线观看一区二区三区| 少妇的逼水好多| 久久久精品94久久精品| 人妻 亚洲 视频| 日本与韩国留学比较| 精品一区二区免费观看| 亚洲av电影在线观看一区二区三区| 91成人精品电影| 各种免费的搞黄视频| 国产精品久久久久久久久免| 国产又爽黄色视频| 色婷婷av一区二区三区视频| 十八禁高潮呻吟视频| 久久久久视频综合| 成年女人在线观看亚洲视频| 久久久精品94久久精品| 亚洲精华国产精华液的使用体验| 久久久久人妻精品一区果冻| 熟女av电影| 一级毛片电影观看| 色5月婷婷丁香| 国产日韩一区二区三区精品不卡| kizo精华| 青青草视频在线视频观看| 黄色配什么色好看| 日韩成人av中文字幕在线观看| 国产在线视频一区二区| 精品少妇久久久久久888优播| 亚洲人成77777在线视频| 国产亚洲精品久久久com| 亚洲av中文av极速乱| 欧美 日韩 精品 国产| 精品一区二区免费观看| 日韩视频在线欧美| 国产片内射在线| 午夜免费观看性视频| www.熟女人妻精品国产 | av片东京热男人的天堂| 国产免费又黄又爽又色| 国产精品免费大片| 女人被躁到高潮嗷嗷叫费观| 观看av在线不卡| 黑人巨大精品欧美一区二区蜜桃 | 中文字幕亚洲精品专区| 观看av在线不卡| 九九爱精品视频在线观看| 极品人妻少妇av视频| 亚洲国产毛片av蜜桃av| 日韩大片免费观看网站| 久久精品熟女亚洲av麻豆精品| www.熟女人妻精品国产 | 视频中文字幕在线观看| 成人免费观看视频高清| 精品国产国语对白av| 天堂中文最新版在线下载| 三级国产精品片| 色吧在线观看| 亚洲精品美女久久av网站| 熟女人妻精品中文字幕| 极品人妻少妇av视频| 伦精品一区二区三区| 亚洲国产精品成人久久小说| 亚洲国产av影院在线观看| 人妻一区二区av| a级毛片黄视频| 亚洲伊人久久精品综合| 99re6热这里在线精品视频| 男人爽女人下面视频在线观看| 草草在线视频免费看| 色5月婷婷丁香| 曰老女人黄片| 亚洲av男天堂| 捣出白浆h1v1| 免费黄色在线免费观看| 爱豆传媒免费全集在线观看| 伦理电影大哥的女人| 人妻系列 视频| 又粗又硬又长又爽又黄的视频| 90打野战视频偷拍视频| 啦啦啦啦在线视频资源| 色哟哟·www| 香蕉丝袜av| 18禁观看日本| 亚洲av电影在线观看一区二区三区| 久久99一区二区三区| 精品亚洲乱码少妇综合久久| 久久久久久久精品精品|