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

    近岸海域典型水體中鎘的水生態(tài)風險評估

    2015-06-05 09:51:45張鈴松王業(yè)耀劉方金小偉孟凡生馬廣文陳平李俊龍丁頁李曌劉喜惠
    生態(tài)毒理學報 2015年2期
    關(guān)鍵詞:遼東灣中鎘海域

    張鈴松,王業(yè)耀,,*,劉方,金小偉,孟凡生,馬廣文,陳平,李俊龍,丁頁,李曌,劉喜惠

    1. 中國環(huán)境科學研究院 環(huán)境基準與風險評估國家重點實驗室,北京 100012 2. 中國環(huán)境監(jiān)測總站,北京 100012

    近岸海域典型水體中鎘的水生態(tài)風險評估

    張鈴松1,王業(yè)耀1,2,*,劉方2,金小偉2,孟凡生1,馬廣文2,陳平2,李俊龍2,丁頁2,李曌2,劉喜惠2

    1. 中國環(huán)境科學研究院 環(huán)境基準與風險評估國家重點實驗室,北京 100012 2. 中國環(huán)境監(jiān)測總站,北京 100012

    為了解我國近岸海域重金屬污染現(xiàn)狀,以鎘為目標污染物,選擇北部灣、渤海灣、杭州灣、遼東灣、閩江口、珠江口、長江口、黃河口等8個典型海灣、河口區(qū)域作為評估目標區(qū)域進行水生態(tài)風險評估。研究結(jié)果顯示:北部灣、渤海灣、杭州灣、閩江口、長江口、珠江口、黃河口等7個海域水體中鎘的暴露水平較低,從2005至2013年其急、慢性暴露的最高潛在影響比例均低于5%,且整體呈下降趨勢。遼東灣水體中鎘從2005至2013年對水生生物的最高潛在影響比例為28.67%(95%置信區(qū)間26.29%~30.06%),5%、10%、15%、20%超過率對應的潛在影響比例分別為18.95%(17.84%~20.23%)、12.51%(11.65%~13.61%)、8.26%(7.40%~9.15%)、5.45%(4.62%~6.16%),其在2005—2010年間呈現(xiàn)略微上升趨勢,2010—2013年間呈現(xiàn)下降趨勢。水體中鎘暴露水平主要受入海河流污染物輸入的影響,其中五里河、大凌河影響最強。

    鎘;近岸海域;水生態(tài)風險;物種敏感度分布

    近岸海域水體承載了豐富的海洋資源,然而近年來,日益加劇的人類活動使沿岸和近海的自然環(huán)境、海洋生物多樣性資源受到一定程度的影響。近岸海域在空間上離陸地更近,受人類活動影響更加強烈,其中污染物的輸入是海洋生態(tài)系統(tǒng)變化的原因之一[1-3],同時也成為人類對海洋資源利用的主要瓶頸之一[4]。

    鎘是人類活動利用、排泄的化學元素之一,其對環(huán)境中的生物具有一定毒性效應。國內(nèi)外研究學者對鎘開展了大量研究工作,結(jié)果表明,鎘能夠造成水生生物生長發(fā)育延緩、畸形率增加、繁殖力下降等,并可引起生物體結(jié)構(gòu)紊亂,造成機體功能障礙,從而使生物物種和群落發(fā)生改變,影響生物多樣性,降低生物資源的利用價值[5-6]。環(huán)境介質(zhì)中鎘含量水平的調(diào)查研究在全世界范圍內(nèi)已廣泛開展,并發(fā)現(xiàn)部分區(qū)域的環(huán)境介質(zhì)中鎘含量相對較高[8-15]。很多國家已將其列為優(yōu)先控制污染物[16],我國環(huán)境保護部門也將鎘列為海洋水質(zhì)必測項目[17-18]。杜建國等[19]基于《中國近海海洋環(huán)境質(zhì)量現(xiàn)狀與背景值研究》[20]中于2006—2007年進行的4次采樣調(diào)查數(shù)據(jù)進行鎘等常見重金屬的水生態(tài)風險評估,結(jié)果發(fā)現(xiàn)黃海北部、大連灣、遼東灣、珠江口受鎘污染較重。但總的來說,當前的文獻報道多集中在環(huán)境介質(zhì)中含量水平的調(diào)查研究方面,而關(guān)于其水生態(tài)風險的報道則很少。Gao等[21]對我國渤海水體、沉積物等環(huán)境介質(zhì)中鎘等重金屬的污染狀況進行綜述,同樣發(fā)現(xiàn)雖然當前已較為全面的開展近岸海域水體中重金屬的監(jiān)測調(diào)查并積累了一定的基礎(chǔ)數(shù)據(jù),然而對于其生態(tài)系統(tǒng)效應仍然不甚了解,進一步的探索研究亟需開展。

    生態(tài)風險評估(ecological risk assessment, ERA)是預測環(huán)境污染物對生態(tài)系統(tǒng)產(chǎn)生有害影響、量化有毒污染物生態(tài)危害的重要方法手段,其最終目的是得出一個風險值,為環(huán)境決策或與其相關(guān)的標準或基準的制定提供參考依據(jù)[22],我國正在大力加強水生態(tài)風險評估研究。海洋水體污染物的水生態(tài)風險評估是環(huán)境管理中的重要組成部分,努力獲得新的評估資料數(shù)據(jù)和結(jié)論可以識別當前的風險水平,為海洋環(huán)境資源的保護和開發(fā)提供有力的支撐,是有關(guān)保護工作的迫切需要,具有較高的實踐意義。因此,本文以鎘為目標污染物,選取北部灣、珠江口、閩江口、長江口、杭州灣、渤海灣、遼東灣8個典型海灣、河口區(qū)域進行水生態(tài)風險評估,以期為我國海洋生態(tài)系統(tǒng)保護提供支持。

    1 研究方法(Methods)

    風險評估一般包括效應評估、暴露評估、風險表征3個步驟[23-25]。物種敏感度分布(species sensitivity distribution, SSD)風險評估法是利用單物種測試中生物受體的毒性數(shù)據(jù)構(gòu)建SSD曲線作為真實水生態(tài)系統(tǒng)物種對污染物敏感性的估計進行效應評估,同時運用環(huán)境實測濃度或模型預測濃度構(gòu)建暴露濃度曲線(exposure concentration distribution, ECD),用于暴露評估;最終,將SSD曲線與ECD曲線進行耦合,進行風險表征[26]。

    1.1 效應評估

    為構(gòu)建SSD曲線,在ECOTOX數(shù)據(jù)庫中檢索鎘對海水水生生物的毒性數(shù)據(jù)。為保證生物受體毒性數(shù)據(jù)的準確性、可靠性、可比性,需要在篩選過程中對毒性實驗的質(zhì)量控制與質(zhì)量保證措施做出嚴格要求。參照文獻中毒性實驗標準方法和數(shù)據(jù)篩選原則[27-30],最終確定毒性實驗應當滿足:

    (1)必須設(shè)置對照實驗,對照組和實驗組的實驗條件一致,對照組中受試生物未出現(xiàn)脅迫和疾病癥狀;

    (2)實驗過程中目標物質(zhì)濃度應以實測為準,其濃度盡量保持恒定;

    (3)實驗過程需控制溫度、pH、DO等水質(zhì)理化參數(shù)恒定;急性毒性的實驗終點為半效應濃度(50% effect concentration, EC50)或半致死濃度(50% lethal concentration, LC50),其測試時間為96 h;慢性毒性的實驗終點為最大無觀察效應濃度(maximum no observed effect concentration, NOEC)或最低觀察效應濃度(lowest observed effect concentration, LOEC),其中魚類測試時間一般應大于28 d,溞類測試時間一般應大于7 d。

    (4)經(jīng)篩選后的受體生物與《中國海洋物種名錄》[31]進行對比,進一步篩選出本地種,并用于構(gòu)建SSD曲線。

    在曲線擬合時,使用目前最為常用的S形logistic模型,如公式(1)所示:

    (1)

    式中y為累積概率;x為生物受體毒性數(shù)據(jù);a、k、xc為模型待估計的參數(shù)。

    1.2 研究區(qū)域及暴露評估

    北部灣、杭州灣、渤海灣、遼東灣、珠江口、閩江口、長江口、黃河口等是自南向北分布于我國近岸海域非常重要的8個海灣、河口,區(qū)域內(nèi)資源豐富,魚、蝦、蟹及貝類等多在此區(qū)域內(nèi)生長繁殖,是我國優(yōu)良的漁場,并且大部分區(qū)域內(nèi)還分布有大量濕地、紅樹林,是亞太地區(qū)候鳥遷徙的重要中轉(zhuǎn)站,具有重要的生態(tài)功能。隨著經(jīng)濟發(fā)展,海域?qū)年懙胤秶鷥?nèi)逐漸形成以“經(jīng)濟開發(fā)區(qū)”等為概念的城市群,承擔區(qū)域的工農(nóng)業(yè)生產(chǎn)、航運等功能,同時該近岸海域水體也成為區(qū)域內(nèi)社會活動所產(chǎn)生的部分污染物的“匯”。以上8個海灣、河口具有典型的代表性,可以作為研究我國近岸海域水生態(tài)風險的模式區(qū)域。

    本文中研究區(qū)共設(shè)置77個監(jiān)測點位,具體坐標設(shè)置如圖1所示,樣品采集及分析測試按照《近岸海域環(huán)境監(jiān)測規(guī)范》[10-11]采集表層水樣,各監(jiān)測點位的監(jiān)測頻次為每年2~3次,采樣在每年4-11月份進行,監(jiān)測數(shù)據(jù)時間跨度為2005—2013年。在暴露評估中將目標區(qū)域內(nèi)全部監(jiān)測點位在2005—2013年內(nèi)所有監(jiān)測數(shù)據(jù)作為樣本,利用S-logistic模型對監(jiān)測數(shù)據(jù)進行擬合構(gòu)建ECD曲線,以此完成暴露評估。

    1.3 風險表征

    概率風險表征通過定性和定量的比較毒性物質(zhì)在環(huán)境中的暴露濃度和生物的劑量-效應濃度的概率分布來表征風險[32],可以更好地表征環(huán)境中毒性物質(zhì)的暴露濃度超過基準值的概率以及可能造成的危害大小而得到廣泛的應用[33-35]。因此,本文采用概率法對生態(tài)風險進行表征。其可以分為安全閾值法(The margin of safety, MOS10)和聯(lián)合概率曲線法(Joint probabilistic curve, JPC)。MOS10法是將物種敏感度分布曲線中累積概率為10%處的濃度值與ECD曲線中累積概率為90%處的濃度值進行比較,以此來量化暴露分布和毒性分布的重疊程度[36]。其具體計算步驟為:

    分別繪制SSD曲線(式(2))與ECD曲線(式(3)),

    y=S(x)

    (2)

    y=E(x)

    (3)

    MOS10值越小,風險越大;MOS10值越大,則表明2條曲線無重疊或重疊很小,風險也較小。

    聯(lián)合概率分布曲線在安全閾值法的基礎(chǔ)上,以毒性數(shù)據(jù)的累積函數(shù)和污染物暴露濃度的反累積函數(shù)繪圖而得。該曲線反映了各潛在影響比例(potential affected fraction, PAF)下暴露濃度超過臨界濃度值的概率(Exceedance),體現(xiàn)了暴露水平與暴露風險之間的關(guān)系[37-38]。其具體步驟為:

    推導ECD曲線的反函數(shù)方程(式(4))與其反累積函數(shù)(式(5))。

    y=E-1(x)

    (4)

    y=1-E(x)

    (5)

    以ECD曲線的反累積函數(shù)作為SSD曲線的自變量對其進行求解(式6)

    y=S[E-1(x)]

    以式(5)、式(6)計算結(jié)果為兩個變量,繪制繪點位置圖,并采用概率分布模型對其進行擬合,則聯(lián)合概率分布曲線可以用式(7)進行表示。

    (7)

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

    2.1 SSD曲線構(gòu)建

    經(jīng)篩選后,獲得的鎘毒性數(shù)據(jù)如表1所示。物種急性毒性數(shù)據(jù)涉及4門5綱8目10科;物種慢性毒性數(shù)據(jù)僅1種。

    由表1可知,物種急性毒性數(shù)據(jù)為4門10科,包括魚類、浮游類、底棲類物種,按照現(xiàn)有水生生物SSD曲線構(gòu)建方法中的要求[51],能夠滿足最小毒性數(shù)據(jù)需求。在所列毒性數(shù)據(jù)中,真寬水蚤(Eurytemora affinis)對水體中鎘的急性毒性最為敏感,其SMAV為147.7 μg·L-1,其次為小頭蟲(Capitella capitata)。而三突柱水虱(Limnoria tripunctata)、鯔魚(Mugil cephalus)、日本對蝦(Penaeus japonicus)、寬葉沙蠶(Nereis grubei)對鎘毒性相對不敏感,其與最敏感物種的SMAV值相差1個數(shù)量級。有研究人員推斷同一分類群內(nèi)物種可能有更加相似的毒性敏感性[52],然而從目前的急性毒性數(shù)據(jù)來看,十足目內(nèi)3個科物種的敏感性相差較大,尚未有類似的規(guī)律,這可能與物種數(shù)量較少有關(guān),此外還可能與不同毒性實驗的實驗條件有一定差異,且缺乏有效的毒性數(shù)據(jù)標準化方法導致的。

    圖1 各海區(qū)監(jiān)測點位設(shè)置圖Fig. 1 Profile of monitoring sites in different sea area

    表1 物種急性毒性數(shù)據(jù)Table 1 The acute toxicity data of saltwater species

    注:SMAV表示種平均急性毒性值。

    Note: SMAV represents species mean acute value.

    表2 物種慢性毒性數(shù)據(jù)Table 2 The chronic toxicity data of saltwater species

    注:SMCV表示種平均慢性毒性值。

    Note: SMCV represents species mean chronic value.

    經(jīng)檢索后,物種慢性毒性數(shù)據(jù)很少,無法利用慢性毒性數(shù)據(jù)直接構(gòu)建SSD曲線,因此,借鑒文獻中急慢性毒性比率法[51]進行外推估計,并以此作為量化慢性毒性效應的標準。經(jīng)模型擬合后分別獲得急性、慢性毒性SSD曲線,如圖2所示。

    2.2 暴露評估

    對各個監(jiān)測站點的水質(zhì)監(jiān)測數(shù)據(jù)采用Kolmogorov-Smirnov、Anderson-Darling檢驗法對其進行正態(tài)分布檢驗,結(jié)果發(fā)現(xiàn)在0.05的置信水平上服從正態(tài)分布或?qū)?shù)正態(tài)分布。按照2.2節(jié)中方法分別對各個海灣、河口水體中鎘的毒性暴露數(shù)據(jù)進行擬合,獲得ECD曲線,其擬合如圖3所示。

    由圖3中8個海區(qū)水體中鎘ECD曲線擬合結(jié)果的調(diào)整可決系數(shù)(AdR2)可知,使用S-logistic模型可以獲得較好的擬合結(jié)果,可以使用該曲線對水體中鎘的暴露水平進行評估。

    2.3 風險表征

    按照2.3節(jié)中方法,經(jīng)計算獲得不同海灣、河口的急性、慢性毒性MOS10值,如表3所示。

    圖2 海洋物種鎘急、慢性毒性SSD曲線Fig. 2 The acute and chronic SSD curve of marine species for Cd

    表3 各海區(qū)水體中鎘的水生態(tài)風險評估MOS10值Table 3 The MOS10 value of Cd in different sea water

    注:MOS10表示安全閾值。

    Note: MOS10 represents the margin of safety.

    由表3可以看出,全部評價區(qū)域的急性、慢性MOS10值均大于1,其中長江口海域的MOS10值最高,分別為1 707.22和132.14,其次為杭州灣,MOS10分別為1 646.14和127.41,說明該海域范圍內(nèi)的水體中基本沒有風險;其他5個海域的MOS10值處于132.69~555.83與10.27~43.02范圍內(nèi),風險也非常低;而遼東灣海域的MOS10值最低,急、慢性評價的MOS10分別為26.01和2.01,該海域范圍內(nèi)的風險需引起一定的關(guān)注。

    為進一步量化各個典型海區(qū)水體中鎘的水生態(tài)風險,運用2.3節(jié)中方法繪制聯(lián)合概率風險曲線,結(jié)果如圖4所示。

    由圖4可知,除遼東灣外,其余7個海域的急、慢性暴露的最高潛在影響比例均低于5%;而遼東灣水體中鎘的最高潛在影響比例為28.67%(95%置信區(qū)間26.29%~30.06%),高于其他7個海域,其5%、10%、15%、20%超過率對應的潛在影響比例分別為18.95%(95%置信區(qū)間17.84%~20.23%)、12.51%(95%置信區(qū)間11.65%~13.61%)、8.26%(95%置信區(qū)間7.40%~9.15%)、5.45% (95%置信區(qū)間4.62%~6.16%)。與杜建國等[19]基于2006年-2007年的調(diào)查數(shù)據(jù)所進行的評估結(jié)果一致。

    聯(lián)合概率曲線法與安全閾值法表征結(jié)果相比,兩者表征結(jié)果基本一致。其中遼東灣的聯(lián)合概率曲線與坐標軸所構(gòu)成的圖形面積最大,反映出較高的潛在影響比例。同時,由表3中各海域的90%EC值、MOS10值及圖4中各海域的聯(lián)合概率曲線形狀可以發(fā)現(xiàn),較高的90%EC值可能未必一定得出較高的風險表征結(jié)果。由于水體中污染物濃度受到輸入強度變化、受納水體的混合擴散稀釋速率等因素的影響而呈現(xiàn)波動性變化,其中個別水體可能呈現(xiàn)小幅度高頻波動,而另一些水體呈現(xiàn)大幅低頻波動。

    圖3 各海區(qū)水體中鎘的暴露濃度曲線Fig. 3 The exposure concentration distribution curve of cadmium in different sea area

    圖4 典型海域水體鎘的水生態(tài)風險聯(lián)合概率曲線Fig. 4 The joint probability curve of Cd ecosystem risk in different seawater

    由圖3中不同海域水體中鎘的暴露分布曲線的形狀可以發(fā)現(xiàn)其暴露濃度值頻數(shù)分布特征。然而經(jīng)曲線擬合后,卻可能獲得相同的90%EC值,但顯然其水生態(tài)風險有一定差異。安全閾值法僅以暴露曲線和敏感度曲線在某一處的重合程度來量化風險大小可能存在較大的不確定性,尤其是當數(shù)據(jù)的樣本容量較小,存在較大的擬合不確定性。聯(lián)合概率曲線則涵蓋了各潛在影響比例水平下暴露濃度超過相應臨界濃度值的概率,反映了水體中污染物濃度的總體頻數(shù)分布情況,能夠更直觀的暴露和效應兩者之間的關(guān)系,在實際的管理中可以為管理者提供更全面的風險評估信息。

    為進一步分析遼東灣海域水體鎘暴露水平的時空變化規(guī)律,繪制其2005—2013年的時間序列圖,如圖5所示。

    圖5 遼東灣水體鎘的暴露水平變化Fig. 5 The box chart and time series of Cd in Liaodong Bay

    由圖5遼東灣水體鎘的暴露水平變化情況可以看出,2005年遼東灣水體中鎘濃度較低,其最高濃度值為0.82 μg·L-1,低于利用ACR法估計的最敏感種真寬水蚤(Eurytemora affinis)的慢性毒性值,因此,其風險也較小。隨后自2006—2010年其峰值增加,并在2010年出現(xiàn)最高值,而后在2010—2013年呈現(xiàn)一定的波動,總體為下降趨勢。

    通常來講,水體中重金屬的外來源主要有陸源輸入、大氣沉降[21]。陸源一般通過入海河流或直排海污染源進入近海水體,其屬于點源范疇;大氣沉降則被劃分為非點源。前者使受納水體呈現(xiàn)較小尺度的空間分布特征,而大氣沉降則很難在較小的空間尺度上呈現(xiàn)區(qū)域特征。由本研究中遼東灣水域水體中鎘濃度的空間分布(圖6)可以發(fā)現(xiàn),遼東灣的9個點位中LD03、LD05、LD08、LD09的鎘濃度相對較高,而其余點位則較低,呈現(xiàn)明顯的空間差異,因此可以判斷遼東灣水體中鎘濃度水平主要受點源影響。進入遼東灣的河流主要有遼河、雙臺子河、大凌河、小凌河、五里河等。其中LD01位于遼河河口處,LD02位于雙臺子河河口處,LD03、LD04、LD05則分別位于大凌河、小凌河、五里河河口處,LD07、LD08、LD09則位于海灣較中心的區(qū)域。因此可以推斷LD03、LD04、LD05點位水體中鎘濃度較高可能與大凌河、小凌河、五里河的污染物輸入有關(guān)。與其他研究人員的研究結(jié)果相比,周秀艷等[53-54]于2003年對遼東灣河口水體、底泥的調(diào)查研究發(fā)現(xiàn)在入遼東灣河口中五里河、大凌河口沉積物中鎘含量較高,而在水體調(diào)查中同時發(fā)現(xiàn)雙臺子河水體中Cd濃度較高,而本研究中靠近雙臺子河的LD02、LD07點位鎘濃度并不高,說明在2005年之后該海域及河口水體中鎘污染空間特征已發(fā)生改變,可能與區(qū)域內(nèi)的鎘排放控制有關(guān)。Wan 等[55]基于2001—2005對遼東灣的調(diào)查數(shù)據(jù)研究發(fā)現(xiàn)五里河可能是遼東灣海域水體中鎘的主要貢獻河流,與本研究所得出的推斷基本一致。比較LD03、LD04、LD05點位的鎘濃度水平可以發(fā)現(xiàn),大凌河、五里河河口水體中鎘濃度相當,均高于小凌河。由遼東灣水體各監(jiān)測斷面的鎘濃度水平空間分布可以發(fā)現(xiàn)水體中鎘暴露水平主要受入海河流污染物輸入的影響,其中五里河、大凌河影響最強,小凌河其次。

    圖6 遼東灣不同監(jiān)測點位鎘暴露水平Fig. 6 Exposure concentration of cadmium in Liaodong Bay

    通過本研究可以得到以下結(jié)論:8個典型海灣、河口中鎘呈現(xiàn)不同的水生態(tài)風險水平,其中北部灣、渤海灣、杭州灣、閩江口、長江口、珠江口、黃河口等7個海域水體中鎘的暴露水平較低,其急、慢性暴露的最高潛在影響比例均低于5%;遼東灣水體中鎘的潛在影響比例高于其他7個區(qū)域,其最高潛在影響比例為28.67%(95%置信區(qū)間26.29%~30.06%),5%、10%、15%、20%超過率對應的潛在影響比例分別為18.95%(17.84%~20.23%)、12.51%(11.65%~13.61%)、8.26%(7.40%~9.15%)、5.45%(4.62%~6.16%)。從遼東灣水體中鎘的暴露水平的時間變化來看,其2005年相對較低,其最高濃度值為0.82 μg·L-1,隨后自2006—2010年其峰值有所增加,并在2010年出現(xiàn)最高值,而后在2010—2013年出現(xiàn)一定的波動,總體呈下降趨勢;從遼東灣水體中鎘的暴露水平空間分布來看,遼東灣水體中鎘暴露水平主要受入海河流污染物輸入的影響,其中五里河、大凌河影響最強、小凌河其次。與其他研究人員的早期研究相比,近年來,遼東灣水體的鎘污染輸入和空間特征發(fā)生變化,因此需要根據(jù)當前情況采取控制措施。

    [1] Lotze H K, Lenihan H S, Bourque B J, et al. Depletion, degradation, and recovery potential of estuaries and coastal seas [J]. Science, 2006, 312(23): 1806-1809

    [2] Pandolfi J M, Bradbury R H, Sala E, et al. Global trajectories of the long-term decline of coral reef ecosystems [J]. Science, 2003, 301(15): 955-958

    [3] Jackson J B C, Kirby M X, Berger W H, et al. Historical overfishing and the recent collapse of coastal ecosystems [J]. Science, 2001, 27(293): 629-637

    [4] 王長友, 王修林, 李克強, 等. 東海陸擾海域銅、鉛、鋅、鎘重金屬排海通量及海洋環(huán)境容量估算[J]. 海洋學報, 2010, 32(4): 62-76

    Wang C Y, Wang X L, Li K Q, et al. The estimation of copper, lead, zinc and cadmium fluxes into the sea area of the East China Sea interferenced by terrigenous matter and their environmental capacities [J]. Acta Oceanologica Sinica, 2010, 32(4): 62-76 (in Chinese)

    [5] 趙元鳳, 呂景才, 宋曉陽, 等. 鎘污染對鰱魚超氧化物歧化酶和過氧化氫酶活性的影響[J]. 農(nóng)業(yè)生物技術(shù)學報, 2002, 10: 267-271

    Zhao Y F, Lu J C, Song X Y, et al. Effect of cadmium on activities of superoxide dismutase and catalase in Aristchthys nobilis [J]. Journal of Agricultural Biotechnology, 2002, 10: 267-271 (in Chinese)

    [6] Weis J S, Smith G M, Zhou T. Altered predator/prey behavior in polluted environments: Implications for fish conservation [J]. Environmental Biology of Fishes, 1999, 55(1-2): 43-51

    [7] Weis J S, Smith G, Zhou T, et al. Effects of contaminants on behavior: Biochemical mechanisms and ecological consequences [J]. BioScience, 2001, 51: 209-217

    [8] Pempkowiak J, Chiffoleau J F, Staniszewski A. The vertical and horizontal distribution of selected trace metals in the Baltic Sea off Poland [J]. Estuarine, Coastal and Shelf Science, 2000, 51(1): 115-125

    [9] Tang D G, Warnken K W, Santschi P H. Distribution and partitioning of trace metals (Cd, Cu, Ni, Pb, Zn) in Galveston Bay waters [J]. Marine Chemistry, 2002, 78(1): 29-45

    [10] Sweeney A, Saudo-Wilhelmy S A. Dissolved metal contamination in the East River-Long Island sound system: Potential biological effects [J]. Marine Pollution Bulletin, 2004, 48(7-8): 663-670

    [11] Gavriil A M, Angelidis M O. Metal and organic carbon distribution in water column of a shallow enclosed bay at the Aegean Sea Archipelago: Kalloni Bay, Island of Lesvos, Greece [J]. Estuarine, Coastal and Shelf Science, 2005, 64(2-3): 200-210

    [12] Li W, Wang N B, Li Q B, et al. Distribution of dissolved metals in seawater of Jinzhou Bay, China [J]. Environmental Toxicology and Chemistry, 2008, 27(1): 43-48

    [13] Huang L L, Pu X M, Pan J F, et al. Heavy metal pollution status in surface sediments of Swan Lake lagoon and Rongcheng Bay in the northern Yellow Sea [J]. Chemosphere, 2013, 93(9): 1957-1964

    [14] Wang Y W, Liang L N, Shi J B, et al. Study on the contamination of heavy metals and their correlations in mollusks collected from coastal sites along the Chinese Bohai Sea [J]. Environment International, 2005, 31(8): 1103-1113

    [15] Thornton G J P, Walsh P D. Heavy metals in the waters of the Nant-y-Fendrod: Change in pollution levels and dynamics associated with the redevelopment of the Lower Swansea Valley, South Wales, UK [J]. The Science of the Total Environment, 2001, 278(123): 45-55

    [16] USEPA. National Recommended Water Quality Criteria [R]. Washington DC: USEPA, 2009

    [17] 環(huán)境保護部. 近岸海域環(huán)境監(jiān)測規(guī)范HJ 442-1998[S]. 北京: 中國環(huán)境科學出版社, 1998

    [18] 環(huán)境保護部. 近岸海域環(huán)境監(jiān)測規(guī)范HJ 442-2008[S]. 北京: 中國環(huán)境科學出版社, 2008

    [19] 杜建國, 趙佳懿, 陳彬, 等. 應用物種敏感性分布評估中國近海和福建主要海灣水體重金屬生態(tài)風險[J]. 生態(tài)毒理學報, 2013, 8(4): 554-560

    Du J G, Zhao J Y, Chen B, et al. Assessing ecological risks of heavy metals to marine organisms in Chinese offshore and Fujian's main bays by species sensitivity distributions [J]. Asian Journal of Ecotoxicology, 2013, 8(4): 554-560 (in Chinese)

    [20] 暨衛(wèi)東. 中國近海海洋環(huán)境質(zhì)量現(xiàn)狀與背景值研究[M]. 北京: 海洋出版社, 2011: 148-158

    [21] Gao X L, Zhou F X, Chen C T A. Pollution status of the Bohai Sea: An overview of the environmental quality assessment related trace metals [J]. Environment International, 2014, 62: 12-30

    [22] 雷炳莉, 黃圣彪, 王子健. 生態(tài)風險評價理論和方法[J]. 化學進展, 2009, 21(1/3): 350-358

    Lei B L, Huang S B, Wang Z J. Theories and methods of ecological risk assessment [J]. Progress in Chemistry, 2009, 21(1/3): 350-358 (in Chinese)

    [23] USEPA. Guidelines for Ecological Risk Assessment [R]. Washington DC: US Environmental Protection Agency, 1998

    [24] Glenn W Suter II. Ecological risk assessment in the United States Environmental Protection Agency: A historical overview [J]. Integrated Environmental Assessment and Management, 2008, 4(3): 285-289

    [25] 孫洪波, 楊桂山, 蘇偉忠, 等. 生態(tài)風險評價研究進展[J]. 生態(tài)學雜志, 2009, 28(2): 335-341

    Sun H B, Yang G S, Su W Z, et al. Research progress on ecological risk assessment [J]. Chinses Journal of Ecology, 2009, 28(2): 335-341 (in Chinese)

    [26] Posthuma L, Glenn W Suter II, Traas T P. Environmental and Ecological Risk Assessment: Species Sensitivity Distributions in Ecotoxicology [M]. Washington D C: Lewis Publishers, 2002

    [27] USEPA. Guidelines for Deriving Numerical National Water Quality Criteria for the Protection of Aquatic Organisms and Their Uses [R]. Washington DC: Office of Research and Development, 1985

    [28] CCME. A Protocol for the Derivation of Water Quality Guidelines for the Protection of Aquatic Life. Canadian Environmental Quality Guidelines [R]. Ottawa: Canadian Council of Ministers of the Environment, 1999

    [29] ASTM. E729-96 (Reapproved 2007) Standard Guide for Conducting Acute Toxicity Tests on Test Materials with Fishes, Macroinvertebrates, and Amphibians [S]. Pennsylvania: American Society for Testing and Materials, 2007

    [30] ASTM. E1193-97 (Reapproved 2004) Standard Guide for Conducting Daphnia magna Life-Cycle Toxicity Tests [S]. Pennsylvania: American Society for Testing and Materials, 2004

    [31] 劉瑞玉. 中國海洋生物名錄[M]. 北京: 科學出版社, 2008: 233-969

    [32] Solomon K R. Overview of recent developments in ecotoxicological risk assessment [J]. Risk Analysis, 1996, 16(5): 627-633

    [33] ECOFRAM. Ecological Committee on FIFRA Risk Assessment Methods Aquatic Reports [R]. Washington DC: ECOFRAM, 1999

    [34] Solomona K, Giesy J, Jonesb P. Probabilistic risk assessment of agrochemicals in the environment [J]. Crop Protection, 2000, 19(8-10): 649-655

    [35] Solomon K R, Baker D B, Richards P R, et al. Ecological risk assessment of atrazine in North American surface waters [J]. Environmental Toxicology and Chemistry, 1996, 15(1): 31-76

    [36] Chen C S. Ecological risk assessment for aquatic species exposed to contaminants in Keelung River, Taiwan [J]. Chemosphere, 2005, 61: 1142-1158

    [37] Wang B, Yu G, Huang J, et al. Tired aquatic ecological risk assessment of organochlorine pesticides and their mixture in Jiangsu reach of Huaihe River, China [J]. Environmental Monitoring and Assessment, 2009, 157: 29-42

    [38] Zolezzi M, Cattaneo C, Tarazona J V. Probabilistic ecological risk assessment of 1,2,4-tricholorobenzene at a former industrial contaminated site [J]. Environmental Science & Technology, 2005, 39(9): 2920-2926

    [39] Reish D J, Martina J M, Piltzb F M, et al. The effect of heavy metals on laboratory populations of two polychaetes with comparisons to the water quality conditions and standards in Southern California marine waters [J]. Water Research, 1976, 10(4): 299-302

    [40] Reish D J, LeMay J A. Toxicity and bioconcentration of metals and organic compounds by polychaeta [J]. Ophelia, 1991, 5(supplement): 653-660

    [41] Hilmy A M, Shabana M B, Daabees A Y. Bioaccumulation of cadmium: Toxicity in Mugil cephalus [J]. Comparative Biochemistry and Physiology Part C: Comparative Pharmacology, 1985, 81(1): 139-143

    [42] Sullivan B K, Buskey E , Miller D C, et al. Effects of copper and cadmium on growth, swimming and predator avoidance in Eurytemora affinis (Copepoda) [J]. Marine Biology, 1983, 77(3): 299-306

    [43] Hong J S, Reish D J. Acute toxicity of cadmium to eight species of marine amphipod and isopod crustaceans from Southern California [J]. Bulletin of Environmental Contamination and Toxicology, 1987, 39(6): 884-888

    [44] Bambang Y, Charmantier G, Thuet P, et al. Effect of cadmium on survival and osmoregulation of various developmental stages of the shrimp Penaeus japonicus (Crustacea: Decapoda) [J]. Marine Biology, 1994, 123(3): 443-450

    [45] Eisler R. Cadmium poisoning in Fundulus heteroclitus (Pisces: Cyprinodontidae) and other marine organisms [J]. Journal of the Fisheries Research Board of Canada, 1971, 28: 1225-1234

    [46] Burton D T, Fisher D J. Acute toxicity of cadmium, copper, zinc, ammonia, 3,3'-dichlorobenzidine, 2,6-dichloro-4-nitroaniline, methylene chloride, and 2,4,6-trichlorophenol to juvenile grass shrimp and killifish [J]. Bulletin of Environmental Contamination and Toxicology, 1990, 44(5): 776-783

    [47] Khan A T, Weis J S, D’Andrea L. Studies of cadmium tolerance in two populations of grass shrimp, Palaemonetes pugio [J]. Bulletin of Environmental Contamination and Toxicology, 1988, 40(5): 30-34

    [48] Eisler R. Acute toxicities of selected heavy metals to the soft-shell clam, Mya arenaria [J]. Bulletin of Environmental Contamination and Toxicology, 1977, 17(2): 137-145

    [49] Eisler R, Hennekey R. Acute toxicities of Cd+2, Cr+6, Hg+2, Ni+2, and Zn+2to estuarine macrofauna [J]. Archives of Environmental Contamination and Toxicology, 1977, 6(1): 315-323

    [50] Ahsanullah M, Negilski D S, MoNey M C. Toxicity of zinc, cadmium and copper to the shrimp Callianasa ustraliensis. I. Effects of individual metals [J]. Marine Biology, 1981, 64: 299-304

    [51] USEPA. Guidelines for Deriving Numerical National Water Quality Criteria for the Protection of Aquatic Organisms and Their Uses [R]. Washington DC: Office of Research and Development, 1985

    [52] 張鈴松. 氨氮水生生物基準推導與風險評估研究[D]. 北京: 北京師范大學, 2013: 26-36

    Zhang L S. Study on aquatic life criteria and ecological assessment of ammonia [D]. Beijing: Beijing Normal University, 2013: 26-36 (in Chinese)

    [53] 周秀艷, 王恩德, 朱恩靜. 遼東灣河口底泥中重金屬的污染評價[J]. 環(huán)境化學, 2004, 23(3): 321-325

    Zhou X Y, Wang E D, Zhu E J. Evaluation on heavy metal pollution in the sediments at the rivermouths around Liaodong Bay [J]. Environmental Chemistry, 2004, 23(3): 321-325 (in Chinese)

    [54] 周秀艷, 王恩德, 劉秀云, 等. 遼東灣河口底質(zhì)重金屬環(huán)境地球化學[J]. 地球化學, 2004, 33(3): 286-290

    Zhou X Y, Wang E D, Liu X Y, et al. Environmental geochemistry of heavy metal in bottom sediments of the river mouths in Liaodong Bay [J]. Geochimica, 2004, 33(3): 286-290 (in Chinese)

    [55] Wan L, Wang N B, Li Q B, et al. Estival distribution of dissolvedmetal concentrations in Liaodong Bay [J]. Bulletin of Environmental Contamination and Toxicology, 2008, 80(4): 311-314

    Ecological Risk Assessment of Cadmium in Chinese Typical Offshore Marine Environment

    Zhang Lingsong1, Wang Yeyao1,2,*, Liu Fang2, Jin Xiaowei2, Meng Fansheng1, Ma Guangwen2, Chen Ping2, Li Junlong2, Ding Ye2, Li Zhao2, Liu Xihui2

    1. State Key Laboratory for Environmental Criteria and Risk Assessment, Chinese Research Academy of Environmental Sciences, Beijing 100012, China 2. China National Environmental Monitoring Center, Beijing 100012, China

    21 October 2014 accepted 17 December 2014

    Probabilistic assessments method was used to assess the ecological risk of cadmium in 8 Chinese typical offshore marine environment including Beibu Bay, Bohai Bay, Hangzhou Bay, Liaodong Bay, Minjiang Estuary, Pearl Estuary, Huanghe Estuary, Yangtze Estuary. The results showed that there is low risk of cadmium in Beibu Bay, Bohai Bay, Hangzhou Bay, Yangtze Estuary, Minjiang Estuary, Pearl River Estuary, Huanghe Estuary. The highest potential affected fractions (PAF) of acute and chronic toxicity were lower than 5%, and overall potential for adverse effects decreased. The PAF of cadmium in Liaodong Bay is relative high, the highest PAF is 28.67% (with 95% confidence interval of 26.29%-30.06%), and the PAF corresponding to exceedance of 5%, 10%, 15%, 20% are 18.95% (17.84%-20.23%), 12.51% (11.65%-13.61%), 8.26% (7.40%-9.15%), 5.45% (4.62%-6.16%), respectively. The overall trend of cadmium concentration in Liaodong Bay is increasing from 2005-2010 and then declining from 2010-2013. It also indicated that the Wuli River and Daling River were main pollution source in the Liaodong Bay.

    cadmium; offshore area; aquatic ecological risk; species sensitivity distribution

    國家水體污染控制與治理科技重大專項(2014ZX07502002);國家自然科學基金(21307165);環(huán)保公益項目(201309008)

    張鈴松(1982-),男,博士,助理研究員,研究方向為水環(huán)境風險評估,E-mail: zlingsong@163.com;

    *通訊作者(Corresponding author), E-mail: yeyaowang@163.com

    10.7524/AJE.1673-5897.20141021001

    2014-10-21 錄用日期:2014-12-17

    1673-5897(2015)2-382-12

    X171.5

    A

    王業(yè)耀(1965-),男,環(huán)境工程博士,研究員,主要研究方向為水污染控制與水生態(tài)風險評估。

    張鈴松, 王業(yè)耀, 劉方, 等. 近岸海域典型水體中鎘的水生態(tài)風險評估[J]. 生態(tài)毒理學報, 2015, 10(2): 382-393

    Zhang L S, Wang Y Y, Liu F, et al. Ecological risk assessment of cadmium in Chinese typical offshore marine environment [J]. Asian Journal of Ecotoxicology, 2015, 10(2): 382-393 (in Chinese)

    猜你喜歡
    遼東灣中鎘海域
    遺落海域
    中國寶玉石(2021年5期)2021-11-18 07:42:32
    埕島海域海上獨立樁拆除方案探討
    四川古藺磺廠工礦廢棄地復墾區(qū)土壤中鎘的賦存形式及生態(tài)風險
    廣東省海域使用統(tǒng)計分析
    酰腙功能化MCM-41 吸附尾礦廢水中鎘離子研究
    水分條件對濕地植物紅蛋吸收鎘及其根際土壤中鎘化學形態(tài)變化的影響
    天然氣分輸站規(guī)劃選址可行性論證初探——以盤錦遼東灣新區(qū)榮興分輸站為例
    遼東灣西部海域表層沉積物重金屬的含量分布與污染評價
    遼東灣坳陷南部旅大10-2構(gòu)造區(qū)源蓋動態(tài)耦合關(guān)系定量評價
    海域評估技術(shù)指引
    美女xxoo啪啪120秒动态图 | 男人的好看免费观看在线视频| 亚洲人与动物交配视频| 男女那种视频在线观看| 黄色日韩在线| www日本黄色视频网| 女生性感内裤真人,穿戴方法视频| 中亚洲国语对白在线视频| 欧美3d第一页| 女生性感内裤真人,穿戴方法视频| 久久久精品欧美日韩精品| 毛片一级片免费看久久久久 | 成人鲁丝片一二三区免费| 窝窝影院91人妻| 成年女人看的毛片在线观看| 精品午夜福利在线看| 亚洲精品久久国产高清桃花| 校园春色视频在线观看| 国产亚洲精品综合一区在线观看| 成人毛片a级毛片在线播放| 91麻豆精品激情在线观看国产| 亚洲aⅴ乱码一区二区在线播放| 麻豆av噜噜一区二区三区| 午夜日韩欧美国产| av天堂在线播放| 女生性感内裤真人,穿戴方法视频| 亚洲欧美日韩高清专用| 日日夜夜操网爽| 两个人的视频大全免费| 国产一区二区在线av高清观看| 国产精华一区二区三区| 怎么达到女性高潮| 国产精品一区二区性色av| 男女床上黄色一级片免费看| 亚洲欧美日韩高清专用| 赤兔流量卡办理| 国产一区二区三区视频了| 黄色女人牲交| 亚洲精品久久国产高清桃花| 怎么达到女性高潮| 一进一出抽搐gif免费好疼| 欧美黑人欧美精品刺激| 国产欧美日韩一区二区精品| 日日摸夜夜添夜夜添小说| 亚洲成人久久性| 一二三四社区在线视频社区8| 亚洲欧美清纯卡通| 淫妇啪啪啪对白视频| 精品不卡国产一区二区三区| 久久久久久久久久黄片| 亚洲av免费高清在线观看| 国产在线精品亚洲第一网站| 国产高潮美女av| 少妇的逼好多水| 欧美色视频一区免费| 欧美一级a爱片免费观看看| 日韩中文字幕欧美一区二区| 一二三四社区在线视频社区8| 久久午夜亚洲精品久久| 亚洲av中文字字幕乱码综合| 免费看a级黄色片| 18禁在线播放成人免费| 中出人妻视频一区二区| 精品福利观看| 97超视频在线观看视频| 成年人黄色毛片网站| 亚洲在线观看片| 久久香蕉精品热| 香蕉av资源在线| 一边摸一边抽搐一进一小说| 热99re8久久精品国产| 可以在线观看毛片的网站| 99久久久亚洲精品蜜臀av| 亚洲熟妇熟女久久| 亚洲成av人片免费观看| 天堂av国产一区二区熟女人妻| 免费在线观看影片大全网站| 一个人免费在线观看的高清视频| 成年免费大片在线观看| 欧美日韩瑟瑟在线播放| 中文字幕人成人乱码亚洲影| 97热精品久久久久久| 日韩欧美国产在线观看| 12—13女人毛片做爰片一| 国产一区二区激情短视频| 国产真实乱freesex| 国产av麻豆久久久久久久| 91字幕亚洲| 90打野战视频偷拍视频| 亚洲av不卡在线观看| 黄色配什么色好看| 久久久久久久久久黄片| 热99re8久久精品国产| 中文资源天堂在线| 色噜噜av男人的天堂激情| 久久亚洲真实| 国产精品电影一区二区三区| 国产精品98久久久久久宅男小说| 久久久久久久久久黄片| 97超级碰碰碰精品色视频在线观看| 日日摸夜夜添夜夜添av毛片 | 国产亚洲欧美在线一区二区| 啪啪无遮挡十八禁网站| 国产美女午夜福利| 动漫黄色视频在线观看| 亚洲 国产 在线| 精品人妻视频免费看| 欧美黄色片欧美黄色片| or卡值多少钱| 丝袜美腿在线中文| 日本撒尿小便嘘嘘汇集6| 天堂影院成人在线观看| 老司机深夜福利视频在线观看| 99久久九九国产精品国产免费| 亚洲人成网站在线播放欧美日韩| 精品人妻1区二区| 亚洲精品色激情综合| 久久6这里有精品| 成人性生交大片免费视频hd| 直男gayav资源| 亚洲成av人片在线播放无| 身体一侧抽搐| 国产精品,欧美在线| 91久久精品国产一区二区成人| 日本撒尿小便嘘嘘汇集6| 国产探花在线观看一区二区| 免费在线观看影片大全网站| 国产欧美日韩精品亚洲av| 成年人黄色毛片网站| 99热这里只有是精品50| 99国产精品一区二区三区| 人人妻人人澡欧美一区二区| 女生性感内裤真人,穿戴方法视频| 日韩有码中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 可以在线观看毛片的网站| 亚洲第一区二区三区不卡| 色吧在线观看| 国产高清有码在线观看视频| 很黄的视频免费| av专区在线播放| av在线蜜桃| 男女下面进入的视频免费午夜| 亚洲无线在线观看| 国产淫片久久久久久久久 | 麻豆av噜噜一区二区三区| 亚洲av成人精品一区久久| 欧美乱色亚洲激情| 日韩欧美国产一区二区入口| 他把我摸到了高潮在线观看| 91在线精品国自产拍蜜月| 欧美日韩福利视频一区二区| 99热这里只有精品一区| 欧美乱色亚洲激情| 久久久久亚洲av毛片大全| 男女床上黄色一级片免费看| 日日干狠狠操夜夜爽| 免费在线观看成人毛片| 成人一区二区视频在线观看| 亚洲精品一卡2卡三卡4卡5卡| 男女下面进入的视频免费午夜| 日日干狠狠操夜夜爽| 欧美极品一区二区三区四区| 观看美女的网站| 99在线视频只有这里精品首页| av在线老鸭窝| 精品人妻偷拍中文字幕| 国产主播在线观看一区二区| 欧美激情在线99| 亚洲天堂国产精品一区在线| 欧美日韩黄片免| 我的老师免费观看完整版| 天堂动漫精品| 精品人妻熟女av久视频| 国产欧美日韩精品一区二区| 每晚都被弄得嗷嗷叫到高潮| 久久精品影院6| 日韩亚洲欧美综合| 久久久国产成人免费| 国产精品久久电影中文字幕| 天堂av国产一区二区熟女人妻| 身体一侧抽搐| 国产蜜桃级精品一区二区三区| 欧美成人性av电影在线观看| 国产精品嫩草影院av在线观看 | 最近中文字幕高清免费大全6 | 久久午夜福利片| 国产精品久久久久久久久免 | 午夜a级毛片| 欧美zozozo另类| 1024手机看黄色片| 可以在线观看毛片的网站| 夜夜看夜夜爽夜夜摸| 国产精品,欧美在线| 国产亚洲av嫩草精品影院| 成熟少妇高潮喷水视频| 国产精品不卡视频一区二区 | 亚洲精品久久国产高清桃花| 国产成人啪精品午夜网站| 亚洲最大成人中文| 欧美在线黄色| 亚洲国产欧洲综合997久久,| 97热精品久久久久久| 欧美另类亚洲清纯唯美| 久久精品人妻少妇| 啪啪无遮挡十八禁网站| 中文资源天堂在线| 亚洲精品色激情综合| 美女高潮的动态| 成年女人看的毛片在线观看| 亚洲aⅴ乱码一区二区在线播放| 看免费av毛片| 国产欧美日韩一区二区三| 中文字幕精品亚洲无线码一区| 成年版毛片免费区| 在线观看美女被高潮喷水网站 | 男人舔奶头视频| 久久午夜福利片| 午夜日韩欧美国产| 美女高潮喷水抽搐中文字幕| 欧美高清性xxxxhd video| 身体一侧抽搐| 男女视频在线观看网站免费| 亚洲午夜理论影院| 黄片小视频在线播放| 久久99热6这里只有精品| 女同久久另类99精品国产91| 精品免费久久久久久久清纯| 日日摸夜夜添夜夜添小说| 热99在线观看视频| 欧美黄色淫秽网站| 国产高清视频在线播放一区| av国产免费在线观看| 亚洲性夜色夜夜综合| 国模一区二区三区四区视频| 亚洲第一电影网av| 在线国产一区二区在线| 国产乱人视频| 黄片小视频在线播放| 简卡轻食公司| 国产成人a区在线观看| h日本视频在线播放| 亚洲综合色惰| 18美女黄网站色大片免费观看| 99久国产av精品| a级一级毛片免费在线观看| 国产高清视频在线播放一区| 国产三级黄色录像| 97超级碰碰碰精品色视频在线观看| 91九色精品人成在线观看| av国产免费在线观看| 女人十人毛片免费观看3o分钟| 日本三级黄在线观看| 麻豆成人av在线观看| 在线观看一区二区三区| 午夜亚洲福利在线播放| 国产精品久久久久久久久免 | 88av欧美| 真人一进一出gif抽搐免费| 麻豆一二三区av精品| 久久人人精品亚洲av| 精品久久久久久久末码| 99热这里只有是精品50| 亚洲 国产 在线| 欧美性感艳星| 国产日本99.免费观看| 看免费av毛片| 人人妻人人澡欧美一区二区| 天美传媒精品一区二区| 又爽又黄无遮挡网站| 国产精品国产高清国产av| 99久久精品热视频| 日韩高清综合在线| 日韩欧美国产在线观看| 一区二区三区高清视频在线| 国产一区二区激情短视频| 两人在一起打扑克的视频| 伊人久久精品亚洲午夜| 五月伊人婷婷丁香| 桃红色精品国产亚洲av| 国产激情偷乱视频一区二区| 午夜福利18| 亚洲黑人精品在线| 又爽又黄a免费视频| 91久久精品电影网| 色综合站精品国产| 少妇人妻精品综合一区二区 | 欧美+日韩+精品| 精品人妻1区二区| 99国产综合亚洲精品| 亚洲人与动物交配视频| 91狼人影院| 无遮挡黄片免费观看| 女人十人毛片免费观看3o分钟| 亚洲男人的天堂狠狠| 一个人看的www免费观看视频| www.色视频.com| 男插女下体视频免费在线播放| 少妇的逼水好多| 91午夜精品亚洲一区二区三区 | 久久久久久国产a免费观看| 国内精品久久久久久久电影| 成人永久免费在线观看视频| 亚洲欧美日韩无卡精品| 精品欧美国产一区二区三| 天天一区二区日本电影三级| 丰满人妻一区二区三区视频av| 女人被狂操c到高潮| 91字幕亚洲| 欧美成人性av电影在线观看| 性欧美人与动物交配| 精品不卡国产一区二区三区| 欧美成人a在线观看| 少妇熟女aⅴ在线视频| 黄色视频,在线免费观看| 日韩大尺度精品在线看网址| 欧美黑人欧美精品刺激| 国产综合懂色| 欧美黄色片欧美黄色片| 亚洲在线自拍视频| 欧美黑人欧美精品刺激| 久久久久久久亚洲中文字幕 | 国产在线男女| a级毛片免费高清观看在线播放| 看免费av毛片| 色哟哟·www| 国产在线精品亚洲第一网站| 亚洲国产欧洲综合997久久,| 成熟少妇高潮喷水视频| 人人妻,人人澡人人爽秒播| 国产真实乱freesex| 日本黄大片高清| 在线观看午夜福利视频| 一级黄色大片毛片| 丝袜美腿在线中文| 国产精品98久久久久久宅男小说| 精品久久久久久久末码| 国产视频一区二区在线看| 性插视频无遮挡在线免费观看| 黄色一级大片看看| 色综合站精品国产| 亚洲成av人片在线播放无| 欧美日韩中文字幕国产精品一区二区三区| a级一级毛片免费在线观看| 床上黄色一级片| 欧美一级a爱片免费观看看| 亚洲av美国av| 欧美最黄视频在线播放免费| 床上黄色一级片| 免费人成在线观看视频色| 久久人人爽人人爽人人片va | 他把我摸到了高潮在线观看| 国产 一区 欧美 日韩| 中文字幕人成人乱码亚洲影| 亚洲国产精品久久男人天堂| 国产免费av片在线观看野外av| 欧美高清性xxxxhd video| 男女下面进入的视频免费午夜| 欧美成人免费av一区二区三区| 美女免费视频网站| 一级av片app| 亚洲真实伦在线观看| 很黄的视频免费| 人人妻,人人澡人人爽秒播| x7x7x7水蜜桃| 婷婷精品国产亚洲av在线| 美女 人体艺术 gogo| 午夜福利高清视频| 久久久久性生活片| 久久国产乱子免费精品| 亚洲男人的天堂狠狠| 国产真实伦视频高清在线观看 | 日韩精品中文字幕看吧| 中亚洲国语对白在线视频| 亚洲专区国产一区二区| 嫁个100分男人电影在线观看| 国产精品久久久久久久电影| 精品一区二区三区视频在线| 亚洲精品亚洲一区二区| xxxwww97欧美| 丁香欧美五月| 精品99又大又爽又粗少妇毛片 | 午夜福利在线观看吧| 国产黄片美女视频| 高清毛片免费观看视频网站| 1000部很黄的大片| 99热6这里只有精品| 亚洲av成人不卡在线观看播放网| 夜夜夜夜夜久久久久| 亚洲自拍偷在线| 欧美丝袜亚洲另类 | 最近视频中文字幕2019在线8| 在线观看舔阴道视频| 亚洲av一区综合| 亚洲最大成人中文| 18美女黄网站色大片免费观看| 国产老妇女一区| 久久久久九九精品影院| 精品一区二区三区人妻视频| 精品人妻熟女av久视频| 午夜影院日韩av| 美女高潮的动态| 亚洲成人中文字幕在线播放| www.熟女人妻精品国产| 亚洲av成人av| 亚洲午夜理论影院| 欧美成人性av电影在线观看| 精品不卡国产一区二区三区| 久9热在线精品视频| 中文字幕av成人在线电影| 久久久久性生活片| 国产一区二区在线观看日韩| 高清在线国产一区| 亚洲18禁久久av| 欧美最黄视频在线播放免费| 亚洲性夜色夜夜综合| 99国产精品一区二区三区| 亚洲国产精品成人综合色| 一级黄色大片毛片| 国产淫片久久久久久久久 | 男人和女人高潮做爰伦理| 日韩免费av在线播放| 久久久久性生活片| 一a级毛片在线观看| 亚洲av电影不卡..在线观看| 国产精品一及| 国产激情偷乱视频一区二区| 日本免费a在线| 免费无遮挡裸体视频| 久久亚洲真实| 亚洲一区二区三区色噜噜| 99精品久久久久人妻精品| 成人国产一区最新在线观看| 99久国产av精品| 国产免费av片在线观看野外av| 国产在线精品亚洲第一网站| 露出奶头的视频| 日本熟妇午夜| 亚洲精华国产精华精| 国产午夜福利久久久久久| 国产淫片久久久久久久久 | 简卡轻食公司| 精品久久久久久,| 欧美另类亚洲清纯唯美| 亚洲自偷自拍三级| 18+在线观看网站| 国产伦人伦偷精品视频| 在线观看66精品国产| 在现免费观看毛片| 久久欧美精品欧美久久欧美| 中出人妻视频一区二区| 美女黄网站色视频| 人人妻人人看人人澡| 少妇裸体淫交视频免费看高清| 日韩精品中文字幕看吧| 欧美乱妇无乱码| 欧美又色又爽又黄视频| 国产精品免费一区二区三区在线| 国产 一区 欧美 日韩| 国产一区二区在线av高清观看| 91九色精品人成在线观看| 免费人成视频x8x8入口观看| 精品国产亚洲在线| 99久久精品国产亚洲精品| 男女下面进入的视频免费午夜| 真实男女啪啪啪动态图| 亚洲狠狠婷婷综合久久图片| 老熟妇乱子伦视频在线观看| 1024手机看黄色片| 国产真实伦视频高清在线观看 | 亚洲成人免费电影在线观看| 中文字幕精品亚洲无线码一区| 欧美在线黄色| 亚洲欧美日韩高清在线视频| 极品教师在线免费播放| 麻豆成人av在线观看| 国产人妻一区二区三区在| 精品人妻视频免费看| 级片在线观看| 国产淫片久久久久久久久 | 久99久视频精品免费| 久久人人精品亚洲av| 直男gayav资源| 国产av不卡久久| 一本综合久久免费| 高清毛片免费观看视频网站| 99国产极品粉嫩在线观看| 蜜桃久久精品国产亚洲av| 亚洲av成人av| 伊人久久精品亚洲午夜| 亚洲av免费高清在线观看| 2021天堂中文幕一二区在线观| 亚洲avbb在线观看| 免费高清视频大片| 亚洲无线观看免费| 俺也久久电影网| 欧美又色又爽又黄视频| 亚洲第一电影网av| 日本精品一区二区三区蜜桃| 午夜精品一区二区三区免费看| 国产高清视频在线观看网站| 国产精品电影一区二区三区| 精品久久国产蜜桃| 国产在线男女| 我要搜黄色片| 精品国产亚洲在线| 首页视频小说图片口味搜索| 亚洲激情在线av| 国产精品人妻久久久久久| 99国产精品一区二区蜜桃av| 99久久久亚洲精品蜜臀av| 亚洲第一欧美日韩一区二区三区| 精品午夜福利在线看| 国产伦一二天堂av在线观看| 18禁裸乳无遮挡免费网站照片| 国产精品久久久久久亚洲av鲁大| 蜜桃久久精品国产亚洲av| 九色成人免费人妻av| 欧美激情国产日韩精品一区| 欧美乱妇无乱码| 日韩成人在线观看一区二区三区| 国产精品一区二区免费欧美| 观看美女的网站| 精品国产亚洲在线| 免费看日本二区| 亚洲自偷自拍三级| 日本熟妇午夜| x7x7x7水蜜桃| 免费av毛片视频| 国产亚洲精品综合一区在线观看| 在线观看午夜福利视频| av在线观看视频网站免费| 中文字幕人妻熟人妻熟丝袜美| 九九热线精品视视频播放| 99久久九九国产精品国产免费| 亚洲人成电影免费在线| avwww免费| 脱女人内裤的视频| 欧美zozozo另类| 午夜福利高清视频| 直男gayav资源| 嫩草影院精品99| 欧美日韩综合久久久久久 | www.熟女人妻精品国产| 国产精品爽爽va在线观看网站| 三级男女做爰猛烈吃奶摸视频| 热99re8久久精品国产| 久久久成人免费电影| 亚洲午夜理论影院| 亚洲五月婷婷丁香| 少妇人妻一区二区三区视频| 看免费av毛片| 欧美黑人巨大hd| 亚洲专区中文字幕在线| 丰满人妻熟妇乱又伦精品不卡| 九九热线精品视视频播放| 精品不卡国产一区二区三区| 美女高潮的动态| www.999成人在线观看| 成人永久免费在线观看视频| 真人一进一出gif抽搐免费| 精品久久久久久久久亚洲 | 欧美又色又爽又黄视频| 最近最新中文字幕大全电影3| 少妇裸体淫交视频免费看高清| 国产一区二区三区视频了| 成熟少妇高潮喷水视频| 亚洲av日韩精品久久久久久密| 免费av不卡在线播放| 99热这里只有是精品在线观看 | 国产伦在线观看视频一区| 在现免费观看毛片| 一区福利在线观看| 午夜免费男女啪啪视频观看 | 亚洲av熟女| 国产爱豆传媒在线观看| 尤物成人国产欧美一区二区三区| 免费大片18禁| 91久久精品电影网| 高清在线国产一区| 在线观看美女被高潮喷水网站 | 男插女下体视频免费在线播放| 亚洲国产精品999在线| 美女被艹到高潮喷水动态| 成人无遮挡网站| 国产午夜精品久久久久久一区二区三区 | 国内精品美女久久久久久| 亚洲av不卡在线观看| 久久久久久大精品| 女同久久另类99精品国产91| 免费看日本二区| 国产亚洲欧美在线一区二区| 国产视频一区二区在线看| 成年女人毛片免费观看观看9| 国产黄a三级三级三级人| 三级男女做爰猛烈吃奶摸视频| 女人被狂操c到高潮| 国产一区二区在线av高清观看| 精品不卡国产一区二区三区| 国产精品伦人一区二区| 久久久久久国产a免费观看| 亚洲成人免费电影在线观看| 色综合欧美亚洲国产小说| 女生性感内裤真人,穿戴方法视频| 欧美区成人在线视频| 婷婷丁香在线五月| 亚洲av日韩精品久久久久久密| 亚洲国产高清在线一区二区三| 午夜福利成人在线免费观看| 一个人免费在线观看的高清视频| 永久网站在线| 日本在线视频免费播放| 国产黄色小视频在线观看| 欧美+亚洲+日韩+国产| 精品久久久久久久久久免费视频|