• <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
    氣體混合爐中氧逸度控制
    69av精品久久久久久| 日本 av在线| 中文亚洲av片在线观看爽| 在线十欧美十亚洲十日本专区| 国产激情偷乱视频一区二区| 级片在线观看| 久久青草综合色| 色尼玛亚洲综合影院| 中文字幕人妻熟女乱码| 这个男人来自地球电影免费观看| 男女之事视频高清在线观看| 亚洲人成77777在线视频| 亚洲第一青青草原| 久久婷婷人人爽人人干人人爱| 亚洲av成人不卡在线观看播放网| 国产视频一区二区在线看| 美女国产高潮福利片在线看| 国产伦在线观看视频一区| 免费观看精品视频网站| 日本三级黄在线观看| 99久久综合精品五月天人人| 国产av又大| 欧美日韩黄片免| 亚洲国产高清在线一区二区三 | 久久午夜综合久久蜜桃| 国产1区2区3区精品| 久久精品亚洲精品国产色婷小说| 757午夜福利合集在线观看| 国语自产精品视频在线第100页| 国产乱人伦免费视频| av电影中文网址| 91九色精品人成在线观看| 久久久久久久久中文| 亚洲午夜理论影院| 欧美乱色亚洲激情| 精品第一国产精品| 午夜成年电影在线免费观看| 日日爽夜夜爽网站| 人人妻人人看人人澡| 色综合婷婷激情| 国产久久久一区二区三区| 国产成人欧美在线观看| 一边摸一边抽搐一进一小说| 69av精品久久久久久| 女性生殖器流出的白浆| 久久香蕉国产精品| 亚洲精品久久成人aⅴ小说| 亚洲自拍偷在线| 久久婷婷人人爽人人干人人爱| 老熟妇仑乱视频hdxx| 午夜福利高清视频| 国产一区二区三区视频了| 亚洲国产精品合色在线| 国产激情偷乱视频一区二区| 午夜久久久久精精品| 成年免费大片在线观看| 99久久国产精品久久久| 欧美中文日本在线观看视频| 青草久久国产| 丝袜人妻中文字幕| 视频在线观看一区二区三区| 一级作爱视频免费观看| 麻豆国产av国片精品| 日韩视频一区二区在线观看| 国产精品一区二区三区四区久久 | 黄频高清免费视频| 日韩大码丰满熟妇| 日本撒尿小便嘘嘘汇集6| 制服人妻中文乱码| 日韩国内少妇激情av| 国语自产精品视频在线第100页| 一区二区三区激情视频| 亚洲午夜精品一区,二区,三区| 女同久久另类99精品国产91| 天天躁夜夜躁狠狠躁躁| 精华霜和精华液先用哪个| 中文亚洲av片在线观看爽| 一进一出抽搐gif免费好疼| 国产在线观看jvid| 日韩欧美一区视频在线观看| 欧美日韩中文字幕国产精品一区二区三区| 免费人成视频x8x8入口观看| 久99久视频精品免费| 一卡2卡三卡四卡精品乱码亚洲| 亚洲欧美激情综合另类| 久久精品成人免费网站| 午夜福利在线观看吧| av超薄肉色丝袜交足视频| 午夜福利视频1000在线观看| 亚洲精品久久国产高清桃花| 99精品久久久久人妻精品| 日韩高清综合在线| 俺也久久电影网| 亚洲精华国产精华精| 国产精品久久久久久人妻精品电影| 两个人视频免费观看高清| 久久精品成人免费网站| 国产精品精品国产色婷婷| 亚洲男人天堂网一区| 亚洲无线在线观看| 亚洲性夜色夜夜综合| 亚洲一码二码三码区别大吗| 在线观看免费日韩欧美大片| 精品熟女少妇八av免费久了| 日本a在线网址| 天天躁狠狠躁夜夜躁狠狠躁| 12—13女人毛片做爰片一| 久久精品亚洲精品国产色婷小说| 色综合婷婷激情| 最新美女视频免费是黄的| 男女那种视频在线观看| 老司机福利观看| 美女 人体艺术 gogo| 国产黄色小视频在线观看| 男女午夜视频在线观看| 亚洲五月色婷婷综合| 女性被躁到高潮视频| 色婷婷久久久亚洲欧美| 国产午夜福利久久久久久| 精品一区二区三区av网在线观看| 制服人妻中文乱码| 欧美国产日韩亚洲一区| 欧美另类亚洲清纯唯美| 久久久久国产一级毛片高清牌| 免费看日本二区| www国产在线视频色| 亚洲av电影不卡..在线观看| √禁漫天堂资源中文www| 午夜成年电影在线免费观看| 日本三级黄在线观看| 久久久久久久久免费视频了| 夜夜躁狠狠躁天天躁| 最近最新中文字幕大全电影3 | 观看免费一级毛片| 精品国产乱子伦一区二区三区| 国产激情欧美一区二区| 成人精品一区二区免费| 香蕉久久夜色| 亚洲午夜精品一区,二区,三区| 又紧又爽又黄一区二区| 国产一级毛片七仙女欲春2 | 在线观看日韩欧美| 日本撒尿小便嘘嘘汇集6| 色av中文字幕| 国产高清videossex| 久久 成人 亚洲| 麻豆一二三区av精品| www日本在线高清视频| 亚洲 欧美 日韩 在线 免费| 每晚都被弄得嗷嗷叫到高潮| 99国产极品粉嫩在线观看| avwww免费| 少妇被粗大的猛进出69影院| 日本精品一区二区三区蜜桃| or卡值多少钱| 久久精品国产亚洲av高清一级| 俺也久久电影网| 99热6这里只有精品| 久久久久久亚洲精品国产蜜桃av| 精品久久久久久久久久久久久 | 久久久久免费精品人妻一区二区 | 国产精品久久久av美女十八| 啦啦啦免费观看视频1| 欧美+亚洲+日韩+国产| 欧美日本亚洲视频在线播放| 国产蜜桃级精品一区二区三区| 国产欧美日韩精品亚洲av| 亚洲狠狠婷婷综合久久图片| 日日干狠狠操夜夜爽| 亚洲精品久久国产高清桃花| 亚洲一区高清亚洲精品| 两个人看的免费小视频| 中文字幕高清在线视频| 亚洲国产精品合色在线| 精品不卡国产一区二区三区| 12—13女人毛片做爰片一| 亚洲成人久久性| 美女大奶头视频| 黄片小视频在线播放| 国产私拍福利视频在线观看| 88av欧美| 长腿黑丝高跟| 亚洲精品美女久久久久99蜜臀| 深夜精品福利| 国产午夜精品久久久久久| 俄罗斯特黄特色一大片| www.www免费av| 后天国语完整版免费观看| 成人手机av| 可以在线观看的亚洲视频| 91麻豆av在线| 一本精品99久久精品77| 亚洲欧美日韩高清在线视频| 久久精品亚洲精品国产色婷小说| 91字幕亚洲| 亚洲专区中文字幕在线| 亚洲欧美日韩无卡精品| 久久亚洲精品不卡| 91成年电影在线观看| 亚洲激情在线av| 十八禁网站免费在线| 久久久水蜜桃国产精品网| 黄频高清免费视频| 别揉我奶头~嗯~啊~动态视频| 成年免费大片在线观看| 成年女人毛片免费观看观看9| e午夜精品久久久久久久| 国产成人啪精品午夜网站| 亚洲 欧美一区二区三区| 国产熟女xx| 麻豆久久精品国产亚洲av| 女人高潮潮喷娇喘18禁视频| 国产精品综合久久久久久久免费| 无限看片的www在线观看| 超碰成人久久| 中文字幕精品免费在线观看视频| 欧美激情久久久久久爽电影| 一边摸一边做爽爽视频免费| 国产成人欧美| 波多野结衣av一区二区av| 巨乳人妻的诱惑在线观看| 波多野结衣高清作品| 99国产精品一区二区蜜桃av| 亚洲人成伊人成综合网2020| 精品久久久久久久久久免费视频| 无限看片的www在线观看| 老司机在亚洲福利影院| 亚洲av成人不卡在线观看播放网| 少妇被粗大的猛进出69影院| 亚洲无线在线观看| 麻豆av在线久日| e午夜精品久久久久久久| 我的亚洲天堂| 成人国产综合亚洲| 日本免费一区二区三区高清不卡| 美女午夜性视频免费| 香蕉丝袜av| 日本一区二区免费在线视频| 观看免费一级毛片| 国产伦在线观看视频一区| 欧美三级亚洲精品| 久久精品影院6| 国产国语露脸激情在线看| 欧美另类亚洲清纯唯美| 人人妻人人看人人澡| 97碰自拍视频| 亚洲男人天堂网一区| 一个人观看的视频www高清免费观看 | 亚洲五月色婷婷综合| 亚洲 国产 在线| 国产免费av片在线观看野外av| 99国产综合亚洲精品| 亚洲人成伊人成综合网2020| 中文在线观看免费www的网站 | 国产单亲对白刺激| 精品日产1卡2卡| 午夜久久久久精精品| 欧美激情久久久久久爽电影| 日韩欧美国产一区二区入口| 天堂√8在线中文| 精品久久久久久久末码| 亚洲一区高清亚洲精品| 免费观看精品视频网站| 老司机在亚洲福利影院| 午夜影院日韩av| 亚洲人成77777在线视频| 亚洲精华国产精华精| 亚洲第一欧美日韩一区二区三区| 成人精品一区二区免费| 男人舔女人的私密视频| 午夜福利免费观看在线| 国产成人影院久久av| 中文字幕人妻熟女乱码| 日韩欧美在线二视频| 韩国av一区二区三区四区| 一级黄色大片毛片| 午夜福利18| av有码第一页| 国产成人精品无人区| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲avbb在线观看| 精品国产乱子伦一区二区三区| 国语自产精品视频在线第100页| 成人av一区二区三区在线看| 久久国产乱子伦精品免费另类| 亚洲av成人av| 亚洲精品久久国产高清桃花| 国产单亲对白刺激| 久久国产亚洲av麻豆专区| 一级a爱视频在线免费观看| 午夜福利18| 黄频高清免费视频| 免费女性裸体啪啪无遮挡网站| 18禁裸乳无遮挡免费网站照片 | 国产单亲对白刺激| 久久青草综合色| 变态另类成人亚洲欧美熟女| 久久久国产欧美日韩av| 亚洲精品美女久久av网站| 69av精品久久久久久| 日本撒尿小便嘘嘘汇集6| 夜夜躁狠狠躁天天躁| 亚洲色图av天堂| 美女大奶头视频| 淫秽高清视频在线观看| 两性夫妻黄色片| 高清在线国产一区| 超碰成人久久| 婷婷精品国产亚洲av在线| 国产av不卡久久| 最近最新免费中文字幕在线| 国产成人欧美在线观看| 成人永久免费在线观看视频| 免费看a级黄色片| a级毛片a级免费在线| aaaaa片日本免费| 怎么达到女性高潮| 国内少妇人妻偷人精品xxx网站 | 亚洲成国产人片在线观看| 免费观看精品视频网站| 成人永久免费在线观看视频| 久久久久久人人人人人| 亚洲精品国产精品久久久不卡| 精品久久久久久久末码| 国产亚洲欧美精品永久| 欧美日韩瑟瑟在线播放| 最好的美女福利视频网| 日韩欧美免费精品| 亚洲第一欧美日韩一区二区三区| 91九色精品人成在线观看| 2021天堂中文幕一二区在线观 | 精品高清国产在线一区| 搞女人的毛片| 精品国产亚洲在线| 他把我摸到了高潮在线观看| 日本 欧美在线| 亚洲avbb在线观看| 色老头精品视频在线观看| 日韩 欧美 亚洲 中文字幕| 99国产综合亚洲精品| 久久草成人影院| 欧美国产精品va在线观看不卡| 人人妻,人人澡人人爽秒播| 女生性感内裤真人,穿戴方法视频| 男人的好看免费观看在线视频 | 色尼玛亚洲综合影院| 国产精品99久久99久久久不卡| 国产免费av片在线观看野外av| 夜夜爽天天搞| 中文在线观看免费www的网站 | 在线观看舔阴道视频| 亚洲成人久久性| 天堂动漫精品| 99久久精品国产亚洲精品| 91在线观看av| 亚洲成人久久性| 嫩草影院精品99| 夜夜夜夜夜久久久久| 久久久久国产精品人妻aⅴ院| 男女床上黄色一级片免费看| 久久性视频一级片| 韩国av一区二区三区四区| 亚洲片人在线观看| 国产精品美女特级片免费视频播放器 | 手机成人av网站| 国产成人系列免费观看| 黄色视频,在线免费观看| 首页视频小说图片口味搜索| 亚洲精品av麻豆狂野| 999久久久国产精品视频| 丝袜美腿诱惑在线| 久久精品国产99精品国产亚洲性色| 欧美国产精品va在线观看不卡| 韩国精品一区二区三区| 国产亚洲精品久久久久5区| 国产亚洲欧美精品永久| 欧美精品亚洲一区二区| 国产精品,欧美在线| 91国产中文字幕| 男男h啪啪无遮挡| 嫩草影视91久久| av在线天堂中文字幕| 中文字幕精品亚洲无线码一区 | 香蕉久久夜色| 波多野结衣av一区二区av| 精品免费久久久久久久清纯| 中出人妻视频一区二区| 成人国产综合亚洲| 他把我摸到了高潮在线观看| 国产极品粉嫩免费观看在线| 99久久99久久久精品蜜桃| 日韩精品青青久久久久久| 午夜亚洲福利在线播放| 正在播放国产对白刺激| 男人舔女人的私密视频| 啪啪无遮挡十八禁网站| 免费在线观看日本一区| 国产精品久久久av美女十八| av在线天堂中文字幕| 最新美女视频免费是黄的| 男女午夜视频在线观看| 国产熟女午夜一区二区三区| 女性被躁到高潮视频| 91成年电影在线观看| 看免费av毛片| avwww免费| 欧美成狂野欧美在线观看| 欧美乱色亚洲激情| 国产一区二区激情短视频| 制服人妻中文乱码| av中文乱码字幕在线| 制服丝袜大香蕉在线| 久久久国产成人精品二区| 白带黄色成豆腐渣| 午夜福利欧美成人| 男人的好看免费观看在线视频 | 91麻豆av在线| 一二三四在线观看免费中文在| 亚洲男人天堂网一区| 成人国语在线视频| 亚洲专区国产一区二区| 亚洲精品粉嫩美女一区| 免费看日本二区| 亚洲国产高清在线一区二区三 | 两性夫妻黄色片| 在线国产一区二区在线| 他把我摸到了高潮在线观看| 99热6这里只有精品| 亚洲熟女毛片儿| 久久久久久久精品吃奶| 精品高清国产在线一区| 少妇被粗大的猛进出69影院| 禁无遮挡网站| 啦啦啦韩国在线观看视频| 色在线成人网| 国产成人影院久久av| 一级毛片高清免费大全| 在线观看一区二区三区| 香蕉国产在线看| 国产熟女xx| 免费av毛片视频| 美女 人体艺术 gogo| 成人国产综合亚洲| 亚洲国产欧洲综合997久久, | 国产亚洲精品久久久久久毛片| 变态另类丝袜制服| 老司机午夜福利在线观看视频| 亚洲国产高清在线一区二区三 | 一级作爱视频免费观看| av电影中文网址| 1024香蕉在线观看| xxxwww97欧美| 熟女电影av网| 日韩视频一区二区在线观看| 在线播放国产精品三级| 韩国精品一区二区三区| 啪啪无遮挡十八禁网站| 丁香六月欧美| 最好的美女福利视频网| 一级毛片女人18水好多| 欧美绝顶高潮抽搐喷水| 岛国视频午夜一区免费看| 婷婷丁香在线五月| ponron亚洲| 桃红色精品国产亚洲av| 国产精品乱码一区二三区的特点| 色综合站精品国产| 观看免费一级毛片| 在线观看一区二区三区| x7x7x7水蜜桃| 亚洲精品一卡2卡三卡4卡5卡| av在线天堂中文字幕| 可以在线观看毛片的网站| 男女做爰动态图高潮gif福利片| xxxwww97欧美| 12—13女人毛片做爰片一| 亚洲aⅴ乱码一区二区在线播放 | 久久国产精品人妻蜜桃| 日韩 欧美 亚洲 中文字幕| bbb黄色大片| 精品不卡国产一区二区三区| 久久中文字幕一级| 欧美久久黑人一区二区| 亚洲成人精品中文字幕电影| 最近在线观看免费完整版| 国产精品1区2区在线观看.| 嫩草影院精品99| 男女那种视频在线观看| 久久国产精品男人的天堂亚洲| 欧美日韩一级在线毛片| 亚洲人成网站在线播放欧美日韩| 听说在线观看完整版免费高清| 淫妇啪啪啪对白视频| 欧美黑人精品巨大| 首页视频小说图片口味搜索| 久久婷婷成人综合色麻豆| 一本久久中文字幕| 777久久人妻少妇嫩草av网站| 91成年电影在线观看| 最近最新中文字幕大全电影3 | 国产精品爽爽va在线观看网站 | 精品不卡国产一区二区三区| 免费女性裸体啪啪无遮挡网站| 亚洲av中文字字幕乱码综合 | 日韩欧美国产在线观看| 国产又色又爽无遮挡免费看| 亚洲人成伊人成综合网2020| 久久久久久免费高清国产稀缺| 搡老妇女老女人老熟妇| 真人做人爱边吃奶动态| 亚洲av美国av| 久久伊人香网站| 久久人人精品亚洲av| 国产精华一区二区三区| 成年版毛片免费区| 欧美激情高清一区二区三区| 男男h啪啪无遮挡| 欧美黑人欧美精品刺激| 午夜久久久在线观看| 日韩成人在线观看一区二区三区| 欧美日韩一级在线毛片| 成人特级黄色片久久久久久久| 日本a在线网址| netflix在线观看网站| 欧美激情极品国产一区二区三区| 91大片在线观看| 少妇被粗大的猛进出69影院| 美女扒开内裤让男人捅视频| 久久国产精品男人的天堂亚洲| 久久久久国产一级毛片高清牌| 亚洲第一av免费看| 亚洲一区高清亚洲精品| 人人妻,人人澡人人爽秒播| 99在线视频只有这里精品首页| 国产亚洲精品久久久久久毛片| 一个人观看的视频www高清免费观看 | 十分钟在线观看高清视频www| 亚洲精品中文字幕在线视频| 亚洲精华国产精华精| 欧美激情久久久久久爽电影| 亚洲人成77777在线视频| 精品国内亚洲2022精品成人| 老汉色av国产亚洲站长工具| 黑人欧美特级aaaaaa片| 国产精品久久久久久精品电影 | 午夜福利18| 亚洲五月婷婷丁香| 岛国在线观看网站| 狠狠狠狠99中文字幕| 色播在线永久视频| 琪琪午夜伦伦电影理论片6080| 看片在线看免费视频| 久久精品影院6| 人妻丰满熟妇av一区二区三区| 国产av一区二区精品久久| 国产成人啪精品午夜网站| 国产又爽黄色视频| 国产单亲对白刺激| 美女大奶头视频| 国产97色在线日韩免费| 波多野结衣高清无吗| 日本撒尿小便嘘嘘汇集6| 精品久久久久久,| 欧美色欧美亚洲另类二区| 成人国产一区最新在线观看| 精品欧美国产一区二区三| 在线观看www视频免费| tocl精华| 日韩免费av在线播放| 人人妻人人澡人人看| 久久精品人妻少妇| 精品一区二区三区四区五区乱码| 999久久久国产精品视频| 可以在线观看的亚洲视频| 18禁观看日本| 国产熟女午夜一区二区三区| 欧美亚洲日本最大视频资源| 午夜福利视频1000在线观看| 一二三四在线观看免费中文在| 最近最新免费中文字幕在线| 国产亚洲av高清不卡| 亚洲aⅴ乱码一区二区在线播放 | 国产精品电影一区二区三区| 国产黄a三级三级三级人| 免费观看精品视频网站| 哪里可以看免费的av片| 国产激情欧美一区二区| 日韩三级视频一区二区三区| 99精品在免费线老司机午夜| 黄片小视频在线播放| 天堂影院成人在线观看| 久久天堂一区二区三区四区| 午夜成年电影在线免费观看| 黄色女人牲交| 成人18禁高潮啪啪吃奶动态图| 亚洲av成人av| 人人妻人人看人人澡| 男女之事视频高清在线观看| 国产精品1区2区在线观看.| 给我免费播放毛片高清在线观看| ponron亚洲| 最新在线观看一区二区三区| 老熟妇仑乱视频hdxx| 1024香蕉在线观看| 日本三级黄在线观看| 国产真人三级小视频在线观看| 欧美日韩亚洲国产一区二区在线观看| 免费av毛片视频| 欧洲精品卡2卡3卡4卡5卡区| 婷婷六月久久综合丁香| 制服人妻中文乱码| 精品久久久久久久毛片微露脸| 好看av亚洲va欧美ⅴa在| 日本 欧美在线| ponron亚洲|