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

    青藏高原納木錯(cuò)流域持久性有機(jī)污染物的多介質(zhì)遷移與歸趨模擬

    2017-10-13 03:06:27任嬌王小萍王傳飛龔平姚檀棟
    生態(tài)毒理學(xué)報(bào) 2017年3期
    關(guān)鍵詞:逸度納木錯(cuò)青藏高原

    任嬌,王小萍,王傳飛,龔平,姚檀棟

    1. 中國科學(xué)院青藏高原研究所,中國科學(xué)院青藏高原環(huán)境變化與地表過程實(shí)驗(yàn)室,北京 1001012. 中國科學(xué)院青藏高原地球科學(xué)卓越創(chuàng)新中心,北京 1001013. 中國科學(xué)院大學(xué),北京100049

    青藏高原納木錯(cuò)流域持久性有機(jī)污染物的多介質(zhì)遷移與歸趨模擬

    任嬌1,3,王小萍1,2,*,王傳飛1,2,龔平1,2,姚檀棟1,2

    1. 中國科學(xué)院青藏高原研究所,中國科學(xué)院青藏高原環(huán)境變化與地表過程實(shí)驗(yàn)室,北京 1001012. 中國科學(xué)院青藏高原地球科學(xué)卓越創(chuàng)新中心,北京 1001013. 中國科學(xué)院大學(xué),北京100049

    南亞排放的持久性有機(jī)污染物(POPs)可隨大氣傳輸?shù)角嗖馗咴?,然而POPs在高原多介質(zhì)間的遷移與分配尚不清晰。本研究利用三級(jí)逸度模型對(duì)4種POPs(六六六α-HCH,滴滴涕p,p′-DDT,菲Phe和苯并芘BaP)在納木錯(cuò)流域的遷移與歸趨進(jìn)行了模擬。結(jié)果表明,大氣沉降是該區(qū)域污染物的主要輸入過程,而降解損失則是主要的輸出途徑。就最終歸趨而言,土壤是POPs在納木錯(cuò)流域的重要儲(chǔ)庫,其存儲(chǔ)了大于50%的POPs。此外,湖水和沉積物分別對(duì)α-HCH和PAHs具有較強(qiáng)的存儲(chǔ)能力。靈敏度分析的結(jié)果表明,環(huán)境溫度、大氣中POPs的濃度及其理化性質(zhì)是影響POPs在環(huán)境中分布的關(guān)鍵參數(shù)。本研究明確了納木錯(cuò)流域不同POPs的遷移方向和歸趨特征,這將為青藏高原生態(tài)安全評(píng)估提供科學(xué)依據(jù)。

    青藏高原;持久性有機(jī)污染物;多介質(zhì)逸度模型;遷移;歸趨

    Received20 January 2017accepted13 March 2017

    Abstract: Persistent organic pollutants (POPs) emitted from South Asia can be transported to the Tibetan Plateau (TP), driven by favorable atmospheric circulation; however, their transfer between different media (air, soil, water, sediment etc.) and their final fate in the TP are unclear. In the present study, a level III fugacity model was applied to simulate the exchange fluxes and concentration distributions of four POPs (α-HCH, p,p′-DDT, Phe and BaP) in Nam Co Basin. The results showed that atmospheric deposition is the main input process for POPs into this region, and degradation is the main elimination pathway. Soil is the major sink of POPs, accounting for more than 50% of their total reserves. In addition, lake water and sediment have a certain storage capacity for α-HCH and polycyclic aromatic hydrocarbons, respectively. Sensitivity analysis showed that the concentrations of POPs in the air, their physicochemical properties, and air temperature, are key parameters affecting the environmental behavior and distribution of POPs. In conclusion, the direction of transfer, the fluxes amounts, and the final fate of different POPs among different media in Nam Co were successfully quantified in the model, thus contributing to a scientific basis for ecological safety assessments within the TP region.

    Keywords: Tibetan Plateau; persistent organic pollutants (POPs); multimedia fugacity model; transfer; fate

    持久性有機(jī)污染物(POPs)因其在環(huán)境中難以降解、對(duì)生物有著高毒性而受到廣泛關(guān)注[1]。由于具有半揮發(fā)性,POPs能夠隨大氣進(jìn)行長距離遷移[2-3],從而在全球范圍內(nèi)分布。即使在遠(yuǎn)離人類活動(dòng)的南極、北極和青藏高原也能監(jiān)測到POPs的蹤跡[4-6]。相比南北極,青藏高原毗鄰印度、巴基斯坦和尼泊爾等污染嚴(yán)重的南亞國家[7],因而更易受到污染源區(qū)的影響。目前,青藏高原的POPs研究已取得初步進(jìn)展,多項(xiàng)證據(jù)都表明南亞排放的有機(jī)氯農(nóng)藥(OCPs)[8-11]和多環(huán)芳烴(PAHs)[12-13]能夠在印度季風(fēng)和西風(fēng)的驅(qū)動(dòng)下傳輸至高原內(nèi)部。青藏高原的土壤[14-15]、草地[16]、冰雪[17-19]、湖水[20]、沉積物[21-22]和魚體[23-24]中都有POPs的相關(guān)報(bào)道,然而這些研究只能說明污染物到達(dá)了青藏高原。POPs在高原各環(huán)境介質(zhì)之間是如何遷移與分配的?其最終歸趨如何?顯然,目前的研究尚無法回答這些問題,只有將多介質(zhì)環(huán)境行為相結(jié)合才能全面認(rèn)識(shí)污染物在青藏高原的來龍去脈、辨析各環(huán)境介質(zhì)在POPs遷移與存儲(chǔ)中的作用,因此青藏高原POPs區(qū)域循環(huán)的系統(tǒng)研究亟待開展。

    青藏高原平均海拔在4 000 m以上,被稱為地球“第三極”,在此特殊地理單元中,污染物的遷移與循環(huán)也必然有其獨(dú)特之處。高原環(huán)境氣候寒冷,低溫會(huì)減弱地表揮發(fā),而加劇POPs的冷凝與沉降作用[25]。然而,近幾十年來,青藏高原正經(jīng)歷著顯著的變暖[26],山地冰川加速消融、徑流量增加[27],這可能會(huì)促進(jìn)污染物的釋放并向下游湖泊輸送。那么,在固有的低溫環(huán)境與全球升溫的耦合作用下,POPs在由大氣-土壤-冰川-湖水-沉積物所組成的區(qū)域系統(tǒng)中經(jīng)歷著怎樣的循環(huán)過程?各介質(zhì)對(duì)POPs的存儲(chǔ)起到多大作用?這些問題顯得尤為重要。

    在污染物多介質(zhì)循環(huán)的研究中,實(shí)測難度較大,而模型的運(yùn)用能夠簡化環(huán)境過程,達(dá)到全面評(píng)估污染物歸趨的目的。近年來,加拿大學(xué)者M(jìn)ackay[28-29]開發(fā)的逸度模型成為研究POPs環(huán)境行為的有力工具,用逸度來描述污染物在多環(huán)境介質(zhì)中的分配、遷移和歸趨是該模型的最大優(yōu)勢。穩(wěn)態(tài)的三級(jí)逸度模型已被運(yùn)用于北京[30]、天津[31-32]、蘭州[33]等地區(qū)研究POPs的跨界面環(huán)境過程,其結(jié)果與實(shí)測吻合較好。除區(qū)域的研究外,Wania和Mackay[34]首次將逸度模型應(yīng)用于全球尺度,揭示了POPs在全球環(huán)境中的遷移規(guī)律??梢?,逸度模型的適用性廣、預(yù)測效果好。因此,本研究擬借助多介質(zhì)逸度模型,選擇大氣、土壤、冰川與湖泊共存的納木錯(cuò)流域作為研究區(qū),模擬POPs在該區(qū)域的界面交換與遷移規(guī)律,以明確其最終歸宿。一方面,這將為正確評(píng)價(jià)青藏高原在POPs全球循環(huán)中的作用提供依據(jù);另一方面,這也是確保高原地區(qū)環(huán)境安全、生態(tài)系統(tǒng)可持續(xù)發(fā)展的重要前提。

    1 研究方法(Methods)

    1.1 研究區(qū)概況

    納木錯(cuò)(30°30′~30°56′N,90°16′~91°03′E, 見圖1)位于青藏高原中部、念青唐古拉山北麓,湖面海拔4 718 m,面積2 015 km2[35]。大氣降水和冰川融水是納木錯(cuò)主要的補(bǔ)給來源。近幾十年來,受氣溫升高和流域內(nèi)冰川加速消融的影響,納木錯(cuò)湖呈擴(kuò)張趨勢[36]。該地區(qū)年均氣溫為0°C,平均風(fēng)速4 m·s-1,年降雨量約400 mm左右,植被覆蓋以高寒草甸和高寒草原為主[36-37]。位于湖泊東南角的中國科學(xué)院納木錯(cuò)多圈層綜合觀測站(簡稱納木錯(cuò)站,圖1)是開展POPs多介質(zhì)研究的重要保障。

    圖1 納木錯(cuò)及采樣點(diǎn)的位置Fig. 1 Location of Nam Co Lake and sampling sites

    1.2 模型框架

    本研究以Mackay[29]的三級(jí)逸度模型為基礎(chǔ),結(jié)合納木錯(cuò)流域的具體環(huán)境構(gòu)建適合高原的穩(wěn)態(tài)多介質(zhì)模型。選取以納木錯(cuò)湖為中心,面積為10 680 km2的流域范圍作為研究區(qū)[38]。將環(huán)境系統(tǒng)分為4個(gè)主相:依次為大氣(1)、水(2)、土壤(3)和沉積物(4)。其中,各主相中又存在著氣、顆粒物、水和生物等子相(見表1)。

    模擬的主要環(huán)境過程包括POPs在大氣、土壤、湖水和沉積物各介質(zhì)之間的界面交換過程。遷移通量為逸度(f, Pa)與相間遷移系數(shù)D值(mol·h-1·Pa-1)的乘積。為考察冰川融水對(duì)污染物輸入的貢獻(xiàn),本研究將徑流作為向湖泊輸入的一個(gè)過程融入到模型中。由于缺乏上風(fēng)向大氣污染物輸入的準(zhǔn)確數(shù)據(jù),本文直接使用當(dāng)?shù)卮髿庵袑?shí)測POPs濃度和徑流向湖泊輸入的POPs通量作為驅(qū)動(dòng)數(shù)據(jù)。依據(jù)穩(wěn)態(tài)假設(shè),各環(huán)境相的質(zhì)量平衡方程可表達(dá)為:

    水:T02t+(D12d+D12r+D12p+D12w)×f1

    +D32×f3+D42×f4=(DR2+D21+D24+D2f)×f2

    土壤:(D13d+D13r+D13p+D13w)×f1

    =(DR3+D31+D32)×f3

    沉積物:D24×f2=(DR4+D42+D4B)×f4式中,T02t為污染物通過徑流向湖泊輸入的速率(mol·h-l);Dij為相間遷移速率,下標(biāo)d, r, p, w分別表示氣態(tài)交換、雨水溶解、顆粒態(tài)干沉降和顆粒態(tài)濕沉降;DRi為污染物在環(huán)境相i中的降解速率;D2f為魚對(duì)污染物的富集速率;D4B為污染物向深層沉積物的掩埋速率。各D值的具體計(jì)算方法參照文獻(xiàn)[29, 39]。

    1.3 參數(shù)識(shí)別

    基于納木錯(cuò)大氣POPs的觀測[40],本研究共選取4種化合物進(jìn)行模擬:六六六(α-HCH)和滴滴涕(p,p′-DDT)為OCPs的代表,主要來源于長距離傳輸;菲(Phe)是大氣中含量最高的PAHs[40],而苯并芘(BaP)的毒性較強(qiáng)且最受關(guān)注,因此將它們分別作為小分子和大分子PAH的代表。通過實(shí)測和文獻(xiàn)調(diào)研,共需收集2套數(shù)據(jù),分別用于模型輸入和模型驗(yàn)證。輸入數(shù)據(jù)包括流域環(huán)境參數(shù)(表1)、介質(zhì)遷移參數(shù)[29]、化合物理化性質(zhì)(表2)和模型驅(qū)動(dòng)數(shù)據(jù)(表2)。表1中各環(huán)境相的面積、深度、各子相的體積分?jǐn)?shù)和有機(jī)碳含量均為實(shí)測數(shù)據(jù)。表2中化合物的理化參數(shù)都經(jīng)過了溫度校正(T=273.15 K),大氣濃度來自于納木錯(cuò)站的觀測結(jié)果[40],徑流輸送量由河水中實(shí)測POPs的濃度與徑流量[38]相乘所得。用于模型驗(yàn)證的數(shù)據(jù)為各介質(zhì)中POPs的實(shí)測濃度值,它們是模擬結(jié)果準(zhǔn)確性評(píng)價(jià)和參數(shù)優(yōu)化的依據(jù)。為使實(shí)測數(shù)據(jù)具有較好的空間代表性,本研究環(huán)納木錯(cuò)湖采集了15個(gè)湖水、12個(gè)沉積物和20個(gè)土壤樣品(圖1)。樣品的前處理與實(shí)驗(yàn)室分析方法詳見文獻(xiàn)[15, 40-41]。表3列出了各介質(zhì)中POPs的實(shí)測濃度數(shù)據(jù)。

    1.4 求解與驗(yàn)證

    基于上述模型框架和參數(shù)收集,利用Matlab編程即可求解得到POPs在土壤、湖水和沉積物3個(gè)主環(huán)境相中的模擬濃度,將其與實(shí)測數(shù)據(jù)進(jìn)行比較,以驗(yàn)證模型的準(zhǔn)確性。為考察輸入?yún)?shù)對(duì)模型結(jié)果的影響,分別進(jìn)行靈敏度和不確定性分析。靈敏度分析旨在檢驗(yàn)單一參數(shù)變化對(duì)模擬結(jié)果的影響,其做法為分別對(duì)各輸入?yún)?shù)作±10%的變動(dòng),計(jì)算輸出結(jié)果Y的相對(duì)變化,常用靈敏度因子S表示[39]:

    S=

    通過比較S的大小可明確對(duì)模型結(jié)果影響較大的關(guān)鍵參數(shù)。在此基礎(chǔ)上,選取S>0.5的高靈敏度參數(shù)進(jìn)行不確定性分析,即利用蒙特卡洛方法考察多參數(shù)同時(shí)變化對(duì)模擬結(jié)果的影響[30]。假設(shè)輸入?yún)?shù)均符合對(duì)數(shù)正態(tài)分布,根據(jù)其均值和標(biāo)準(zhǔn)差產(chǎn)生隨機(jī)數(shù),運(yùn)行10 000次,可得模擬濃度的波動(dòng)范圍和變異系數(shù)(CV)。

    2 結(jié)果與討論(Results and discussion)

    2.1 模型驗(yàn)證

    經(jīng)過模擬計(jì)算,獲得了4種POPs在各子相中的模擬濃度,如圖2所示。為驗(yàn)證模型的準(zhǔn)確性,對(duì)模擬結(jié)果與實(shí)測值進(jìn)行了比較,發(fā)現(xiàn)各化合物在各子相的模擬與實(shí)測濃度都吻合較好,差異小于3倍(圖2)。這表明本研究構(gòu)建的多介質(zhì)逸度模型能夠較為準(zhǔn)確地反映當(dāng)?shù)氐奈廴厩闆r,適用于納木錯(cuò)地區(qū)POPs的遷移與歸趨模擬。

    2.2 多介質(zhì)遷移與最終歸趨

    基于模擬濃度可進(jìn)一步計(jì)算得到污染物在介質(zhì)間的遷移通量和各介質(zhì)中的儲(chǔ)量,其結(jié)果如圖3所示。從交換方向上看,α-HCH, p,p′-DDT, Phe和BaP都呈現(xiàn)出大氣向土壤、大氣向湖水、湖水向沉積物的凈沉降,這表明納木錯(cuò)的土壤和湖泊是大氣中污染物的匯。其中,以大氣向土壤的遷移通量最大,4種化合物的氣-土交換通量分別為1.1×10-4, 2.0×10-4, 2.0×10-1和3.9×10-3kg·h-1(圖3)。就湖水-沉積物界面而言,α-HCH和p,p′-DDT的沉降和釋放通量相差較小,而Phe和BaP向沉積物的積累通量遠(yuǎn)大于釋放通量(圖3),這可能會(huì)導(dǎo)致PAHs更多地存儲(chǔ)于湖泊沉積物中。除介質(zhì)間遷移外,降解損失為各介質(zhì)中污染物的主要輸出途徑(圖3)。

    表2 化合物的理化性質(zhì)與模型驅(qū)動(dòng)參數(shù)Table 2 Physicochemical properties and input data of the simulated chemicals

    表3 POPs在納木錯(cuò)各介質(zhì)中的實(shí)測濃度Table 3 Measured concentrations of POPs in the multimedia environment of Nam Co

    注:C22、C23、C2f、C33和C43分別表示污染物在湖水、湖水中懸浮顆粒物、魚體、土壤和沉積物中的濃度。ww為濕重;dw為干重。

    Note: C22, C23, C2f, C33and C43are the concentrations in lake water, suspended particles, fish, soil and sediment, respectively. ww: wet weight; dw: dry weight.

    圖2 模擬濃度與實(shí)測濃度的比較Fig. 2 Comparison between the predicted and measured concentrations

    近年來,多項(xiàng)研究都指出冰川消融已成為一種重要的“二次源”在向高山湖泊輸入POPs[42-43],然而支持這一觀點(diǎn)的實(shí)測數(shù)據(jù)較為缺乏。青藏高原在20世紀(jì)以來氣候快速變暖,升溫率是全球同期平均值的2倍,達(dá)到每10年0.3~0.4°C[44]。這導(dǎo)致高原冰川快速消融,促進(jìn)了徑流量的增加,并可能引起污染物的釋放。為探明冰川消融對(duì)POPs輸出和循環(huán)的影響,本研究將冰川融水徑流作為向湖泊輸入的一個(gè)過程融入多介質(zhì)模擬中。其結(jié)果顯示,對(duì)于4種目標(biāo)化合物,徑流輸入的貢獻(xiàn)都很小,即使對(duì)于水溶性較好的α-HCH,也僅有3%是來自于徑流輸入。這說明大氣沉降仍是湖泊中污染物的主要貢獻(xiàn)者(圖3)。這與Nellier等[45]在歐洲阿爾卑斯山冰川補(bǔ)給湖的研究結(jié)果一致。

    由于來源和理化性質(zhì)的差異,4種模擬化合物在納木錯(cuò)流域的儲(chǔ)量和最終歸趨表現(xiàn)出不同規(guī)律。α-HCH和p,p′-DDT的總儲(chǔ)量分別為3.5 kg和28.4 kg,它們主要來源于大氣的長距離傳輸[40]。相比之下,由于PAHs存在本地源的貢獻(xiàn)[46],Phe和BaP在納木錯(cuò)環(huán)境中的儲(chǔ)量比OCPs高2~3個(gè)數(shù)量級(jí),分別為8.2 t和3.9 t。這表明PAHs是當(dāng)?shù)刂饕奈廴疚?。就最終歸趨而言,湖水是α-HCH的歸宿之一,可占其總量的46%(圖3),這與α-HCH具有較好的水溶性有關(guān)。而疏水性強(qiáng)的p,p′-DDT不易在湖水中分配,大于88%的DDT都最終歸趨于土壤(圖3)。這說明DDT類化合物對(duì)青藏高原的陸地生態(tài)系統(tǒng)影響更大。同樣地,土壤也是PAHs類化合物的主要?dú)w宿,分別占Phe和BaP總儲(chǔ)量的49%和63%(圖3)。此外,PAHs在湖泊沉積物中也有一定的積累(約占35%,圖3),沉積物是其第二大儲(chǔ)庫,這可能對(duì)水生生態(tài)系統(tǒng)存在潛在風(fēng)險(xiǎn)。

    2.3 靈敏度分析

    靈敏度分析能夠反映模型輸出結(jié)果對(duì)參數(shù)變化的敏感程度,其結(jié)果對(duì)于模型優(yōu)化有著重要的指導(dǎo)意義。表4以S>0.5為標(biāo)準(zhǔn),列出了本研究中對(duì)模擬結(jié)果影響較大的高靈敏度參數(shù)。結(jié)果顯示,大氣中污染物的濃度(C11和C13)對(duì)所有化合物都有較大影響,S值接近于1,表現(xiàn)出大氣沉降對(duì)該地區(qū)污染物輸入的決定性作用。環(huán)境溫度(T)的S值為負(fù),且絕對(duì)值最大,這表明隨著溫度的降低,環(huán)境中POPs的濃度增加,該結(jié)果與POPs的冷捕集效應(yīng)一致[47]。由此可知,青藏高原溫度較低的地區(qū)必然富集了更多的污染物,而溫度的升高會(huì)使青藏高原“匯”的作用減弱。同時(shí),溫度的高靈敏度也證實(shí)了模擬中化合物理化性質(zhì)參數(shù)進(jìn)行溫度校正的必要性。

    如表4所示,化合物的亨利常數(shù)(H)、過冷飽和蒸汽壓(PL)和正辛醇-水分配系數(shù)(Kow)對(duì)模擬結(jié)果也有很大影響,部分S的絕對(duì)值遠(yuǎn)大于1。對(duì)介質(zhì)遷移類參數(shù)而言,僅氣-水界面遷移系數(shù)(K12)、干沉降速率(Kp)、降水速率(Kw)和濕沉降清除效率(Sr)的靈敏度系數(shù)較高。此外,還有部分參數(shù)的靈敏度具有專一性,即只對(duì)某一相的濃度有影響。例如,顆粒物中的有機(jī)碳含量(O23)僅對(duì)湖水懸浮顆粒物中POPs濃度有較大影響;土壤中降解速率(t3)和土壤相的厚度(h3)對(duì)土壤相濃度有較大影響;而沉積物中降解速率(t4)、沉積物的有機(jī)碳含量(O43)和水中顆粒物的沉降速率(KS)都只對(duì)沉積物相濃度有影響。

    可見,大氣中的濃度、環(huán)境溫度和化合物的性質(zhì)是影響POPs在環(huán)境中分布的關(guān)鍵因素,這些參數(shù)的準(zhǔn)確性是模擬結(jié)果可靠性的保證。

    2.4 不確定性分析

    模擬濃度的變異系數(shù)CV值越大,則表明其不確定性越大。圖4為本研究模型的不確定性分析結(jié)果,湖水和土壤中污染物濃度的CV值為100%左右,不確定較??;而沉積物的不確定性最大,CV為130%以上。這與影響沉積物濃度的參數(shù)較多、不確定性較大有關(guān)。與其他研究相比,本研究中模型的不確定性與北京地區(qū)的研究結(jié)果相當(dāng)(CV<150%)[30],而遠(yuǎn)小于蘭州地區(qū)多介質(zhì)模擬的誤差(CV達(dá)400%)[33]。這也證實(shí)了本研究中多介質(zhì)模型的準(zhǔn)確性。

    圖3 各環(huán)境介質(zhì)中污染物的儲(chǔ)量和介質(zhì)間的遷移通量Fig. 3 Mass balance diagram showing the transfer fluxes and reserves in different media of Nam Co Basin

    表4 參數(shù)靈敏度分析Table 4 Sensitivity analysis of the parameters in the model

    注: 表示參數(shù)的靈敏度因子小于0.5。

    Note: means the sensitivity coefficient (S) is lower than 0.5.

    圖4 模型模擬濃度的變異系數(shù)Fig. 4 The coefficient of variation (CV) of POPs concentrations in sub-phases

    綜上所述,多介質(zhì)環(huán)境逸度模型被成功應(yīng)用于青藏高原納木錯(cuò)流域,4種典型POPs在該流域大氣-土壤-湖泊-沉積物系統(tǒng)中的遷移和歸趨得以量化。研究結(jié)果顯示,大氣沉降是該地區(qū)污染物向土壤和湖泊輸入的主要過程,土壤和沉積物是POPs的最終歸宿。明確不同POPs的環(huán)境歸趨也為今后的研究指明了一定方向。首先,應(yīng)當(dāng)著重關(guān)注DDT沿土壤-草地-陸生生物的傳遞與風(fēng)險(xiǎn)評(píng)價(jià)。其次,青藏高原水生生物中PAHs類化合物的研究相對(duì)缺乏,其含量、富集規(guī)律和健康風(fēng)險(xiǎn)都有待研究。作為對(duì)POPs在青藏高原歸趨模擬的初步探索,本研究證實(shí)了逸度模型在高原的適用性,未來應(yīng)在此基礎(chǔ)上開展更加全面和深入的模擬研究。一方面,應(yīng)當(dāng)考慮環(huán)境和遷移過程的季節(jié)性差異,著重提高模型的時(shí)空分辨率;另一方面,利用四級(jí)逸度模型的長期預(yù)測優(yōu)勢,開展未來氣候變化情景下污染物在青藏高原的遷移與歸趨變化具有重要意義。

    [1] 余剛, 周隆超, 黃俊, 等. 持久性有機(jī)污染物和《斯德哥爾摩公約》履約[J]. 環(huán)境保護(hù), 2010, 23: 13-15

    Yu G, Zhou L C, Huang J, et al. Persistent organic pollutants and Stockholm Convention[J]. Environmental Protection, 2010, 23: 13-15 (in Chinese)

    [2] Beyer A, Mackay D, Matthies M, et al. Assessing long-range transport potential of persistent organic pollutants [J]. Environmental Science & Technology, 2000, 34(4): 699-703

    [3] Wang R, Tao S, Wang B, et al. Sources and pathways of polycyclic aromatic hydrocarbons transported to Alert, the Canadian High Arctic [J]. Environmental Science & Technology, 2010, 44(3): 1017-1022

    [4] MacDonald R W, Barrie L A, Bidleman T F, et al. Contaminants in the Canadian Arctic: 5 years of progress in understanding sources, occurrence and pathways [J]. Science of the Total Environment, 2000, 254(2-3): 93-234

    [5] Kallenborn R, Breivik K, Eckhardt S, et al. Long-term monitoring of persistent organic pollutants (POPs) at the Norwegian Troll station in Dronning Maud Land, Antarctica [J]. Atmospheric Chemistry and Physics, 2013, 13(14): 6983-6992

    [6] Wang X P, Gong P, Yao T D, et al. Passive air sampling of organochlorine pesticides, polychlorinated biphenyls, and polybrominated diphenyl ethers across the Tibetan Plateau [J]. Environmental Science & Technology, 2010, 44(8): 2988-2993

    [7] Ali U, Syed J H, Malik R N, et al. Organochlorine pesticides (OCPs) in South Asian region: A review [J]. Science of the Total Environment, 2014, 476: 705-717

    [8] Sheng J J, Wang X P, Gong P, et al. Monsoon-driven transport of organochlorine pesticides and polychlorinated biphenyls to the Tibetan Plateau: Three year atmospheric monitoring study [J]. Environmental Science & Technology, 2013, 47(7): 3199-3208

    [9] Gong P, Wang X P, Sheng J J, et al. Variations of organochlorine pesticides and polychlorinated biphenyls in atmosphere of the Tibetan Plateau: Role of the monsoon system [J]. Atmospheric Environment, 2010, 44(21-22): 2518-2523

    [10] Gong P, Wang X P, Xue Y G, et al. Influence of atmospheric circulation on the long-range transport of organochlorine pesticides to the western Tibetan Plateau [J]. Atmospheric Research, 2015, 166: 157-164

    [11] Wang X P, Ren J, Gong P, et al. Spatial distribution of the persistent organic pollutants across the Tibetan Plateau and its linkage with the climate systems: A 5-year air monitoring study [J]. Atmospheric Chemistry and Physics, 2016, 16(11): 6901-6911

    [12] Wang X P, Gong P, Sheng J J, et al. Long-range atmospheric transport of particulate polycyclic aromatic hydrocarbons and the incursion of aerosols to the southeast Tibetan Plateau [J]. Atmospheric Environment, 2015, 115: 124-131

    [13] Ma W L, Qi H, Baidron S, et al. Implications for long-range atmospheric transport of polycyclic aromatic hydrocarbons in Lhasa, China [J]. Environmental Science and Pollution Research, 2013, 20(8): 5525-5533

    [14] Tao S, Wang W T, Liu W X, et al. Polycyclic aromatic hydrocarbons and organochlorine pesticides in surface soils from the Qinghai-Tibetan Plateau [J]. Journal of Environmental Monitoring, 2011, 13(1): 175-181

    [15] Wang X P, Sheng J J, Gong P, et al. Persistent organic pollutants in the Tibetan surface soil: Spatial distribution, air-soil exchange and implications for global cycling [J]. Environmental Pollution, 2012, 170: 145-151

    [16] Wang C F, Wang X P, Yuan X H, et al. Organochlorine pesticides and polychlorinated biphenyls in air, grass and yak butter from Namco in the central Tibetan Plateau [J]. Environmental Pollution, 2015, 201: 50-57

    [17] Wang X P, Halsall C, Codling G, et al. Accumulation of perfluoroalkyl compounds in Tibetan mountain snow: Temporal patterns from 1980 to 2010 [J]. Environmental Science & Technology, 2014, 48(1): 173-181

    [18] Wang X P, Gong P, Zhang Q G, et al. Impact of climate fluctuations on deposition of DDT and hexachlorocyclohexane in mountain glaciers: Evidence from ice core records [J]. Environmental Pollution, 2010, 158(2): 375-380

    [19] Kang J H, Sungdeuk C, Hyokeun P, et al. Atmospheric deposition of persistent organic pollutants to the East Rongbuk Glacier in the Himalayas [J]. Science of the Total Environment, 2009, 408(1): 57-63

    [20] 張偉玲, 張干, 祁士華, 等. 西藏錯(cuò)鄂湖和羊卓雍湖水體及沉積物中有機(jī)氯農(nóng)藥的初步研究[J]. 地球化學(xué), 2003, 32(4): 363-367

    Zhang W L, Zhang G, Qi S H, et al. A preliminary study of organochlorinepesticides in water and sediments from two Tibetan Lakes[J]. Geochimica, 2003, 32(4): 363-367 (in Chinese)

    [21] Cheng H, Lin T, Zhang G, et al. DDTs and HCHs in sediment cores from the Tibetan Plateau [J]. Chemosphere, 2014, 94: 183-189

    [22] Yang R Q, Xie T, Li A, et al. Sedimentary records of polycyclic aromatic hydrocarbons (PAHs) in remote lakes across the Tibetan Plateau [J]. Environmental Pollution, 2016, 214: 1-7

    [23] Shi Y L, Pan Y Y, Yang R Q, et al. Occurrence of perfluorinated compounds in fish from Qinghai-Tibetan Plateau [J]. Environment International, 2010, 36(1): 46-50

    [24] Yang R Q, Wang Y W, Li A, et al. Organochlorine pesticides and PCBs in fish from lakes of the Tibetan Plateau and the implications [J]. Environmental Pollution, 2010, 158(6): 2310-2316

    [25] Daly G L,Wania F. Organic contaminants in mountains [J]. Environmental Science & Technology, 2005, 39(2): 385-398

    [26] Liu X, Chen B. Climatic warming in the Tibetan Plateau during recent decades [J]. International Journal of Climatology, 2000, 20(14): 1729-1742

    [27] Yang K, Wu H, Qin J, et al. Recent climate changes over the Tibetan Plateau and their impacts on energy and water cycle: A review [J]. Global and Planetary Change, 2014, 112: 79-91

    [28] Mackay D. Finding fugacity feasible [J]. Environmental Science & Technology, 1979, 13(10): 1218-1223

    [29] Mackay D. Multimedia Environmental Models: The Fugacity Approach [M]. Second edition. Boca Raton, FL:Lewis Publishers, 2001: 1-183

    [30] 曹紅英, 梁濤, 陶澍. 北京地區(qū)有機(jī)氯農(nóng)藥的跨界面遷移與歸趨[J]. 應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào), 2004, 12(3): 249-258

    Cao H Y, Liang T, Tao S. Simulating the transfer flux and fate of organochlorine pesticide in Beijing [J]. Journal of Basic Science and Engineering, 2004, 12(3): 249-258 (in Chinese)

    [31] Tao S, Cao H Y, Liu W, et al. Fate modeling of phenanthrene with regional variation in Tianjin, China [J]. Environmental Science & Technology, 2003, 37(11): 2453-2459

    [32] 曹紅英, 龔鐘明, 曹軍, 等. 估算天津環(huán)境中γ-HCH歸宿的逸度模型[J]. 環(huán)境科學(xué), 2003, 24(2): 77-81

    Cao H Y, Gong Z M, Cao J, et al. Evaluating the fate of γ-HCH using fugacity model in Tianjin environment [J]. Environmental Science, 2003, 24(2): 77-81 (in Chinese)

    [33] 高宏, 董繼元, 吳軍年. 蘭州地區(qū)HCHs的跨界面遷移與歸趨[J]. 中國環(huán)境科學(xué), 2008, 28(5): 407-411

    Gao H, Dong J Y, Wu J N. Transfer and fate of HCHs in Lanzhou region [J]. China Environmental Science, 2008, 28(5): 407-411 (in Chinese)

    [34] Wania F, Mackay D. A global distribution model for persistent organic-chemicals [J]. Science of the Total Environment, 1995, 160-161: 211-232

    [35] 王君波, 朱立平, 鞠建廷, 等. 西藏納木錯(cuò)水深分布及現(xiàn)代湖沼學(xué)特征初步分析[J]. 湖泊科學(xué), 2009, 21(1): 128-134

    Wang J B, Zhu L P, Ju J T, et al. Bathymetric survey and modern limnological parameters of Nam Co, central Tibet [J].Journal of Lake Sciences, 2009, 21(1): 128-134(in Chinese)

    [36] 朱立平, 謝曼平, 吳艷紅. 西藏納木錯(cuò)1971~2004年湖泊面積變化及其原因的定量分析[J]. 科學(xué)通報(bào), 2010, 55(18): 1789-1798

    Zhu L P, Xie M P, Wu Y H. Quantitative analysis of lake area variations and the influence factors from 1971 to 2004 in the Nam Co Basin of the Tibetan Plateau. [J]. Chinese Science Bulletin, 2010, 55(18): 1789-1798(in Chinese)

    [37] 游慶龍, 康世昌, 李潮流, 等. 青藏高原納木錯(cuò)氣象要素變化特征[J]. 氣象, 2007, 33(3): 54-60

    You Q L, Kang S C, Li C L, et al. Variation features of meteorological elements at Namco station, Tibetan Plateau [J]. Meteorological Monthly, 2007, 33(3): 54-60 (in Chinese)

    [38] Zhou S Q, Kang S C, Chen F, et al. Water balance observations reveal significant subsurface water seepage from Lake Nam Co, south-central Tibetan Plateau [J]. Journal of Hydrology, 2013, 491: 89-99

    [39] Xu F L, Qin N, Zhu Y, et al. Multimedia fate modeling of polycyclic aromatic hydrocarbons (PAHs) in Lake Small Baiyangdian, Northern China [J]. Ecological Modelling, 2013, 252: 246-257

    [40] Ren J, Wang X P, Wang C F, et al. Atmospheric processes of organic pollutants over a remote lake of the central Tibetan Plateau: Implications for regional cycling [J]. Atmospheric Chemistry and Physics, 2017, 17: 1401-1415

    [41] Ren J, Wang X P, Wang C F, et al. Biomagnification of persistent organic pollutants along a high-altitude aquatic food chain in the Tibetan Plateau: Processes and mechanisms [J]. Environmental Pollution, 2017, 220: 636-643

    [42] Bogdal C, Schmid P, Zennegg M, et al. Blast from the past: Melting glaciers as a relevant source for persistent organic pollutants [J]. Environmental Science & Technology, 2009, 43(21): 8173-8177

    [43] Schmid P, Bogdal C, Blüthgen N, et al. The missing piece: Sediment records in remote mountain lakes confirm glaciers being secondary sources of persistent organic pollutants [J]. Environmental Science & Technology, 2010, 45(1): 203-208

    [44] 陳德亮, 徐柏青, 姚檀棟, 等. 青藏高原環(huán)境變化科學(xué)評(píng)估:過去、現(xiàn)在與未來[J]. 科學(xué)通報(bào), 2015, 60(32): 3025-3035

    Chen D L, Xu B Q, Yao T D, et al. Assessment of past, present and future environmental changes on the Tibetan Plateau[J]. Chinese Science Bulletin, 2015, 60(32): 3025-3035(in Chinese)

    [45] Nellier Y M, Perga M E, Cottin N, et al. Mass budget in two high altitude lakes reveals their role as atmospheric PCB sinks [J]. Science of the Total Environment, 2015, 511: 203-213

    [46] Li C L, Kang S C, Chen P F, et al. Characterizations of particle-bound trace metals and polycyclic aromatic hydrocarbons (PAHs) within Tibetan tents of south Tibetan Plateau, China [J]. Environmental Science and Pollution Research, 2012, 19(5): 1620-1628

    [47] Scheringer M, Wegmann F, Fenner K, et al. Investigation of the cold condensation of persistent organic pollutants with a global multimedia fate model [J]. Environmental Science & Technology, 2000, 34(9): 1842-1850

    MultimediaFateModelingofPersistentOrganicPollutantsinNamCoBasin,TibetanPlateau

    Ren Jiao1,3, Wang Xiaoping1,2,*, Wang Chuanfei1,2, Gong Ping1,2, Yao Tandong1,2

    1. Key Laboratory of Tibetan Environment Changes and Land Surface Processes, Institute of Tibetan Plateau Research, Chinese Academy of Sciences, Beijing 100101, China2. CAS Center for Excellence in Tibetan Plateau Earth Sciences, Beijing 100101, China3. University of Chinese Academy of Sciences, Beijing 100049, China

    10.7524/AJE.1673-5897.20170120002

    2017-01-20錄用日期2017-03-13

    1673-5897(2017)3-170-10

    X171.5

    A

    王小萍(1976-),女,博士,研究員,主要研究方向?yàn)榍嗖馗咴h(huán)境污染與變化。

    國家自然科學(xué)基金項(xiàng)目(41222010, 41671480);中國科學(xué)院青年創(chuàng)新促進(jìn)會(huì)項(xiàng)目(CAS2011067)

    任嬌(1990-),女,博士研究生,研究方向?yàn)榍嗖馗咴h(huán)境污染,E-mail: renjiao@itpcas.ac.cn

    *通訊作者(Corresponding author), E-mail: wangxp@itpcas.ac.cn

    任嬌, 王小萍, 王傳飛, 等. 青藏高原納木錯(cuò)流域持久性有機(jī)污染物的多介質(zhì)遷移與歸趨模擬[J]. 生態(tài)毒理學(xué)報(bào),2017, 12(3): 170-179

    Ren J, Wang X P, Wang C F, et al. Multimedia fate modeling of persistent organic pollutants in Nam Co Basin, Tibetan Plateau [J]. Asian Journal of Ecotoxicology, 2017, 12(3): 170-179 (in Chinese)

    猜你喜歡
    逸度納木錯(cuò)青藏高原
    納木錯(cuò)
    青藏高原上的“含羞花”
    地幔氧逸度的時(shí)空變化
    輕輕松松聊漢語 納木錯(cuò)湖
    金橋(2020年9期)2020-10-27 02:00:06
    不同溫度、壓強(qiáng)、氧逸度條件下斜方輝石含水性的實(shí)驗(yàn)研究
    為了讓青藏高原的天更藍(lán)、水更綠、草原更美
    幔源巖漿氧化還原狀態(tài)及對(duì)巖漿礦床成礦的制約*
    納木錯(cuò)
    納木錯(cuò)美麗的天仙
    青年歌聲(2017年5期)2017-03-15 01:21:48
    氣體混合爐中氧逸度控制
    午夜精品国产一区二区电影| 黄片大片在线免费观看| 桃色一区二区三区在线观看| 免费在线观看日本一区| 亚洲精品国产区一区二| 久久天堂一区二区三区四区| 亚洲欧美一区二区三区黑人| 熟女少妇亚洲综合色aaa.| 午夜福利影视在线免费观看| 每晚都被弄得嗷嗷叫到高潮| 天堂动漫精品| 别揉我奶头~嗯~啊~动态视频| 久久久久精品国产欧美久久久| 亚洲va日本ⅴa欧美va伊人久久| 午夜视频精品福利| 欧美激情极品国产一区二区三区| 日本一区二区免费在线视频| 久久久久久久久中文| 妹子高潮喷水视频| 国产精品二区激情视频| 电影成人av| 两人在一起打扑克的视频| 看片在线看免费视频| 国产一区二区三区视频了| 一进一出好大好爽视频| 啪啪无遮挡十八禁网站| 多毛熟女@视频| 日韩免费高清中文字幕av| 欧美一区二区精品小视频在线| 欧美老熟妇乱子伦牲交| 欧美成狂野欧美在线观看| 黄色丝袜av网址大全| 亚洲国产精品合色在线| 巨乳人妻的诱惑在线观看| 亚洲性夜色夜夜综合| 国产成人av教育| 人妻丰满熟妇av一区二区三区| 中文亚洲av片在线观看爽| 日韩精品免费视频一区二区三区| 亚洲国产欧美一区二区综合| www.精华液| 久久久久国产精品人妻aⅴ院| 在线视频色国产色| 久久亚洲真实| 91大片在线观看| 侵犯人妻中文字幕一二三四区| 午夜a级毛片| 高清黄色对白视频在线免费看| 国产成人欧美在线观看| a级毛片黄视频| 日本欧美视频一区| 99国产精品一区二区蜜桃av| 高清av免费在线| 女性被躁到高潮视频| 国产成人系列免费观看| 国产精品 国内视频| 黑人猛操日本美女一级片| 成年版毛片免费区| 欧美激情 高清一区二区三区| 一边摸一边抽搐一进一小说| 亚洲欧美日韩高清在线视频| 99久久人妻综合| 久久久精品欧美日韩精品| 国产色视频综合| 成熟少妇高潮喷水视频| 欧洲精品卡2卡3卡4卡5卡区| 一级作爱视频免费观看| 欧美乱妇无乱码| www.www免费av| 免费人成视频x8x8入口观看| 亚洲精品粉嫩美女一区| 国产精品九九99| av福利片在线| 91在线观看av| 日韩精品免费视频一区二区三区| 中文欧美无线码| 91老司机精品| 免费在线观看影片大全网站| 天天添夜夜摸| e午夜精品久久久久久久| av电影中文网址| 村上凉子中文字幕在线| 亚洲国产欧美一区二区综合| 久久青草综合色| 99久久久亚洲精品蜜臀av| 啦啦啦 在线观看视频| av福利片在线| 不卡av一区二区三区| 一a级毛片在线观看| 欧美在线一区亚洲| 精品久久久久久久久久免费视频 | 亚洲片人在线观看| av网站免费在线观看视频| 免费av中文字幕在线| 亚洲专区国产一区二区| 国产一区二区激情短视频| 三级毛片av免费| 黑人操中国人逼视频| 亚洲自拍偷在线| 亚洲av第一区精品v没综合| 色综合欧美亚洲国产小说| 97人妻天天添夜夜摸| 国产精品久久久人人做人人爽| 这个男人来自地球电影免费观看| 黄色丝袜av网址大全| 精品国产亚洲在线| 俄罗斯特黄特色一大片| 亚洲精品国产色婷婷电影| 日本黄色视频三级网站网址| 91成年电影在线观看| 久久久国产成人精品二区 | 亚洲五月天丁香| 中文字幕最新亚洲高清| 国产深夜福利视频在线观看| 久久亚洲真实| 脱女人内裤的视频| 成人三级黄色视频| 国产精品偷伦视频观看了| 怎么达到女性高潮| 在线视频色国产色| 久久天躁狠狠躁夜夜2o2o| 丰满迷人的少妇在线观看| 少妇粗大呻吟视频| 日日干狠狠操夜夜爽| 欧洲精品卡2卡3卡4卡5卡区| 超碰97精品在线观看| 久久久久久大精品| 亚洲精品国产色婷婷电影| 丝袜在线中文字幕| 性色av乱码一区二区三区2| 法律面前人人平等表现在哪些方面| 变态另类成人亚洲欧美熟女 | 两性夫妻黄色片| 日日爽夜夜爽网站| 日韩 欧美 亚洲 中文字幕| 一级,二级,三级黄色视频| 黑人欧美特级aaaaaa片| 成人国产一区最新在线观看| 日韩视频一区二区在线观看| 国产片内射在线| 亚洲精品国产区一区二| 免费女性裸体啪啪无遮挡网站| 午夜久久久在线观看| 欧美成人午夜精品| 黄色视频不卡| 欧美另类亚洲清纯唯美| 丝袜在线中文字幕| 日本a在线网址| 伦理电影免费视频| 久久久国产欧美日韩av| 后天国语完整版免费观看| 精品卡一卡二卡四卡免费| 欧美在线黄色| 后天国语完整版免费观看| 午夜91福利影院| 精品一品国产午夜福利视频| 精品一区二区三区视频在线观看免费 | 国内毛片毛片毛片毛片毛片| 国产在线精品亚洲第一网站| 搡老岳熟女国产| 久久精品91无色码中文字幕| 午夜日韩欧美国产| 国产高清激情床上av| 成人手机av| 午夜日韩欧美国产| 国产91精品成人一区二区三区| 欧美日韩福利视频一区二区| 色老头精品视频在线观看| 黄片小视频在线播放| 一进一出抽搐动态| 亚洲国产精品合色在线| 国产一区二区激情短视频| 久久久久国产一级毛片高清牌| 国产亚洲精品久久久久久毛片| 五月开心婷婷网| 成熟少妇高潮喷水视频| 国产深夜福利视频在线观看| 亚洲欧美一区二区三区黑人| 日日爽夜夜爽网站| 国产不卡一卡二| 大香蕉久久成人网| 国产精品久久久久成人av| 国产欧美日韩一区二区三| 别揉我奶头~嗯~啊~动态视频| 老司机深夜福利视频在线观看| 欧美午夜高清在线| x7x7x7水蜜桃| 久久久久久免费高清国产稀缺| 极品教师在线免费播放| 最近最新中文字幕大全电影3 | 看片在线看免费视频| 在线视频色国产色| 别揉我奶头~嗯~啊~动态视频| 不卡av一区二区三区| 88av欧美| 精品欧美一区二区三区在线| 成人国产一区最新在线观看| 久久久久九九精品影院| 中文字幕精品免费在线观看视频| 亚洲精品国产精品久久久不卡| 亚洲第一欧美日韩一区二区三区| 午夜视频精品福利| 18禁裸乳无遮挡免费网站照片 | 变态另类成人亚洲欧美熟女 | 亚洲精品在线美女| 伊人久久大香线蕉亚洲五| 在线观看免费视频网站a站| 国产一区二区在线av高清观看| 成人av一区二区三区在线看| 宅男免费午夜| 最近最新免费中文字幕在线| 两个人看的免费小视频| 每晚都被弄得嗷嗷叫到高潮| 成年人免费黄色播放视频| 在线观看午夜福利视频| 久久国产精品人妻蜜桃| 妹子高潮喷水视频| 国产黄a三级三级三级人| 人人妻人人澡人人看| 国产激情欧美一区二区| 丰满饥渴人妻一区二区三| 久久精品影院6| 国产成人av教育| 好看av亚洲va欧美ⅴa在| 啦啦啦免费观看视频1| 国产精品 国内视频| 天堂影院成人在线观看| 精品少妇一区二区三区视频日本电影| 少妇被粗大的猛进出69影院| 中亚洲国语对白在线视频| 国产日韩一区二区三区精品不卡| 999久久久精品免费观看国产| 男女床上黄色一级片免费看| 中文字幕人妻熟女乱码| 麻豆一二三区av精品| 三上悠亚av全集在线观看| 在线免费观看的www视频| 18禁裸乳无遮挡免费网站照片 | 成人国产一区最新在线观看| 成人亚洲精品一区在线观看| 无人区码免费观看不卡| 国产亚洲精品综合一区在线观看 | 深夜精品福利| 十八禁人妻一区二区| 国产在线精品亚洲第一网站| 身体一侧抽搐| 亚洲精品美女久久av网站| 久99久视频精品免费| 99国产精品99久久久久| 婷婷精品国产亚洲av在线| 午夜福利在线免费观看网站| x7x7x7水蜜桃| 午夜福利影视在线免费观看| 午夜影院日韩av| 欧美日韩亚洲综合一区二区三区_| 最新在线观看一区二区三区| 精品久久久精品久久久| 精品久久久久久,| av网站免费在线观看视频| 精品第一国产精品| 免费在线观看影片大全网站| 19禁男女啪啪无遮挡网站| 在线观看日韩欧美| 中文字幕av电影在线播放| 99久久人妻综合| 日本a在线网址| 亚洲国产中文字幕在线视频| 老熟妇乱子伦视频在线观看| 80岁老熟妇乱子伦牲交| 十八禁人妻一区二区| 黄色a级毛片大全视频| 国产高清视频在线播放一区| 国产精品爽爽va在线观看网站 | 国产一区二区三区视频了| 国产三级黄色录像| 99久久精品国产亚洲精品| 一本大道久久a久久精品| 国产在线观看jvid| 18禁裸乳无遮挡免费网站照片 | 美女午夜性视频免费| 久久精品成人免费网站| 亚洲欧美日韩高清在线视频| 欧美+亚洲+日韩+国产| 神马国产精品三级电影在线观看 | 久久精品国产亚洲av高清一级| 免费久久久久久久精品成人欧美视频| 一级a爱视频在线免费观看| 一进一出好大好爽视频| 香蕉国产在线看| 欧美最黄视频在线播放免费 | 久久中文字幕一级| 亚洲一区高清亚洲精品| 国产精品国产高清国产av| 美女午夜性视频免费| 亚洲一区二区三区欧美精品| 亚洲 欧美一区二区三区| 啦啦啦免费观看视频1| 欧美乱妇无乱码| 亚洲午夜理论影院| 亚洲欧美日韩无卡精品| 欧美在线一区亚洲| 欧美日韩亚洲国产一区二区在线观看| 侵犯人妻中文字幕一二三四区| 亚洲精品美女久久av网站| 两人在一起打扑克的视频| 欧美日韩av久久| 99热只有精品国产| a在线观看视频网站| 91老司机精品| 99国产极品粉嫩在线观看| 69av精品久久久久久| 91老司机精品| 成人国产一区最新在线观看| 久久人妻熟女aⅴ| 夫妻午夜视频| 色哟哟哟哟哟哟| videosex国产| 老司机午夜福利在线观看视频| 一二三四在线观看免费中文在| 亚洲久久久国产精品| 99久久国产精品久久久| 国产精品乱码一区二三区的特点 | 亚洲精品粉嫩美女一区| 老司机午夜十八禁免费视频| 欧美日韩瑟瑟在线播放| 深夜精品福利| 桃色一区二区三区在线观看| 一二三四在线观看免费中文在| www.999成人在线观看| 亚洲人成伊人成综合网2020| 老熟妇乱子伦视频在线观看| 亚洲久久久国产精品| 色婷婷av一区二区三区视频| 手机成人av网站| 1024香蕉在线观看| 人人妻人人澡人人看| 老司机亚洲免费影院| 国产av一区二区精品久久| 亚洲色图 男人天堂 中文字幕| 在线观看舔阴道视频| 欧美激情 高清一区二区三区| 免费在线观看黄色视频的| 老司机深夜福利视频在线观看| 91麻豆精品激情在线观看国产 | 丝袜在线中文字幕| 天天添夜夜摸| 美女扒开内裤让男人捅视频| 日本黄色视频三级网站网址| 国产av精品麻豆| 中文字幕最新亚洲高清| 伊人久久大香线蕉亚洲五| 老鸭窝网址在线观看| 窝窝影院91人妻| 在线观看一区二区三区激情| 超色免费av| 中文字幕高清在线视频| 精品一品国产午夜福利视频| 国产欧美日韩一区二区精品| 日韩免费高清中文字幕av| 国产91精品成人一区二区三区| 日本vs欧美在线观看视频| 黑人巨大精品欧美一区二区mp4| 91成年电影在线观看| 级片在线观看| 悠悠久久av| 久久精品人人爽人人爽视色| 97超级碰碰碰精品色视频在线观看| 欧美人与性动交α欧美软件| av视频免费观看在线观看| 18禁黄网站禁片午夜丰满| 日韩精品中文字幕看吧| 亚洲精品美女久久av网站| 欧美成人免费av一区二区三区| 动漫黄色视频在线观看| 国产有黄有色有爽视频| 两个人看的免费小视频| 亚洲自拍偷在线| 欧美日韩视频精品一区| 高清欧美精品videossex| 日本精品一区二区三区蜜桃| 在线av久久热| 欧美日韩黄片免| 黄频高清免费视频| 看免费av毛片| 国产不卡一卡二| 无人区码免费观看不卡| 久久国产精品影院| 露出奶头的视频| 久久人人爽av亚洲精品天堂| 精品久久久久久,| 午夜精品在线福利| 757午夜福利合集在线观看| 国产黄a三级三级三级人| 国产成人啪精品午夜网站| 嫩草影视91久久| 人成视频在线观看免费观看| 一级毛片精品| 男女下面进入的视频免费午夜 | 欧美日韩瑟瑟在线播放| 女生性感内裤真人,穿戴方法视频| 日本三级黄在线观看| 国产av在哪里看| 欧美色视频一区免费| 曰老女人黄片| 国产成人精品在线电影| 色综合欧美亚洲国产小说| 国产高清视频在线播放一区| 国产精品 国内视频| 精品卡一卡二卡四卡免费| 他把我摸到了高潮在线观看| 亚洲免费av在线视频| 亚洲成人免费电影在线观看| 国产有黄有色有爽视频| 一级毛片女人18水好多| 波多野结衣av一区二区av| 欧美色视频一区免费| 亚洲欧美精品综合一区二区三区| 国产成人精品在线电影| 中文字幕最新亚洲高清| 久久久久久人人人人人| av超薄肉色丝袜交足视频| 国产伦一二天堂av在线观看| 在线免费观看的www视频| 99久久综合精品五月天人人| 亚洲国产欧美日韩在线播放| 亚洲黑人精品在线| 成人永久免费在线观看视频| 亚洲熟妇中文字幕五十中出 | 免费在线观看日本一区| 国产男靠女视频免费网站| 亚洲精华国产精华精| 热re99久久国产66热| 精品国产美女av久久久久小说| 色婷婷av一区二区三区视频| 亚洲精品在线美女| 日韩av在线大香蕉| 国产熟女xx| 午夜免费观看网址| 色综合欧美亚洲国产小说| 嫩草影视91久久| 国产色视频综合| 午夜激情av网站| 久久中文字幕一级| 精品一区二区三区av网在线观看| 一a级毛片在线观看| 美女高潮喷水抽搐中文字幕| 中文字幕人妻熟女乱码| 变态另类成人亚洲欧美熟女 | 中文字幕另类日韩欧美亚洲嫩草| 国产精品电影一区二区三区| 久久亚洲精品不卡| 成年女人毛片免费观看观看9| 国产高清视频在线播放一区| 在线观看午夜福利视频| 免费在线观看视频国产中文字幕亚洲| 国产一区二区激情短视频| 两个人看的免费小视频| 久久国产精品男人的天堂亚洲| 精品无人区乱码1区二区| 久久婷婷成人综合色麻豆| 老鸭窝网址在线观看| 99国产精品一区二区三区| 免费高清在线观看日韩| 高潮久久久久久久久久久不卡| 另类亚洲欧美激情| av网站在线播放免费| 波多野结衣av一区二区av| 99久久国产精品久久久| 免费在线观看影片大全网站| 久久精品成人免费网站| 岛国在线观看网站| 久久中文字幕人妻熟女| 中文字幕av电影在线播放| 岛国视频午夜一区免费看| 亚洲男人的天堂狠狠| 午夜福利欧美成人| 很黄的视频免费| 三上悠亚av全集在线观看| 亚洲一区二区三区欧美精品| 老司机亚洲免费影院| 日韩欧美免费精品| 精品国产国语对白av| 91成年电影在线观看| 国产精品电影一区二区三区| www.999成人在线观看| 亚洲精品国产色婷婷电影| 在线观看免费视频网站a站| 精品午夜福利视频在线观看一区| 涩涩av久久男人的天堂| 国产精品美女特级片免费视频播放器 | 叶爱在线成人免费视频播放| 美女福利国产在线| 免费高清视频大片| 精品一区二区三区av网在线观看| 他把我摸到了高潮在线观看| 夫妻午夜视频| 国产亚洲精品第一综合不卡| 一进一出抽搐动态| 午夜视频精品福利| 波多野结衣一区麻豆| 999久久久精品免费观看国产| 十八禁人妻一区二区| 久久国产精品影院| 视频在线观看一区二区三区| 欧美不卡视频在线免费观看 | 国产极品粉嫩免费观看在线| 不卡一级毛片| 一边摸一边抽搐一进一小说| 国产精品一区二区精品视频观看| 一夜夜www| 亚洲国产精品合色在线| 大型av网站在线播放| 欧美日韩瑟瑟在线播放| 国产免费现黄频在线看| 中文字幕av电影在线播放| 美女高潮到喷水免费观看| 99精品在免费线老司机午夜| 亚洲熟妇熟女久久| 亚洲在线自拍视频| 美女高潮喷水抽搐中文字幕| 成人18禁高潮啪啪吃奶动态图| 少妇被粗大的猛进出69影院| 法律面前人人平等表现在哪些方面| 亚洲片人在线观看| 午夜久久久在线观看| 黄色女人牲交| 91大片在线观看| 美女高潮喷水抽搐中文字幕| 淫妇啪啪啪对白视频| 国产精品二区激情视频| 人人妻人人添人人爽欧美一区卜| 亚洲男人天堂网一区| 久久人人爽av亚洲精品天堂| 岛国视频午夜一区免费看| 黄片播放在线免费| 国产一区二区三区综合在线观看| 精品国产乱子伦一区二区三区| 亚洲av成人av| 日日摸夜夜添夜夜添小说| 久久久久久久午夜电影 | 日韩三级视频一区二区三区| 窝窝影院91人妻| 美女高潮到喷水免费观看| 一区二区三区国产精品乱码| 大陆偷拍与自拍| 精品卡一卡二卡四卡免费| 免费高清在线观看日韩| 咕卡用的链子| 欧美激情高清一区二区三区| 欧美日韩瑟瑟在线播放| 夜夜爽天天搞| 中文亚洲av片在线观看爽| 亚洲色图 男人天堂 中文字幕| 国产高清视频在线播放一区| 一进一出抽搐动态| 久久欧美精品欧美久久欧美| 夜夜夜夜夜久久久久| 日本欧美视频一区| 国产精品久久电影中文字幕| av视频免费观看在线观看| 国产aⅴ精品一区二区三区波| 欧美黑人欧美精品刺激| 亚洲在线自拍视频| 国产三级黄色录像| 亚洲欧美一区二区三区久久| 51午夜福利影视在线观看| 色婷婷久久久亚洲欧美| 亚洲五月天丁香| 中文亚洲av片在线观看爽| 成年人黄色毛片网站| 欧美一区二区精品小视频在线| 免费在线观看亚洲国产| 99久久人妻综合| 亚洲第一欧美日韩一区二区三区| 日韩精品免费视频一区二区三区| 在线播放国产精品三级| 欧美黄色片欧美黄色片| 久久国产亚洲av麻豆专区| 嫩草影院精品99| 国产精品乱码一区二三区的特点 | 午夜福利在线观看吧| 成人特级黄色片久久久久久久| 男人的好看免费观看在线视频 | 国产欧美日韩一区二区三| 精品久久久久久久久久免费视频 | 热99re8久久精品国产| 女人高潮潮喷娇喘18禁视频| 中文字幕色久视频| 一区二区日韩欧美中文字幕| 久久影院123| 国产精品九九99| www日本在线高清视频| 人成视频在线观看免费观看| 日本黄色视频三级网站网址| 在线观看免费午夜福利视频| 99久久人妻综合| 成人三级做爰电影| 丰满的人妻完整版| 亚洲 欧美一区二区三区| 久久中文字幕一级| 曰老女人黄片| tocl精华| 成人国产一区最新在线观看| 亚洲成国产人片在线观看| 精品人妻在线不人妻| 午夜免费鲁丝| 国产精品自产拍在线观看55亚洲| 一边摸一边抽搐一进一出视频| 一区二区三区精品91| 精品电影一区二区在线| www国产在线视频色| 一进一出好大好爽视频| 久久国产亚洲av麻豆专区| 亚洲欧美一区二区三区黑人| 女同久久另类99精品国产91|