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

    湛江東海島紅樹林濕地表層土壤重金屬空間分布特征及生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)

    2018-09-03 03:48:56羅松英王嘉琦周敏葉嘉雯柯思茵陳東平賴曉清
    關(guān)鍵詞:紅樹林金屬元素站位

    羅松英,王嘉琦,周敏,葉嘉雯,柯思茵,陳東平,賴曉清

    嶺南師范學(xué)院地理系,廣東 湛江 524048

    紅樹林是分布在海陸交替過渡地帶的濕地生態(tài)系統(tǒng),受陸源包括人類活動(dòng)以及海域環(huán)境信息的疊加影響(李瑞利等,2012;Liu et al.,2014)。由于紅樹林生境潮流速度較低,土壤沉積物具有粘粒含量高、有機(jī)質(zhì)含量高、硫化物含量高、還原性強(qiáng)等特征,能夠大量接受來自潮汐和地表徑流所攜帶的重金屬污染物,往往成為潛在的重金屬富集區(qū)(程皓等,2009;Bosire et al.,2005)。重金屬污染具有隱蔽性、毒性高、不易分解等特點(diǎn),可通過食物鏈的放大作用最終危害人體健康,因此重金屬污染備受關(guān)注(Nath et al.,2013)。對(duì)于整個(gè)紅樹林生態(tài)系統(tǒng)而言,土壤沉積物是重金屬污染物的主要貯存庫(kù),絕大部分的重金屬沉淀固定在土壤沉積物中(戴紀(jì)翠等,2009)。紅樹林土壤重金屬含量的高低與其生境理化性質(zhì)、土壤重金屬背景值高低和受附近城市排放重金屬污染的程度等密切相關(guān)(劉景春,2006;李柳強(qiáng)等,2008)。國(guó)內(nèi)外對(duì)不同海岸紅樹林重金屬污染的研究結(jié)果表明,不同紅樹林區(qū)重金屬污染程度不同,受陸域環(huán)境影響較大的紅樹林濕地重金屬潛在生態(tài)風(fēng)險(xiǎn)有逐漸加強(qiáng)的趨勢(shì)(郭菊蘭等,2013;Sandilyan et al.,2014;Li et al.,2015;Miola et al.,2016)。

    東海島位于廣東省湛江市東南部,地勢(shì)平坦,島上紅樹林濕地面積約1475 hm2,屬于湛江紅樹林國(guó)家級(jí)自然保護(hù)區(qū)范圍(王中丙,2013;朱小丹,2016)。自1992年東海島開發(fā)試驗(yàn)區(qū)設(shè)立以來,島上填海項(xiàng)目逐漸增多,潮間帶區(qū)域逐漸被吞食,紅樹林濕地面積大幅減少(李曉敏,2008)。隨著鋼鐵、化工、石化等項(xiàng)目相繼落戶東海島,海島周邊生態(tài)環(huán)境壓力不容忽視。同時(shí)由于東海島及周邊地區(qū)城市化和工業(yè)化進(jìn)程的急劇加速,工業(yè)廢水、城市生活污水、養(yǎng)殖排污、船舶污染和農(nóng)業(yè)面源污染對(duì)生態(tài)環(huán)境造成持續(xù)性壓力,生長(zhǎng)于潮間帶環(huán)境的紅樹林濕地重金屬污染問題日益凸顯(劉加飛,2013)。目前,針對(duì)東海島紅樹林濕地土壤重金屬的研究也少見報(bào)道,僅見劉加飛(2013)對(duì)東海島潮間帶表層沉積物重金屬污染進(jìn)行過研究,其研究對(duì)象并未專門針對(duì)紅樹林濕地生態(tài)系統(tǒng),且僅對(duì)Cu、Zn、Cd、Pb、Cr這5種重金屬元素進(jìn)行污染評(píng)價(jià),研究不夠系統(tǒng)全面。

    本研究通過分析東海島紅樹林濕地土壤重金屬含量及其空間分布特征,采用潛在生態(tài)危害指數(shù)法和多元統(tǒng)計(jì)法對(duì)紅樹林濕地土壤重金屬污染現(xiàn)狀進(jìn)行評(píng)價(jià)及來源解析,以期為東海島紅樹林環(huán)境污染監(jiān)測(cè)與生態(tài)保護(hù)提供科學(xué)參考。

    1 材料與方法

    1.1 研究區(qū)概況

    湛 江 東 海 島 ( 20°55'~21°05'N , 110°11'~110°33'E)位于湛江灣與雷州灣之間,全島面積約為289.49 km2,海岸線為126 km,是中國(guó)第五大島嶼(王中丙,2013)。屬于北熱帶季風(fēng)海洋性氣候區(qū),年平均氣溫23 ℃,年均降水雨量1800~2000 mm,長(zhǎng)夏無寒冬,熱帶氣旋是影響該地區(qū)的主要災(zāi)害性天氣。潮汐類型為不規(guī)則半日潮,波浪以風(fēng)浪為主,多年平均波高為0.9 m(陳則實(shí),1999)。紅樹林主要分布在東海島南岸、西岸以及北岸靠西段,主要紅樹植物為白骨壤(Avicennia marina),少量的秋茄(Kandelia candel)、無瓣海桑(Sonneratia apetala)和紅海欖(Rhizophora stylosa)。西岸紅樹林常見分布于河口,周圍有大量的水產(chǎn)養(yǎng)殖池,養(yǎng)殖廢水直接排入紅樹林;南岸紅樹林主要分布于內(nèi)灣或背風(fēng)區(qū),面向雷州灣,有少量水產(chǎn)養(yǎng)殖池或鹽場(chǎng)。北岸有少量紅樹林分布,位于東參碼頭和對(duì)岸小島上。自東海島成立開發(fā)實(shí)驗(yàn)區(qū)以來,該地的交通、通訊、電力、海水養(yǎng)殖、鋼鐵、石化、農(nóng)副產(chǎn)品加工、飼料、城郊型農(nóng)業(yè)等臨港工業(yè)和資源型工業(yè)的規(guī)模和效益不斷擴(kuò)大和提高,已經(jīng)成為開發(fā)區(qū)經(jīng)濟(jì)的規(guī)模產(chǎn)業(yè)(劉加飛,2013)。據(jù)調(diào)查資料顯示,湛江灣和雷州灣的主要污染來源有城市生活污水、工業(yè)廢水、養(yǎng)殖排污、船舶污染和農(nóng)業(yè)面源污染(劉芳文等,2015)。

    1.2 樣品采集與研究方法

    1.2.1 樣品采集

    綜合考慮東海島紅樹林分布面積、潮位變化及人類活動(dòng)影響等因素,于2017年6月沿著海岸線在紅樹林典型分布區(qū)布設(shè)了7個(gè)采樣站位,從東南向西北方向分別是YC、SYC、JJC、XWC、TMHHK、DCMT和DTSD。采樣站位分布如圖1所示,每個(gè)站位采集土壤樣品2~3個(gè),共15個(gè)樣品。采樣時(shí)先用塑料小鏟刮去表層枯枝殘葉,然后使用直徑為6 cm的PVC管壓入淤泥中采集表層0~15 cm土壤樣品,每個(gè)樣點(diǎn)采用梅花采樣法取5個(gè)土樣混合,采集約1 kg土樣,所采集的樣品裝入密封的乙烯封口袋中,帶回實(shí)驗(yàn)室。

    1.2.2 樣品處理及測(cè)試方法

    野外采集的土壤樣品經(jīng)過自然風(fēng)干、研磨、過篩(100目)處理后,測(cè)定砷(As)、鎘(Cd)、鉻(Cr)、銅(Cu)汞(Hg)、鎳(Ni)、鉛(Pb)、鋅(Zn)8種重金屬元素的質(zhì)量分?jǐn)?shù)。測(cè)試地點(diǎn)為廣州澳實(shí)礦物實(shí)驗(yàn)室,其測(cè)試過程為:稱取兩份樣品,經(jīng)過王水消解和四酸(高氯酸、硝酸、氫氟酸和鹽酸)消解處理后,分別采用電感耦合等離子體發(fā)射光譜(ICP-AES)和電感耦合等離子體質(zhì)譜儀(ICP-MS)進(jìn)行測(cè)試,綜合取值為最后的檢測(cè)結(jié)果(陳碧珊等,2017)。對(duì)每3個(gè)樣品進(jìn)行平行測(cè)定,同時(shí)測(cè)定空白樣。其中ICP-AES型號(hào)是Varian VISTA,ICP-MS型號(hào)是Agilent 7700x。微量分析的標(biāo)樣準(zhǔn)確度RSD小于10%,回收率高于90%。

    1.2.3 研究方法

    本研究采用潛在生態(tài)危害指數(shù)法對(duì)東海島紅樹林土壤重金屬進(jìn)行污染及潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià),該方法相對(duì)快速簡(jiǎn)便,綜合考慮了重金屬的毒性、遷移轉(zhuǎn)化規(guī)律、區(qū)域敏感性以及重金屬區(qū)域背景值的差異,可以綜合反映沉積物中重金屬對(duì)生態(tài)環(huán)境的影響潛力,并給出危害程度定量劃分方法,是適合異源沉積物之間評(píng)價(jià)的有效手段(鄧?yán)龋?014)。潛在生態(tài)危害指數(shù)法計(jì)算公式如下:

    (1)單因子污染指數(shù)

    式中,Ci為重金屬i的污染系數(shù);Ci為重金屬

    fi i的實(shí)測(cè)濃度含量; Cn為重金屬的評(píng)價(jià)參比值。東海島典型地帶性土壤為磚紅壤性紅壤,故本研究采用廣東省土壤的環(huán)境背景值中磚紅壤中的重金屬含量作為背景值。

    (2)某一區(qū)域重金屬i的潛在生態(tài)危害系數(shù)

    i反映了重金屬的毒性水平及生物對(duì)重金屬污染的敏感程度,其毒性響應(yīng)系數(shù)見表 1(徐爭(zhēng)啟等,2008)。

    表1 重金屬毒性系數(shù)Table 1 The toxicity factor of heavy metals

    (3)多種重金屬潛在生態(tài)風(fēng)險(xiǎn)指數(shù)

    式中,RI(the potential ecological risk index)為多種重金屬潛在生態(tài)風(fēng)險(xiǎn)指數(shù);為金屬i的潛在生態(tài)危害系數(shù);為重金屬毒性響應(yīng)系數(shù),反映重金屬的毒性水平及生物對(duì)金屬污染的敏感程度。重金屬污染潛在生態(tài)風(fēng)險(xiǎn)指標(biāo)和分級(jí)關(guān)系如表2所示(黃建國(guó),2007)。

    1.2.4 數(shù)據(jù)處理

    多元統(tǒng)計(jì)分析通過SPSS 22.0軟件完成,采用相關(guān)性和主成分分析法進(jìn)行重金屬來源分析。相關(guān)性分析是指對(duì)兩個(gè)或多個(gè)具備相關(guān)性的變量元素進(jìn)行分析,從而衡量?jī)蓚€(gè)變量因素的相關(guān)程度。通過重金屬間的相關(guān)系數(shù)可更加準(zhǔn)確地表明其來源途徑的相似程度(雷凌明等,2014)。主成分分析是研究變量之間相互關(guān)系的多變量分析技術(shù),通過對(duì)原變量的線性變換,形成了一系列新的相互獨(dú)立變量,稱之為主成分(孫欽幫等,2017)。在進(jìn)行主成分分析之前,采用KMO和Bartlett法對(duì)原始數(shù)據(jù)集進(jìn)行因子分析適宜性檢驗(yàn)(李玉等,2006)。KMO值為0.700,大于0.5;Bartlett的球形度檢驗(yàn)相伴概率為 0.000,小于顯著性水平 0.05,表明原始數(shù)據(jù)集適合進(jìn)行因子分析(姚曉飛等,2011)。

    表2 重金屬污染潛在生態(tài)風(fēng)險(xiǎn)指標(biāo)和分級(jí)分析Table 2 Indices and grades of potential ecological risk index of heavy metal pollution

    表3 東海島紅樹林濕地表層土壤重金屬元素質(zhì)量分?jǐn)?shù)Table 3 The mass fraction of heavy metals in the surface soil of mangrove wetland in Donghai Island

    2 結(jié)果與討論

    2.1 土壤重金屬質(zhì)量分?jǐn)?shù)特征及空間分布

    由表 3 可知,As、Cd、Cr、Cu、Hg、Ni、Pb、Zn 等 8種元素重金屬質(zhì)量分?jǐn)?shù)平均值超過雷州半島土壤環(huán)境背景值(許煉烽等,1996);8種元素重金屬質(zhì)量分?jǐn)?shù)平均值在《土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB15618—1995)一級(jí)標(biāo)準(zhǔn)范圍內(nèi),但YC、XWC和TMHHK站位采樣區(qū)的個(gè)別元素,如As、Hg、Pb質(zhì)量分?jǐn)?shù)超過《土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)》一級(jí)標(biāo)準(zhǔn)值。從重金屬質(zhì)量分?jǐn)?shù)范圍看(表 3),8種重金屬元素質(zhì)量分?jǐn)?shù)范圍差異較大,As為 6.800~26.600 mg?kg-1,Cd 為 0.023~0.068 mg?kg-1,Cr為 21.000~74.000 mg?kg-1,Cu 為 5.700~21.000 mg?kg-1,Hg 為0.037~0.200 mg?kg-1,Ni為 8.700~ 32.000 mg?kg-1,Pb 為 16.000~45.800 mg?kg-1,Zn 為 32.000~112.000 mg·kg-1。變異系數(shù)可以反映各重金屬元素的離散程度,且變異系數(shù)越大說明土壤重金屬污染受人為活動(dòng)干擾越強(qiáng)烈(趙秀峰等,2010)。從變異系數(shù)上看(表3),東海島紅樹林濕地表層土壤中重金屬元素Pb的變異系數(shù)最小,為33.91%,離散程度較低,反映該元素分布相對(duì)較均勻,推測(cè)目前受人為活動(dòng)干擾不明顯;As、Cd、Cr、Cu、Ni、Zn的變異系數(shù)均在40%以上,分布較不均勻,說明受到一定程度的人為活動(dòng)干擾;Hg的變異系數(shù)最大,為62.16%,分布最為不均,推測(cè)其受人為活動(dòng)干擾也較強(qiáng)烈。

    從重金屬質(zhì)量分?jǐn)?shù)空間分布上看,7個(gè)采樣站位重金屬質(zhì)量分?jǐn)?shù)差異較大,選取土壤中超標(biāo)狀況比較嚴(yán)重的4種重金屬元素As、Hg、Pb和Zn進(jìn)行分析,結(jié)果如圖2所示。其中As質(zhì)量分?jǐn)?shù)平均值在7個(gè)采樣站位中均高于雷州半島土壤環(huán)境背景值,以西部靠近水產(chǎn)養(yǎng)殖地的XWC和TMHHK 2個(gè)采樣站位為最高,超出國(guó)家土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)值范圍;Hg含量平均值在 7個(gè)采樣站位中也均高于雷州半島土壤環(huán)境背景值,以西南岸XWC采樣站位為最高,超出國(guó)家標(biāo)準(zhǔn)值范圍;Pb除了在SYC、JJC和DCMT這3個(gè)采樣站位外,在其他4個(gè)采樣站位的質(zhì)量分?jǐn)?shù)平均值均接近或高于雷州半島土壤環(huán)境背景值,其中YC和TMHHK采樣站位的Pb質(zhì)量分?jǐn)?shù)平均值均超出國(guó)家標(biāo)準(zhǔn)范圍;Zn質(zhì)量分?jǐn)?shù)平均值在7個(gè)采樣站位中均高于雷州半島土壤環(huán)境背景值,以東南岸YC采樣站位為最高,但未超出國(guó)家標(biāo)準(zhǔn)值范圍。

    2.2 表層土壤重金屬潛在生態(tài)風(fēng)險(xiǎn)分析及空間分布

    2.2.1 表層土壤重金屬污染分析及空間分布

    東海島紅樹林濕地表層土壤不同重金屬元素單因子污染評(píng)價(jià)結(jié)果如表 4所示,As、Cd、Cr、Cu、Hg、Ni、Pb、Zn 8種金屬元素均達(dá)到中等污染。8種重金屬元素污染程度排序如下:Zn>Hg>As>Ni>Cu>Cr>Pb>Cd。其中,西南岸 XWC 站位 8種重金屬元素的單因

    表4 東海島紅樹林濕地土壤重金屬單因子污染指數(shù)Table 4 The single factor pollution index of heavy metals in the surface soil of mangrove wetland in Donghai Island

    圖2 東海島7個(gè)采樣站位主要超標(biāo)重金屬元素質(zhì)量分?jǐn)?shù)與國(guó)家土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)及環(huán)境背景值的對(duì)比Fig. 2 Comparison of the mass fraction of the excessive heavy metals in seven sampling stations in Donghai Island with national soil environment quality level standard and environmental background values

    子污染指數(shù)最高,Hg為很強(qiáng)污染程度,As和Zn為強(qiáng)污染程度;Cd、Cr、Cu、Ni和Pb 5種元素為中等污染程度。東南岸JJC和SYC站位8種元素的單因子污染指數(shù)普遍較低,這兩個(gè)區(qū)域污染程度相近,除As、Hg、Zn 3種元素達(dá)到中等污染程度級(jí)別外,其余5種元素均表現(xiàn)為輕微污染級(jí)別。

    從超標(biāo)狀況比較嚴(yán)重的As、Hg、Pb、Zn 4種重金屬元素污染程度空間分布看(圖3),其污染程度與元素質(zhì)量分?jǐn)?shù)的高低表現(xiàn)一致。其中,As在東海島西部采樣站位XWC和TMHHK污染最嚴(yán)重,為強(qiáng)污染,其他5個(gè)采樣站位為中等污染;Hg以西南岸XWC站位最嚴(yán)重,為很強(qiáng)污染,其他6個(gè)采樣站位為中等污染;Pb除在SYC、JJC和DCMT這3個(gè)采樣站位為輕微污染外,其他4個(gè)采樣站位為中等污染;Zn在所有紅樹林采樣站位中的污染程度普遍最高,其中YC、XWC和TMHHK為強(qiáng)污染,其他4個(gè)采樣站位為中等污染。

    通過野外實(shí)地調(diào)查及對(duì)比土壤環(huán)境背景值發(fā)現(xiàn),東海島紅樹林濕地土壤重金屬污染程度與土壤重金屬背景值高低及周圍人為活動(dòng)排放重金屬程度密切相關(guān)。采樣站位XWC和TMHHK中這4種重金屬元素污染均較為明顯,主要與臨近居民點(diǎn)和養(yǎng)殖魚塘有關(guān);XWC站位野外現(xiàn)場(chǎng)可見有多處生活垃圾傾倒點(diǎn)及養(yǎng)殖魚塘排污口,而TMHHK站位則位于入海河口,周圍分布著大量養(yǎng)殖魚塘,且東海島西部以農(nóng)業(yè)和水產(chǎn)養(yǎng)殖業(yè)為主,入海河口成為各類污染物的匯集處(李柳強(qiáng)等,2008)。采樣站位YC的Zn為強(qiáng)污染,其余3種元素為中等污染,主要與該采樣區(qū)農(nóng)業(yè)種植和水產(chǎn)養(yǎng)殖業(yè)有關(guān)。由此推測(cè),東海島紅樹林土壤重金屬污染主要來自水產(chǎn)養(yǎng)殖排污、農(nóng)業(yè)面源污染及生活排污,這與劉芳文等(2015)研究結(jié)果一致。東海島由養(yǎng)殖排污、農(nóng)業(yè)面源污染及生活排污引起的重金屬污染問題應(yīng)當(dāng)引起重視,應(yīng)注意避免人為直接排污及合理施用農(nóng)藥化肥。

    2.2.2 表層土壤重金屬潛在生態(tài)風(fēng)險(xiǎn)特征及空間分布

    圖3 東海島紅樹林濕地表層土壤主要超標(biāo)重金屬元素污染程度空間分布Fig. 3 Spatial distribution of excessive heavy metals pollution in the surface soil of mangrove wetland in Donghai Island

    表5 東海島紅樹林濕地土壤重金屬潛在生態(tài)風(fēng)險(xiǎn)參數(shù)和潛在生態(tài)風(fēng)險(xiǎn)指數(shù)Table 5 The potential ecological risk facter and potential ecological risk index of heavy metals in the surface soil of the mangrove wetland in Donghai Island

    東海島紅樹林濕地表層土壤重金屬生態(tài)風(fēng)險(xiǎn)危害指數(shù)計(jì)算結(jié)果如表5所示。從單個(gè)重金屬元素的潛在生態(tài)風(fēng)險(xiǎn)指數(shù)看,8種重金屬元素中 Hg潛在生態(tài)風(fēng)險(xiǎn)指數(shù)最高,平均值在 80~160之間,屬于強(qiáng)潛在生態(tài)風(fēng)險(xiǎn);其余7種重金屬元素屬于輕微潛在生態(tài)風(fēng)險(xiǎn)。各重金屬元素的平均潛在生態(tài)風(fēng)險(xiǎn)程度為 Hg>Cd>As>Ni>Cu>Pb>Zn>Cr。多種重金屬生態(tài)風(fēng)險(xiǎn)指數(shù)值在 112.002~351.690之間,平均值為 190.070,達(dá)到中等潛在生態(tài)風(fēng)險(xiǎn)。東海島紅樹林濕地表層土壤重金屬總潛在生態(tài)風(fēng)險(xiǎn)程度空間分布如圖4所示,不同站位的紅樹林濕地土壤重金屬污染的潛在生態(tài)風(fēng)險(xiǎn)指數(shù)有所差異。其中,西南岸XWC站位的總潛在生態(tài)風(fēng)險(xiǎn)程度最高,潛在生態(tài)風(fēng)險(xiǎn)指數(shù)平均值為351.69,為強(qiáng)生態(tài)風(fēng)險(xiǎn);YC、DCMT、DTSD和TMHHK的總潛在生態(tài)風(fēng)險(xiǎn)程度次之,潛在生態(tài)風(fēng)險(xiǎn)指數(shù)均值超過 150,為中等生態(tài)風(fēng)險(xiǎn);東南岸SYC和JJC站位的總潛在生態(tài)風(fēng)險(xiǎn)程度最低,為輕微生態(tài)風(fēng)險(xiǎn)。綜合潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果看,東海島各站位紅樹林土壤重金屬污染的潛在生態(tài)風(fēng)險(xiǎn)屬于輕微-強(qiáng)生態(tài)風(fēng)險(xiǎn)。其中西南岸XWC站位的強(qiáng)潛在生態(tài)風(fēng)險(xiǎn)主要來自 Hg、Cd和As污染,這3種重金屬為生物非必需元素,具有強(qiáng)毒性。結(jié)合XWC站位周圍環(huán)境,可以推測(cè)該區(qū)域強(qiáng)生態(tài)風(fēng)險(xiǎn)主要來自生活垃圾排放以及水產(chǎn)養(yǎng)殖業(yè)中使用的藥物、飼料等,如飼料添加劑的廣泛使用或含Hg、As農(nóng)用化學(xué)品的使用會(huì)使土壤受到污染(邵學(xué)新等,2007);而研究區(qū)也常見船舶往來,船舶排污、潤(rùn)滑油泄漏、機(jī)械部件磨損等會(huì)導(dǎo)致Cd含量的升高(冉延平,2012)。

    2.3 土壤重金屬來源分析

    2.3.1 相關(guān)性分析

    在研究區(qū)內(nèi),紅樹林濕地表層土壤中各重金屬元素含量及其之間的比率具有相對(duì)的穩(wěn)定性,當(dāng)沉積物來源相似或相同時(shí),其各元素之間具有顯著的相關(guān)性(韓卓汝等,2013)。不同站位間,重金屬元素的相關(guān)性是否顯著,可以反映各站位沉積環(huán)境的相似性和受人為影響程度的強(qiáng)弱(陳守莉等,2007)。從8種土壤重金屬之間的相關(guān)系數(shù)(表6)可以看出,除Hg、Cd以外,其他6種重金屬元素(As、Cr、Cu、Ni、Pb和 Zn)之間表現(xiàn)為強(qiáng)相關(guān)性,相關(guān)系數(shù)均在0.7以上。尤其是Cr-Ni、Cr-Pb、Cu-Ni、Cu-Zn、Ni-Pb這5組元素相關(guān)性極為顯著,相關(guān)系數(shù)分別達(dá)到 0.983、0.977、0.971、0.985、0.978(P<0.01),說明這些重金屬元素可能具有相近或相同的污染源。而Hg、Cd與其他重金屬元素相關(guān)性不明顯,一方面反映了Hg、Cd表現(xiàn)出與其他重金屬元素不同的地球化學(xué)行為特征,另一方面也反映了其與其他重金屬元素之間有著不同來源。

    表6 東海島紅樹林濕地土壤不同重金屬元素間的相關(guān)系數(shù)Table 6 Correlation coefficient between different heavy metal elements in the surface soil of mangrove wetland in Donghai Island

    2.3.2 主成分分析

    通過主成分分析(表7),紅樹林濕地表層土壤沉積物中8種重金屬污染物(8個(gè)變量)可由2個(gè)主成分(特征值:5.976+1.095=7.071個(gè)變量)反映88.381%的信息,即對(duì)前 2個(gè)主成分進(jìn)行分析大致能夠反映數(shù)據(jù)的大部分信息。

    表7 主成分分析主要計(jì)算結(jié)果Table 7 The main calculated results of principal component analysis (PCA)

    第一主成分的貢獻(xiàn)率為68.988%,特點(diǎn)表現(xiàn)為因子變量在As、Cr、Cu、Ni、Pb、Zn的濃度上有較高的正載荷,說明這些重金屬元素有相同或相似的物質(zhì)來源。野外實(shí)地調(diào)查發(fā)現(xiàn),XWC、TMHHK和DTSD站位周圍有大量的水產(chǎn)養(yǎng)殖池塘,YC站位附近主要為大面積的農(nóng)業(yè)種植區(qū)和少量的養(yǎng)殖魚塘,這4個(gè)站位均可見養(yǎng)殖排污口和農(nóng)藥化肥使用的殘留痕跡;從圖3和表4中也可以看出,As、Pb、Zn分別在YC、XWC、TMHHK和DTSD這4站位的污染程度都比較高。從對(duì)第一主成分的貢獻(xiàn)率上可知,As、Cr、Cu、Ni、Pb、Zn的貢獻(xiàn)率都超過90%,同時(shí)這6種元素之間表現(xiàn)為強(qiáng)相關(guān)性。因此,可把As、Cr、Cu、Ni、Pb、Zn這6種重金屬元素理解為第一主成分所代表的實(shí)際意義,即反映了農(nóng)業(yè)面源污染和養(yǎng)殖排污對(duì)東海島紅樹林濕地表層土壤的污染(韓卓汝等,2013;蔡龍炎,2010)。

    第二主成分的貢獻(xiàn)率為19.393%,該因子包含的變量有Cd和Hg,相應(yīng)的因子載荷值為0.795和0.619,但相關(guān)性分析顯示Cd和Hg之間的相關(guān)性不明顯,說明Cd和Hg的來源可能不同。野外實(shí)地調(diào)查發(fā)現(xiàn),YC和 XWC站位周圍常有漁船來往,而DCMT和DTSD則位于碼頭附近,從表4中可以看出,這4個(gè)站位的Cd污染程度比其他3個(gè)站位的高,故推測(cè)Cd主要可能來源于船舶排污。這與冉延平(2012)研究結(jié)果相吻合,其研究認(rèn)為靠近港口的研究區(qū)常有船舶往來,船舶排污、漏機(jī)械部件磨損、漏油等會(huì)導(dǎo)致 Cd污染。一般認(rèn)為 Hg可能與東海島的鋼鐵冶煉有關(guān)(雷凌明等,2014),但從圖3可以看出,DCMT和DTSD這兩個(gè)站位比其他 5個(gè)站位更靠近鋼鐵廠(位于東海島北部),而DCMT和DTSD的Hg污染程度卻不是最嚴(yán)重的,因此推測(cè)目前紅樹林濕地土壤Hg元素的污染可能與鋼鐵冶煉無關(guān)。通過野外實(shí)地調(diào)查發(fā)現(xiàn),YC、SYC、JJC、XWC、TMHHK、DCMT和DTSD這7個(gè)站位共同特征為可見生活垃圾自然排放,故推測(cè)Hg污染可能來自于生活排污(銀燕春等,2015)。另外Cd和Hg在第一主成分中的因子載荷為0.347和0.467,說明部分的Cd和Hg還受農(nóng)業(yè)面源污染、養(yǎng)殖排污或其他類型物源的影響。

    3 結(jié)論

    本研究對(duì)湛江東海島紅樹林濕地表層土壤重金屬元素進(jìn)行質(zhì)量分?jǐn)?shù)測(cè)定,采用潛在生態(tài)危害指數(shù)法進(jìn)行污染及潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià),結(jié)合相關(guān)性分析和主成分分析方法對(duì)污染源進(jìn)行分析,得到以下結(jié)論:

    (1)東海島紅樹林濕地表層土壤中8種重金屬質(zhì)量分?jǐn)?shù)平均值均超過雷州半島土壤環(huán)境背景值,未超過國(guó)家土壤環(huán)境質(zhì)量一級(jí)標(biāo)準(zhǔn)限值。7個(gè)采樣站位重金屬質(zhì)量分?jǐn)?shù)空間分布差異較大,各金屬變異程度為 Hg>Ni>As>Cd>Zn>Cr>Cu>Pb。

    (2)單因子污染指數(shù)顯示,東海島紅樹林濕地表層土壤各重金屬元素總體上屬于中等污染。從空間分布特征看,表層土壤重金屬污染程度與重金屬元素質(zhì)量分?jǐn)?shù)的高低較為一致。

    (3)不同站位的紅樹林濕地表層土壤重金屬潛在生態(tài)風(fēng)險(xiǎn)有所差異,屬于輕微至強(qiáng)生態(tài)風(fēng)險(xiǎn);重金屬污染總體上達(dá)到中等潛在生態(tài)風(fēng)險(xiǎn)。其中XWC站位重金屬污染為強(qiáng)潛在生態(tài)風(fēng)險(xiǎn),主要污染元素為Hg、Cd和As,推測(cè)污染主要來自生活排污、船舶污染和養(yǎng)殖排污。

    (4)統(tǒng)計(jì)分析結(jié)果表明,東海島紅樹林濕地表層土壤中 As、Cr、Cu、Ni、Pb和Zn這 6種重金屬元素之間具有強(qiáng)相關(guān)性,說明這些重金屬元素可能具有相近或相同的污染源,結(jié)合實(shí)地調(diào)查結(jié)果推測(cè)其來源主要與水產(chǎn)養(yǎng)殖排污及農(nóng)業(yè)面源污染有關(guān);而Cd和Hg之間的相關(guān)性不明顯,推測(cè)其來源分別與船舶污染和生活排污有關(guān)。

    猜你喜歡
    紅樹林金屬元素站位
    藏著寶藏的紅樹林
    固體廢物中金屬元素檢測(cè)方法研究進(jìn)展
    云南化工(2021年8期)2021-12-21 06:37:08
    提高政治站位 對(duì)標(biāo)國(guó)內(nèi)一流
    建黨百年說“站位”
    水上消防(2021年3期)2021-08-21 03:12:00
    提升站位講政治 創(chuàng)新?lián)?dāng)爭(zhēng)出彩
    微波消解-ICP-MS法同時(shí)測(cè)定牛蒡子中8種重金屬元素
    中成藥(2018年11期)2018-11-24 02:57:28
    神奇的紅樹林
    走過紅樹林
    歌海(2018年4期)2018-05-14 12:46:15
    環(huán)境監(jiān)測(cè)中重金屬元素分析方法探析
    河南科技(2014年8期)2014-02-27 14:08:07
    微波消解-HR-CS GFAAS法快速順序測(cè)定秀珍菇中金屬元素
    在线播放无遮挡| 男女国产视频网站| 18禁裸乳无遮挡免费网站照片| 婷婷色av中文字幕| 内射极品少妇av片p| 22中文网久久字幕| 国产一区亚洲一区在线观看| videossex国产| 欧美成人精品欧美一级黄| 亚洲人成网站在线观看播放| 亚洲激情五月婷婷啪啪| 中文精品一卡2卡3卡4更新| 日本av手机在线免费观看| 日本爱情动作片www.在线观看| 99热6这里只有精品| 夜夜爽夜夜爽视频| 麻豆久久精品国产亚洲av| 免费电影在线观看免费观看| 一边摸一边抽搐一进一小说| 国产亚洲av片在线观看秒播厂 | 亚洲av中文字字幕乱码综合| 一区二区三区高清视频在线| 欧美精品国产亚洲| 亚洲国产精品成人久久小说| 黄色欧美视频在线观看| 成人毛片60女人毛片免费| 日韩欧美三级三区| 久久久色成人| 久久这里有精品视频免费| 国产精品一区二区性色av| 亚洲欧洲国产日韩| 97热精品久久久久久| 嘟嘟电影网在线观看| 天堂影院成人在线观看| 丝袜喷水一区| 精品久久久久久久人妻蜜臀av| 国产精品.久久久| 亚洲中文字幕日韩| 男的添女的下面高潮视频| 国产真实伦视频高清在线观看| 亚洲国产精品sss在线观看| 久久精品国产自在天天线| 看非洲黑人一级黄片| 色综合亚洲欧美另类图片| 国产精品人妻久久久影院| 中文字幕精品亚洲无线码一区| 久久精品夜色国产| 超碰97精品在线观看| 蜜桃久久精品国产亚洲av| 亚洲人成网站高清观看| 成人高潮视频无遮挡免费网站| 亚洲精品aⅴ在线观看| 麻豆久久精品国产亚洲av| 男女啪啪激烈高潮av片| 欧美性猛交黑人性爽| 久久久色成人| 国产精品久久久久久久久免| 日本一本二区三区精品| 日韩在线高清观看一区二区三区| 伦理电影大哥的女人| 在线观看av片永久免费下载| 亚洲成人中文字幕在线播放| 亚洲欧洲国产日韩| 长腿黑丝高跟| 欧美性猛交╳xxx乱大交人| 99久久无色码亚洲精品果冻| 大香蕉97超碰在线| 国产麻豆成人av免费视频| 久久久国产成人精品二区| av线在线观看网站| 久久久精品大字幕| 国产伦一二天堂av在线观看| 国产精品一区二区在线观看99 | 国产在线男女| 久久久午夜欧美精品| 成人特级av手机在线观看| 国产麻豆成人av免费视频| 久久久久久久久中文| 亚洲精品乱码久久久久久按摩| 精品久久久久久久久av| a级毛色黄片| 亚洲国产日韩欧美精品在线观看| av国产免费在线观看| 久久综合国产亚洲精品| 超碰97精品在线观看| 国产午夜精品论理片| 国产亚洲5aaaaa淫片| 亚洲最大成人av| 久久精品国产亚洲网站| 麻豆成人av视频| 亚洲欧美日韩无卡精品| 国产精品久久久久久av不卡| 99热这里只有是精品50| 天天躁夜夜躁狠狠久久av| 97在线视频观看| 亚洲成人久久爱视频| 亚洲人成网站在线观看播放| 婷婷色综合大香蕉| 欧美成人a在线观看| 久久精品国产亚洲av天美| 三级经典国产精品| 亚洲av福利一区| 亚洲欧美中文字幕日韩二区| 97超碰精品成人国产| 只有这里有精品99| 18禁在线无遮挡免费观看视频| 精品人妻一区二区三区麻豆| 91午夜精品亚洲一区二区三区| 一级毛片电影观看 | 中文亚洲av片在线观看爽| 在线免费观看的www视频| 春色校园在线视频观看| 亚洲成av人片在线播放无| 2021天堂中文幕一二区在线观| 国产黄色小视频在线观看| 免费黄网站久久成人精品| 亚洲成色77777| 久久久久久伊人网av| 男人和女人高潮做爰伦理| 可以在线观看毛片的网站| 日韩精品青青久久久久久| 国产精品精品国产色婷婷| 亚洲欧美精品综合久久99| av福利片在线观看| 国国产精品蜜臀av免费| 色视频www国产| 免费看光身美女| 亚洲欧美成人综合另类久久久 | 波多野结衣高清无吗| 乱人视频在线观看| 欧美成人午夜免费资源| 22中文网久久字幕| 丝袜美腿在线中文| 亚洲人成网站高清观看| 亚洲真实伦在线观看| 久久精品久久精品一区二区三区| 成人美女网站在线观看视频| 亚洲精品aⅴ在线观看| 国产高清国产精品国产三级 | 精品午夜福利在线看| 日韩 亚洲 欧美在线| 久久久久久久久久成人| 日韩三级伦理在线观看| 干丝袜人妻中文字幕| 岛国毛片在线播放| 亚洲熟妇中文字幕五十中出| 91精品国产九色| 国产高清不卡午夜福利| or卡值多少钱| 中文字幕人妻熟人妻熟丝袜美| 麻豆久久精品国产亚洲av| 国产高清国产精品国产三级 | 又粗又硬又长又爽又黄的视频| 亚洲人与动物交配视频| 欧美不卡视频在线免费观看| 男女边吃奶边做爰视频| 一级av片app| 成人亚洲欧美一区二区av| 长腿黑丝高跟| 成人鲁丝片一二三区免费| 国产精品国产高清国产av| 国产私拍福利视频在线观看| 干丝袜人妻中文字幕| 亚洲成人av在线免费| 床上黄色一级片| 麻豆av噜噜一区二区三区| 国产白丝娇喘喷水9色精品| 国产精品人妻久久久影院| 99热全是精品| 亚洲av二区三区四区| 亚洲国产欧洲综合997久久,| 国产91av在线免费观看| 午夜福利网站1000一区二区三区| 欧美区成人在线视频| 中文精品一卡2卡3卡4更新| 美女高潮的动态| 色网站视频免费| 国产黄a三级三级三级人| 成人二区视频| 1000部很黄的大片| 久久韩国三级中文字幕| 精品久久久久久久人妻蜜臀av| 黄片无遮挡物在线观看| 国内揄拍国产精品人妻在线| 美女黄网站色视频| 三级国产精品欧美在线观看| 深夜a级毛片| 青春草亚洲视频在线观看| 亚洲国产色片| 观看免费一级毛片| 免费看光身美女| 日韩一本色道免费dvd| 七月丁香在线播放| 精品久久久久久久久av| 男人的好看免费观看在线视频| 精品人妻一区二区三区麻豆| 九草在线视频观看| 美女大奶头视频| 久久99热这里只频精品6学生 | 水蜜桃什么品种好| 成年av动漫网址| 大香蕉久久网| 国产不卡一卡二| 久久午夜福利片| 国产成人精品一,二区| 在线观看66精品国产| 亚洲欧美成人综合另类久久久 | 三级男女做爰猛烈吃奶摸视频| 美女黄网站色视频| 乱系列少妇在线播放| 99久久精品一区二区三区| 亚洲中文字幕一区二区三区有码在线看| 免费不卡的大黄色大毛片视频在线观看 | 黑人高潮一二区| 超碰av人人做人人爽久久| 久99久视频精品免费| 精品99又大又爽又粗少妇毛片| 久久久久国产网址| 波多野结衣高清无吗| 丰满乱子伦码专区| 国产精品一区www在线观看| 国产大屁股一区二区在线视频| 国产午夜精品论理片| 男的添女的下面高潮视频| 少妇的逼好多水| 国模一区二区三区四区视频| 成人二区视频| 久久鲁丝午夜福利片| 亚洲精品日韩av片在线观看| 亚洲av电影在线观看一区二区三区 | 国产老妇伦熟女老妇高清| 黑人高潮一二区| 亚洲av日韩在线播放| 亚洲av男天堂| 一个人观看的视频www高清免费观看| 如何舔出高潮| 国产免费男女视频| 免费看a级黄色片| 中文字幕熟女人妻在线| 久久精品夜色国产| 久久久久免费精品人妻一区二区| 日韩欧美精品v在线| 超碰av人人做人人爽久久| 99热全是精品| 18+在线观看网站| 婷婷色麻豆天堂久久 | 国产精品久久电影中文字幕| 看十八女毛片水多多多| 我的女老师完整版在线观看| 精品一区二区三区视频在线| 91狼人影院| 草草在线视频免费看| 六月丁香七月| 乱码一卡2卡4卡精品| 波多野结衣巨乳人妻| 久久久久精品人妻al黑| 国产精品秋霞免费鲁丝片| 青青草视频在线视频观看| 精品久久久精品久久久| 女的被弄到高潮叫床怎么办| 午夜免费观看性视频| 一级片免费观看大全| 黄色怎么调成土黄色| 国产深夜福利视频在线观看| 中文字幕人妻丝袜制服| 欧美国产精品一级二级三级| 亚洲精品中文字幕在线视频| 欧美最新免费一区二区三区| 一本—道久久a久久精品蜜桃钙片| 亚洲精品av麻豆狂野| 久久久久久人妻| 亚洲国产精品一区二区三区在线| 一级片'在线观看视频| 一级黄片播放器| 欧美最新免费一区二区三区| 大片电影免费在线观看免费| 婷婷色av中文字幕| 免费观看在线日韩| 国产一区二区在线观看av| 婷婷色综合www| 少妇精品久久久久久久| 老司机影院毛片| 久久久欧美国产精品| 人人妻人人澡人人爽人人夜夜| 多毛熟女@视频| 久热久热在线精品观看| 久久鲁丝午夜福利片| 日韩中字成人| 午夜91福利影院| 精品人妻在线不人妻| 99香蕉大伊视频| 亚洲高清免费不卡视频| 自线自在国产av| 国产 一区精品| 啦啦啦中文免费视频观看日本| 高清毛片免费看| 波野结衣二区三区在线| 日日撸夜夜添| 女性被躁到高潮视频| 欧美变态另类bdsm刘玥| 日本与韩国留学比较| 日本爱情动作片www.在线观看| 久久精品久久久久久久性| 国产探花极品一区二区| 美女主播在线视频| 精品国产国语对白av| 亚洲三级黄色毛片| 欧美亚洲 丝袜 人妻 在线| 亚洲精品av麻豆狂野| 午夜福利视频在线观看免费| 在线看a的网站| 欧美另类一区| 一边摸一边做爽爽视频免费| 美女大奶头黄色视频| 婷婷色综合www| 亚洲成色77777| 色视频在线一区二区三区| 久久狼人影院| av免费在线看不卡| 午夜激情av网站| 亚洲经典国产精华液单| 少妇被粗大猛烈的视频| 黑人欧美特级aaaaaa片| 亚洲国产最新在线播放| 飞空精品影院首页| 亚洲人成77777在线视频| 久久久久视频综合| 国产片特级美女逼逼视频| 国产又色又爽无遮挡免| 母亲3免费完整高清在线观看 | 久久鲁丝午夜福利片| 美女国产视频在线观看| 中文乱码字字幕精品一区二区三区| 午夜免费观看性视频| 亚洲精品久久久久久婷婷小说| xxxhd国产人妻xxx| 久久久国产欧美日韩av| 免费人妻精品一区二区三区视频| 美女福利国产在线| 亚洲国产毛片av蜜桃av| 丝袜人妻中文字幕| 少妇 在线观看| 美女大奶头黄色视频| 我要看黄色一级片免费的| 成人亚洲精品一区在线观看| 啦啦啦视频在线资源免费观看| 亚洲在久久综合| 精品一区在线观看国产| 日韩视频在线欧美| 99久久中文字幕三级久久日本| 久久久久久久大尺度免费视频| 亚洲精品,欧美精品| 亚洲人与动物交配视频| 1024视频免费在线观看| 欧美xxxx性猛交bbbb| 国产成人精品久久久久久| 丰满少妇做爰视频| av.在线天堂| 欧美日韩国产mv在线观看视频| 看十八女毛片水多多多| 久久午夜福利片| 欧美最新免费一区二区三区| 啦啦啦视频在线资源免费观看| 国产一区二区三区综合在线观看 | 精品国产一区二区三区久久久樱花| 亚洲 欧美一区二区三区| 久久av网站| 丰满迷人的少妇在线观看| 欧美老熟妇乱子伦牲交| 国产片内射在线| 亚洲精品久久成人aⅴ小说| 亚洲综合色网址| 国产色爽女视频免费观看| 亚洲伊人久久精品综合| av线在线观看网站| 色吧在线观看| 免费人成在线观看视频色| 18禁裸乳无遮挡动漫免费视频| 亚洲国产精品国产精品| 一本—道久久a久久精品蜜桃钙片| 午夜福利视频在线观看免费| 各种免费的搞黄视频| 高清在线视频一区二区三区| 七月丁香在线播放| 一级爰片在线观看| 国产男人的电影天堂91| 久久精品熟女亚洲av麻豆精品| 日韩av免费高清视频| 日韩精品免费视频一区二区三区 | 欧美成人午夜免费资源| 精品一区二区免费观看| 免费少妇av软件| 大片电影免费在线观看免费| 午夜福利乱码中文字幕| 搡老乐熟女国产| 国产毛片在线视频| xxxhd国产人妻xxx| 亚洲中文av在线| 9191精品国产免费久久| 亚洲欧美精品自产自拍| 亚洲国产精品成人久久小说| 亚洲激情五月婷婷啪啪| 亚洲成人一二三区av| 日本av免费视频播放| 久热久热在线精品观看| 亚洲伊人久久精品综合| 精品亚洲乱码少妇综合久久| av片东京热男人的天堂| 欧美精品高潮呻吟av久久| 亚洲综合色网址| 亚洲三级黄色毛片| √禁漫天堂资源中文www| 午夜视频国产福利| 激情五月婷婷亚洲| 日韩成人伦理影院| 看免费av毛片| 人体艺术视频欧美日本| 国产亚洲精品久久久com| 免费观看在线日韩| 最近2019中文字幕mv第一页| 午夜av观看不卡| 香蕉精品网在线| 视频区图区小说| 91成人精品电影| 亚洲综合色网址| 在线观看免费视频网站a站| 男女高潮啪啪啪动态图| 十八禁网站网址无遮挡| 亚洲 欧美一区二区三区| 90打野战视频偷拍视频| 女性被躁到高潮视频| 日韩人妻精品一区2区三区| 最近的中文字幕免费完整| 欧美亚洲 丝袜 人妻 在线| 香蕉精品网在线| 国产麻豆69| 亚洲欧美色中文字幕在线| 国产男女内射视频| 午夜福利,免费看| 日日爽夜夜爽网站| 你懂的网址亚洲精品在线观看| 在现免费观看毛片| 一级,二级,三级黄色视频| a级毛色黄片| 美女主播在线视频| 精品亚洲成国产av| 午夜久久久在线观看| 久久毛片免费看一区二区三区| 亚洲情色 制服丝袜| 国产成人精品福利久久| 午夜福利影视在线免费观看| 久久免费观看电影| 女人被躁到高潮嗷嗷叫费观| 最近最新中文字幕免费大全7| 成年动漫av网址| 国产视频首页在线观看| 亚洲中文av在线| 日韩免费高清中文字幕av| 22中文网久久字幕| 黄色 视频免费看| 亚洲av男天堂| 80岁老熟妇乱子伦牲交| 91在线精品国自产拍蜜月| 亚洲成人av在线免费| 国产片特级美女逼逼视频| 只有这里有精品99| 丝袜在线中文字幕| av网站免费在线观看视频| 青春草国产在线视频| 国产 一区精品| 黑人巨大精品欧美一区二区蜜桃 | 亚洲欧美色中文字幕在线| 欧美少妇被猛烈插入视频| 国产精品人妻久久久久久| 两个人看的免费小视频| 免费av不卡在线播放| 男人爽女人下面视频在线观看| 午夜日本视频在线| 建设人人有责人人尽责人人享有的| 街头女战士在线观看网站| 自线自在国产av| xxx大片免费视频| av视频免费观看在线观看| xxxhd国产人妻xxx| 免费女性裸体啪啪无遮挡网站| 国产精品一区www在线观看| 成人手机av| 亚洲国产精品国产精品| 高清av免费在线| 国产精品不卡视频一区二区| 王馨瑶露胸无遮挡在线观看| 午夜福利网站1000一区二区三区| 天堂俺去俺来也www色官网| 亚洲国产av新网站| 男女无遮挡免费网站观看| 一级毛片电影观看| 伦理电影大哥的女人| 国产亚洲精品久久久com| 一级,二级,三级黄色视频| 晚上一个人看的免费电影| 亚洲成国产人片在线观看| 日韩电影二区| 天堂8中文在线网| 精品国产一区二区三区久久久樱花| 精品一区二区三区视频在线| av播播在线观看一区| 男女午夜视频在线观看 | 高清视频免费观看一区二区| 一级毛片我不卡| 免费女性裸体啪啪无遮挡网站| 国产精品一区www在线观看| av免费观看日本| 久久久国产欧美日韩av| 国产亚洲最大av| 少妇高潮的动态图| a 毛片基地| 制服诱惑二区| 国产高清三级在线| 色视频在线一区二区三区| 黑丝袜美女国产一区| 亚洲人与动物交配视频| 99热国产这里只有精品6| 精品一区二区三卡| 国产日韩一区二区三区精品不卡| 99视频精品全部免费 在线| 水蜜桃什么品种好| 国产国语露脸激情在线看| 秋霞伦理黄片| 91在线精品国自产拍蜜月| 亚洲欧美精品自产自拍| 亚洲国产欧美日韩在线播放| 国产精品国产三级国产av玫瑰| 熟女av电影| 人妻人人澡人人爽人人| 视频在线观看一区二区三区| 国产在线一区二区三区精| 九色成人免费人妻av| 高清视频免费观看一区二区| 久久国内精品自在自线图片| 波多野结衣一区麻豆| 老女人水多毛片| 欧美性感艳星| 免费观看av网站的网址| 99热全是精品| 色吧在线观看| 国产毛片在线视频| 一本久久精品| 制服诱惑二区| 国产高清三级在线| 人人妻人人添人人爽欧美一区卜| 精品熟女少妇av免费看| 激情视频va一区二区三区| 18禁裸乳无遮挡动漫免费视频| 欧美激情国产日韩精品一区| 又黄又粗又硬又大视频| 免费观看a级毛片全部| 国产一区有黄有色的免费视频| 亚洲av电影在线观看一区二区三区| 香蕉丝袜av| a级毛色黄片| 亚洲av男天堂| 91精品三级在线观看| 五月天丁香电影| 少妇被粗大猛烈的视频| 午夜视频国产福利| 美女福利国产在线| 欧美亚洲 丝袜 人妻 在线| 日日啪夜夜爽| 美女xxoo啪啪120秒动态图| 日本免费在线观看一区| 麻豆乱淫一区二区| 久久久久久伊人网av| 亚洲,欧美,日韩| 99久国产av精品国产电影| 99香蕉大伊视频| 午夜免费观看性视频| 久久 成人 亚洲| 这个男人来自地球电影免费观看 | 亚洲成人一二三区av| 中文精品一卡2卡3卡4更新| 岛国毛片在线播放| 国产 精品1| 日韩制服骚丝袜av| kizo精华| 国产精品久久久久成人av| 一区在线观看完整版| 日韩视频在线欧美| 美女福利国产在线| 久久99精品国语久久久| 免费黄频网站在线观看国产| 免费观看性生交大片5| 永久网站在线| 欧美成人精品欧美一级黄| 久久久精品94久久精品| 插逼视频在线观看| 如何舔出高潮| 午夜久久久在线观看| 黑丝袜美女国产一区| 亚洲精品456在线播放app| 大片电影免费在线观看免费| 成人亚洲精品一区在线观看| 久久久精品区二区三区| 中文精品一卡2卡3卡4更新| 成人亚洲精品一区在线观看| 亚洲精品成人av观看孕妇| av又黄又爽大尺度在线免费看| 99国产精品免费福利视频| 五月天丁香电影| 亚洲伊人久久精品综合| 欧美激情国产日韩精品一区| 男男h啪啪无遮挡| 最近最新中文字幕大全免费视频 | 国产深夜福利视频在线观看| 亚洲国产精品专区欧美| 夜夜爽夜夜爽视频|