• <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在线观看网站| 婷婷色av中文字幕| 99热这里只有是精品50| 国产午夜精品一二区理论片| 高清视频免费观看一区二区| 搞女人的毛片| 99精国产麻豆久久婷婷| 在线看a的网站| 亚洲av国产av综合av卡| 国产精品国产三级专区第一集| 3wmmmm亚洲av在线观看| 大片电影免费在线观看免费| 日本av手机在线免费观看| 国模一区二区三区四区视频| 国产亚洲精品久久久com| 狂野欧美激情性bbbbbb| 在线观看一区二区三区激情| 国产欧美亚洲国产| 国产精品福利在线免费观看| 亚洲内射少妇av| 国产成人a∨麻豆精品| 熟女人妻精品中文字幕| 国产亚洲午夜精品一区二区久久 | 一级av片app| 国产精品麻豆人妻色哟哟久久| 看黄色毛片网站| av福利片在线观看| 中文欧美无线码| 亚洲精品色激情综合| 女人十人毛片免费观看3o分钟| 久久久久九九精品影院| 亚洲av不卡在线观看| 欧美潮喷喷水| 永久网站在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 天天躁夜夜躁狠狠久久av| 亚洲国产av新网站| 国产成人免费无遮挡视频| 国产成人aa在线观看| 亚洲不卡免费看| 中国美白少妇内射xxxbb| 寂寞人妻少妇视频99o| 久久99蜜桃精品久久| 欧美精品国产亚洲| 男插女下体视频免费在线播放| 久久精品国产亚洲网站| 亚洲人成网站高清观看| 精品国产一区二区三区久久久樱花 | 久久久国产一区二区| 国产成人免费无遮挡视频| 秋霞在线观看毛片| 高清在线视频一区二区三区| 三级经典国产精品| 哪个播放器可以免费观看大片| 少妇的逼好多水| 亚洲精品久久午夜乱码| 乱系列少妇在线播放| 精品人妻一区二区三区麻豆| 日韩在线高清观看一区二区三区| 久久99蜜桃精品久久| 一级毛片久久久久久久久女| 亚洲不卡免费看| 亚洲成色77777| 成人国产麻豆网| 六月丁香七月| 黄色欧美视频在线观看| 国产精品国产三级专区第一集| 久久ye,这里只有精品| 国产精品国产三级国产专区5o| 免费av观看视频| freevideosex欧美| 久久97久久精品| 国产精品久久久久久久久免| 波野结衣二区三区在线| xxx大片免费视频| 美女内射精品一级片tv| 精品久久久噜噜| 我的老师免费观看完整版| 亚洲av福利一区| 一区二区三区乱码不卡18| av黄色大香蕉| 国内少妇人妻偷人精品xxx网站| 色哟哟·www| 久久精品夜色国产| 天天一区二区日本电影三级| 亚洲伊人久久精品综合| 老司机影院成人| 成人一区二区视频在线观看| 久久99精品国语久久久| 在线观看免费高清a一片| 久久精品国产亚洲av天美| 岛国毛片在线播放| 免费不卡的大黄色大毛片视频在线观看| 美女cb高潮喷水在线观看| 国产一区有黄有色的免费视频| 免费观看在线日韩| 91久久精品国产一区二区三区| 国产中年淑女户外野战色| 欧美高清成人免费视频www| 日产精品乱码卡一卡2卡三| 激情 狠狠 欧美| 国产精品一二三区在线看| 精品少妇黑人巨大在线播放| 久久久久久伊人网av| 亚洲精品国产av成人精品| 天堂网av新在线| 午夜福利高清视频| 嫩草影院新地址| 亚洲久久久久久中文字幕| 波野结衣二区三区在线| 久久久久久久午夜电影| 欧美极品一区二区三区四区| 国产女主播在线喷水免费视频网站| 99久久中文字幕三级久久日本| 亚洲最大成人手机在线| 男女无遮挡免费网站观看| 男女啪啪激烈高潮av片| 丝袜美腿在线中文| 国产色婷婷99| 国产真实伦视频高清在线观看| 免费观看性生交大片5| 日本三级黄在线观看| 国产熟女欧美一区二区| 少妇的逼好多水| 国产乱来视频区| www.av在线官网国产| 亚洲伊人久久精品综合| 免费在线观看成人毛片| 岛国毛片在线播放| 又粗又硬又长又爽又黄的视频| 最近2019中文字幕mv第一页| 99久久九九国产精品国产免费| 国产伦精品一区二区三区视频9| 久久久久九九精品影院| 伦精品一区二区三区| av在线app专区| 精品久久久久久久人妻蜜臀av| 亚洲av中文av极速乱| 嫩草影院新地址| 97在线人人人人妻| 国内少妇人妻偷人精品xxx网站| 国产精品精品国产色婷婷| 女的被弄到高潮叫床怎么办| 美女脱内裤让男人舔精品视频| 中文精品一卡2卡3卡4更新| 国产黄色视频一区二区在线观看| 看黄色毛片网站| 五月伊人婷婷丁香| 亚洲欧美清纯卡通| 成人一区二区视频在线观看| 精品人妻视频免费看| 卡戴珊不雅视频在线播放| 精品久久久久久电影网| 秋霞在线观看毛片| 国产有黄有色有爽视频| 国产综合精华液| 亚洲精品乱码久久久v下载方式| 你懂的网址亚洲精品在线观看| 一本色道久久久久久精品综合| 国产精品久久久久久精品电影小说 | a级毛色黄片| 九草在线视频观看| 大片电影免费在线观看免费| 在线亚洲精品国产二区图片欧美 | 亚洲精品第二区| 十八禁网站网址无遮挡 | 日韩人妻高清精品专区| 好男人视频免费观看在线| 亚洲图色成人| 男女边摸边吃奶| 午夜免费鲁丝| 男女那种视频在线观看| 精品人妻视频免费看| 丰满少妇做爰视频| 人妻 亚洲 视频| 少妇人妻一区二区三区视频| 春色校园在线视频观看| 伦精品一区二区三区| 亚洲国产精品专区欧美| 国产成人aa在线观看| 久久鲁丝午夜福利片| 国产成人免费观看mmmm| 久久精品久久久久久久性| 国产精品一区二区在线观看99| 国产av不卡久久| 性色av一级| 成人亚洲欧美一区二区av| 国产成人福利小说| av线在线观看网站| 日产精品乱码卡一卡2卡三| 美女内射精品一级片tv| 麻豆成人av视频| 久久久成人免费电影| 国产片特级美女逼逼视频| 亚洲国产精品999| 亚洲怡红院男人天堂| 好男人在线观看高清免费视频| 尾随美女入室| 国产亚洲91精品色在线| 亚洲人成网站在线播| 亚洲av欧美aⅴ国产| 久久6这里有精品| 大又大粗又爽又黄少妇毛片口| 亚洲欧美日韩无卡精品| 国产淫片久久久久久久久| 日韩国内少妇激情av| 国产一区亚洲一区在线观看| 成人亚洲精品av一区二区| 三级国产精品欧美在线观看| 一个人观看的视频www高清免费观看| 青春草视频在线免费观看| 欧美日韩在线观看h| 18禁在线无遮挡免费观看视频| 日韩一区二区视频免费看| 婷婷色综合大香蕉| 午夜激情福利司机影院| 高清午夜精品一区二区三区| 国产av国产精品国产| 日韩电影二区| av免费在线看不卡| 欧美 日韩 精品 国产| 永久免费av网站大全| 日本av手机在线免费观看| 免费黄频网站在线观看国产| 国产真实伦视频高清在线观看| 毛片女人毛片| 九草在线视频观看| 久久精品综合一区二区三区| 国产一区有黄有色的免费视频| 欧美日韩视频精品一区| 亚洲欧美一区二区三区国产| 国产熟女欧美一区二区| 国产成人免费无遮挡视频| 国产伦精品一区二区三区四那| 一区二区av电影网| 日本黄色片子视频| 亚洲精品乱码久久久久久按摩| 性插视频无遮挡在线免费观看| 久久这里有精品视频免费| 免费播放大片免费观看视频在线观看| 热re99久久精品国产66热6| 嫩草影院新地址| 亚洲人与动物交配视频| 身体一侧抽搐| 嫩草影院入口| 国产伦在线观看视频一区| 亚洲av免费在线观看| 婷婷色av中文字幕| 汤姆久久久久久久影院中文字幕| 热99国产精品久久久久久7| 免费看光身美女| 2018国产大陆天天弄谢| 国产在线一区二区三区精| 国产极品天堂在线| 亚洲va在线va天堂va国产| 久久这里有精品视频免费| 免费av毛片视频| 亚洲欧美日韩另类电影网站 | 1000部很黄的大片| 欧美成人a在线观看| 国产精品国产av在线观看| 精品久久久久久电影网| 婷婷色av中文字幕| 天天躁夜夜躁狠狠久久av| 最后的刺客免费高清国语| 亚洲无线观看免费| 亚洲精品国产成人久久av| 亚洲国产精品国产精品| 亚洲成人中文字幕在线播放| 久久久久国产精品人妻一区二区| 免费观看av网站的网址| 国产精品久久久久久av不卡| 国产精品人妻久久久久久| 五月玫瑰六月丁香| 亚洲图色成人| 日韩不卡一区二区三区视频在线| a级毛色黄片| 国产免费福利视频在线观看| 婷婷色综合www| 亚洲最大成人中文| 777米奇影视久久| 亚洲精品国产av蜜桃| 亚洲激情五月婷婷啪啪| 亚洲欧美中文字幕日韩二区| 全区人妻精品视频| 亚洲精品影视一区二区三区av| 人人妻人人看人人澡| xxx大片免费视频| 久久久久久久久久久丰满| 亚洲欧美日韩无卡精品| 干丝袜人妻中文字幕| 欧美精品一区二区大全| 99久国产av精品国产电影| 伦理电影大哥的女人| av在线天堂中文字幕| 男人添女人高潮全过程视频| 一级毛片 在线播放| 欧美 日韩 精品 国产| 国产极品天堂在线| 少妇熟女欧美另类| 永久网站在线| 能在线免费看毛片的网站| 一个人看视频在线观看www免费| 丰满乱子伦码专区| 日韩一区二区视频免费看| 97超视频在线观看视频| 日日啪夜夜爽| 日本与韩国留学比较| 美女国产视频在线观看| 国产v大片淫在线免费观看| 久久久亚洲精品成人影院| 禁无遮挡网站| 欧美日韩综合久久久久久| 嫩草影院新地址| 日韩成人伦理影院| 国产黄a三级三级三级人| 国产大屁股一区二区在线视频| 日本色播在线视频| 亚洲精品456在线播放app| 蜜臀久久99精品久久宅男| 能在线免费看毛片的网站| 五月开心婷婷网| 综合色丁香网| 国产精品99久久久久久久久| 丰满乱子伦码专区| 99久久精品一区二区三区| 亚洲欧美成人综合另类久久久| 久久精品久久久久久久性| 日韩av在线免费看完整版不卡| 99久国产av精品国产电影| 熟女电影av网| 一本久久精品| 成年女人看的毛片在线观看| 国产亚洲av嫩草精品影院| 波野结衣二区三区在线| 丰满少妇做爰视频| 亚洲国产高清在线一区二区三| videossex国产| 少妇高潮的动态图| 少妇猛男粗大的猛烈进出视频 | 18禁裸乳无遮挡免费网站照片| 成人免费观看视频高清| 国产av国产精品国产| 日韩国内少妇激情av| 亚洲av一区综合| 91狼人影院| 日韩强制内射视频| 亚洲精华国产精华液的使用体验| 久久ye,这里只有精品| 80岁老熟妇乱子伦牲交| 噜噜噜噜噜久久久久久91| 成人欧美大片| 精品一区在线观看国产| 国产精品国产av在线观看| 一级毛片电影观看| 18+在线观看网站| 久久久久久国产a免费观看| 欧美潮喷喷水| 日韩强制内射视频| 欧美人与善性xxx| 青青草视频在线视频观看| 亚洲欧美成人精品一区二区| 久久综合国产亚洲精品| 天天一区二区日本电影三级| 国产免费视频播放在线视频| 亚洲色图av天堂| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产日韩一区二区| 免费观看性生交大片5| 亚洲在久久综合| 日本av手机在线免费观看| 日韩国内少妇激情av| 日本av手机在线免费观看| 国产探花极品一区二区| 国产黄片视频在线免费观看| 色综合色国产| 一区二区三区免费毛片| 成年女人看的毛片在线观看| 在线免费十八禁| 国产女主播在线喷水免费视频网站| 精品一区二区三区视频在线| 亚洲精品乱码久久久久久按摩| 久久精品夜色国产| 亚洲自拍偷在线| 超碰97精品在线观看| 国产精品无大码| 女人被狂操c到高潮| 精品人妻视频免费看| 如何舔出高潮| 亚洲国产精品专区欧美| tube8黄色片| 18禁裸乳无遮挡动漫免费视频 | 成人漫画全彩无遮挡| 97超碰精品成人国产| 我要看日韩黄色一级片| 国产一区二区三区综合在线观看 | 久久久久久久国产电影| 美女国产视频在线观看| 亚洲在线观看片| 2018国产大陆天天弄谢| 97超碰精品成人国产| 一级二级三级毛片免费看| 特级一级黄色大片| 毛片一级片免费看久久久久| av.在线天堂| 国产精品熟女久久久久浪| 边亲边吃奶的免费视频| 少妇人妻精品综合一区二区| 亚洲精品,欧美精品| 成人二区视频| 亚洲国产av新网站| 亚洲在线观看片| 色哟哟·www| 国产成人午夜福利电影在线观看| 日韩成人av中文字幕在线观看| 亚洲av国产av综合av卡| 嫩草影院入口| 日本黄大片高清| 亚洲欧美一区二区三区国产| 亚洲成人精品中文字幕电影| 亚洲人成网站在线播| 免费av毛片视频| 99热这里只有是精品50| 91久久精品国产一区二区三区| av在线亚洲专区| 一区二区三区四区激情视频| 国产精品一区二区在线观看99| 亚洲精品国产av成人精品| 我的老师免费观看完整版| 狂野欧美白嫩少妇大欣赏| 婷婷色麻豆天堂久久| 国产69精品久久久久777片| av免费观看日本| 亚洲av不卡在线观看| 看非洲黑人一级黄片| 别揉我奶头 嗯啊视频| 亚洲精品乱久久久久久| 赤兔流量卡办理| 欧美日韩视频高清一区二区三区二| 日韩,欧美,国产一区二区三区| 天天躁日日操中文字幕| 免费观看av网站的网址| 少妇人妻精品综合一区二区| 美女主播在线视频| 免费黄频网站在线观看国产| 欧美zozozo另类| av在线app专区| 欧美日韩精品成人综合77777| 久久久久久久国产电影| 免费电影在线观看免费观看| 国产成人a∨麻豆精品| 另类亚洲欧美激情| 精品少妇久久久久久888优播| 精品酒店卫生间| 性色avwww在线观看| 欧美精品一区二区大全| 国产亚洲av嫩草精品影院| 中文字幕av成人在线电影| 国产亚洲一区二区精品| 一区二区av电影网| 男人添女人高潮全过程视频| 一级二级三级毛片免费看| 国产精品.久久久| 最近手机中文字幕大全| 午夜日本视频在线| 18禁动态无遮挡网站| 精品一区二区免费观看| 国产探花在线观看一区二区| 美女xxoo啪啪120秒动态图| 成年版毛片免费区| 一本一本综合久久| 亚洲欧美精品专区久久| 亚洲国产成人一精品久久久| 亚州av有码| 国产伦精品一区二区三区四那| 超碰av人人做人人爽久久| 久久久久久久久久久免费av| 老司机影院成人| 久久久午夜欧美精品| 亚洲精品国产av蜜桃| 亚州av有码| 日韩在线高清观看一区二区三区| 毛片一级片免费看久久久久| 91久久精品电影网| 亚洲精品中文字幕在线视频 | 简卡轻食公司| 日本午夜av视频| 真实男女啪啪啪动态图| 日韩强制内射视频| 波野结衣二区三区在线| 国产 精品1| 国产亚洲5aaaaa淫片| 午夜精品国产一区二区电影 | 久久精品国产亚洲av天美| 久久久国产一区二区| 伦精品一区二区三区| 一区二区三区免费毛片| 日韩av在线免费看完整版不卡| 欧美日韩亚洲高清精品| 欧美精品国产亚洲| 免费看光身美女| 禁无遮挡网站| 国产午夜精品久久久久久一区二区三区| 欧美高清性xxxxhd video| 亚洲av福利一区| 日韩成人伦理影院| 一级毛片我不卡| 男女边摸边吃奶| 亚洲久久久久久中文字幕| 2018国产大陆天天弄谢| 免费大片黄手机在线观看| 国产精品熟女久久久久浪| 国产精品福利在线免费观看| 国产大屁股一区二区在线视频| 91久久精品国产一区二区三区| 国产精品.久久久| 在线免费十八禁| 精品久久久久久久久亚洲| 免费少妇av软件| 观看美女的网站| 免费观看在线日韩| 啦啦啦啦在线视频资源| 丰满人妻一区二区三区视频av| 国产精品人妻久久久久久| 精品久久久久久久末码| 久久久午夜欧美精品| 精品亚洲乱码少妇综合久久| 亚洲精品国产色婷婷电影| 黄色日韩在线| 国产黄色视频一区二区在线观看| 亚洲欧洲国产日韩| 美女被艹到高潮喷水动态| 日本爱情动作片www.在线观看| 国产成人精品久久久久久| 成人鲁丝片一二三区免费| 97超视频在线观看视频| 99re6热这里在线精品视频| 亚洲最大成人手机在线| 亚洲av中文av极速乱| 久久精品国产鲁丝片午夜精品| 丰满乱子伦码专区| 成人亚洲欧美一区二区av| 韩国av在线不卡| 七月丁香在线播放| 亚洲性久久影院| 日日摸夜夜添夜夜爱| 高清午夜精品一区二区三区| 精品一区二区三区视频在线| 亚洲av国产av综合av卡| 如何舔出高潮| 国产精品一及| 亚洲三级黄色毛片| 激情五月婷婷亚洲| 亚洲欧美日韩东京热| 国产国拍精品亚洲av在线观看| 欧美潮喷喷水| 综合色av麻豆| 国产伦精品一区二区三区四那| 欧美国产精品一级二级三级 | 亚洲av一区综合| 久久久久国产网址| 偷拍熟女少妇极品色| 中国美白少妇内射xxxbb| 色视频在线一区二区三区| 色综合色国产| 成人无遮挡网站| 美女脱内裤让男人舔精品视频| 国产精品熟女久久久久浪| 国产高潮美女av| 久久久久国产精品人妻一区二区| 天天躁夜夜躁狠狠久久av| 久久99热6这里只有精品| 国产午夜精品久久久久久一区二区三区| 在线亚洲精品国产二区图片欧美 | 成年人午夜在线观看视频| 欧美亚洲 丝袜 人妻 在线| 国产久久久一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91 | 国产老妇伦熟女老妇高清| 日产精品乱码卡一卡2卡三| 狠狠精品人妻久久久久久综合| 亚洲av男天堂| av专区在线播放| 日韩成人av中文字幕在线观看| 99精国产麻豆久久婷婷| 午夜爱爱视频在线播放| 国产 一区精品| av黄色大香蕉| 久久久精品欧美日韩精品| 精品国产露脸久久av麻豆| 啦啦啦啦在线视频资源| 日韩av免费高清视频| 国产毛片a区久久久久| 在线免费观看不下载黄p国产| 亚洲国产精品专区欧美| 在线精品无人区一区二区三 | 可以在线观看毛片的网站| 久久精品久久精品一区二区三区| 国产乱人视频| 亚洲天堂国产精品一区在线| 久久精品久久精品一区二区三区| 男女边摸边吃奶| 久久久a久久爽久久v久久| 久久6这里有精品| 亚洲精品自拍成人| av黄色大香蕉| 午夜亚洲福利在线播放| 亚洲成人av在线免费| 精品人妻熟女av久视频| 亚洲天堂av无毛| 国产色爽女视频免费观看| 国产精品女同一区二区软件| 国产老妇女一区| 亚洲欧美日韩另类电影网站 |