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

    2013年春夏季萊州灣海水環(huán)境要素特征和富營(yíng)養(yǎng)化評(píng)估

    2016-06-09 08:54:05徐艷東李佳蕙吳興偉馬元慶
    關(guān)鍵詞:評(píng)價(jià)

    徐艷東,魏 瀟,李佳蕙,吳興偉,馬元慶,孫 偉

    1.山東省海洋資源與環(huán)境研究院,山東省海洋生態(tài)修復(fù)重點(diǎn)實(shí)驗(yàn)室, 山東 煙臺(tái) 264006 2.中國(guó)海洋大學(xué)環(huán)境科學(xué)與工程學(xué)院, 山東 青島 266100

    2013年春夏季萊州灣海水環(huán)境要素特征和富營(yíng)養(yǎng)化評(píng)估

    徐艷東1,2,魏 瀟1,李佳蕙1,吳興偉1,馬元慶1,孫 偉1

    1.山東省海洋資源與環(huán)境研究院,山東省海洋生態(tài)修復(fù)重點(diǎn)實(shí)驗(yàn)室, 山東 煙臺(tái) 264006 2.中國(guó)海洋大學(xué)環(huán)境科學(xué)與工程學(xué)院, 山東 青島 266100

    根據(jù)2013年5(春季)、8月(夏季)萊州灣海水環(huán)境要素的調(diào)查資料,采用富營(yíng)養(yǎng)化指數(shù)、潛在性富營(yíng)養(yǎng)化評(píng)價(jià)模式和灰色聚類分析方法研究環(huán)境要素特征和評(píng)估海水富營(yíng)養(yǎng)化狀況。結(jié)果表明,無(wú)機(jī)氮是萊州灣水質(zhì)的主要污染要素,春夏季的N/P平均值分別為100.76、117.84,潛在性富營(yíng)養(yǎng)化評(píng)價(jià)模式結(jié)果表明,春夏季各站位的營(yíng)養(yǎng)級(jí)均只包括ⅣP、ⅥP兩類,磷限制為萊州灣的營(yíng)養(yǎng)鹽結(jié)構(gòu)特征;富營(yíng)養(yǎng)化指數(shù)評(píng)價(jià)結(jié)果表明,春季和夏季E>1站位比例分別為65%、20%;灰色聚類分析結(jié)果表明,春季Ⅱ級(jí)、Ⅲ級(jí)的站位比例分別為95%、5%,夏季Ⅱ、Ⅲ級(jí)的站位比例分別為70%、25%,Ⅱ級(jí)中的部分站位具有較大潛在富營(yíng)養(yǎng)化風(fēng)險(xiǎn)。

    富營(yíng)養(yǎng)化;環(huán)境要素;灰色聚類分析;氮磷比;萊州灣

    萊州灣位于山東半島西北部、渤海南部,是渤海的三大海灣之一,是中國(guó)典型的半封閉性陸架海灣,灣內(nèi)水深較淺,周邊分布有黃河、小清河、膠萊河等10余條河流,是黃渤海漁業(yè)生物的主要產(chǎn)卵場(chǎng)、棲息地和多種漁業(yè)的傳統(tǒng)漁場(chǎng)[1-2]。但隨著周邊地區(qū)經(jīng)濟(jì)的快速發(fā)展,萊州灣生態(tài)環(huán)境遭到破壞,海水氮磷比失衡現(xiàn)象較為明顯,赤潮時(shí)有發(fā)生[3]。目前,關(guān)于萊州灣富營(yíng)養(yǎng)化狀況已有不少報(bào)道[4-9],但主要局限于2010年之前的調(diào)查研究,對(duì)近年來萊州灣海水環(huán)境要素的系統(tǒng)調(diào)查和基于潛在性富營(yíng)養(yǎng)化評(píng)價(jià)模式和灰色聚類分析富營(yíng)養(yǎng)化的研究未見報(bào)道。本文以萊州灣2013年春夏季2個(gè)航次海水中氮磷營(yíng)養(yǎng)鹽及相關(guān)環(huán)境要素的調(diào)查數(shù)據(jù)為基礎(chǔ),采用富營(yíng)養(yǎng)化指數(shù)、潛在性富營(yíng)養(yǎng)化評(píng)價(jià)模式、灰色聚類分析方法研究環(huán)境要素特征和評(píng)估海水富營(yíng)養(yǎng)化狀況,為掌握該海域海水富營(yíng)養(yǎng)化形成機(jī)理和控制營(yíng)養(yǎng)鹽的總量、指導(dǎo)海水營(yíng)養(yǎng)鹽監(jiān)測(cè)方案優(yōu)化和指標(biāo)篩選提供技術(shù)支撐,對(duì)保護(hù)及修復(fù)脆弱的萊州灣海洋生態(tài)環(huán)境具有重要意義。

    1 材料和方法

    1.1 采樣站位、樣品采集和分析方法

    山東省海洋資源與環(huán)境研究院在萊州灣布設(shè)20個(gè)站位,2013年萊州灣水質(zhì)采樣站位見圖1。

    圖1 2013年萊州灣水質(zhì)采樣站位

    2013年5月(春季)、8月(夏季)開展了兩個(gè)航次的現(xiàn)場(chǎng)調(diào)查,指標(biāo)主要為溶解氧(DO)、化學(xué)需氧量(COD)、硝酸鹽氮(NO3-N)、亞硝酸鹽氮(NO2-N)、氨氮(NH3-N)、活性磷酸鹽(PO4-P)、葉綠素a(Chl-a)等。Chl-a只采集表層水樣,其他指標(biāo)在采樣水深大于10 m時(shí)采集表、底層水樣。樣品采集、貯存和運(yùn)輸按照《海洋監(jiān)測(cè)規(guī)范》第3部分:樣品采集、貯存和運(yùn)輸[10]的方法進(jìn)行。海水分析樣品的制備、消化和分析均按《海洋監(jiān)測(cè)規(guī)范》(GB 17378.4—2007)第4部分海水分析[11]的方法進(jìn)行。采樣過程中同時(shí)采集現(xiàn)場(chǎng)空白樣和

    現(xiàn)場(chǎng)平行樣,實(shí)驗(yàn)室分析采用平行樣、加標(biāo)回收樣、質(zhì)控樣分析等進(jìn)行質(zhì)量控制,以保證監(jiān)測(cè)數(shù)據(jù)的可靠性。

    1.2 數(shù)據(jù)分析和評(píng)價(jià)方法

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

    采用EXCEL、SPSS 20.0進(jìn)行數(shù)據(jù)處理和顯著性分析。取表、底層樣品的平均濃度代表該站位的實(shí)測(cè)濃度。各營(yíng)養(yǎng)鹽濃度季節(jié)間均值差異的顯著性檢驗(yàn)采用配對(duì)樣本t檢驗(yàn)方法,P<0.05表示差異顯著[12]。富營(yíng)養(yǎng)化空間分布圖的繪制采用Surfer 8.0軟件。

    1.2.2 富營(yíng)養(yǎng)化指數(shù)

    富營(yíng)養(yǎng)化指數(shù)是目前國(guó)內(nèi)最常用的一種富營(yíng)養(yǎng)化評(píng)價(jià)方法,其計(jì)算公式[13]為:

    (1)

    式中:E為富營(yíng)養(yǎng)化指數(shù),CCOD、CDIN、CDIP分別為COD、無(wú)機(jī)氮(DIN)、無(wú)機(jī)磷(DIP)的濃度,單位均為mg/L,其中DIP采用PO4-P計(jì)算。當(dāng)E≥1時(shí),表明水體呈現(xiàn)富營(yíng)養(yǎng)化,E值越大,富營(yíng)養(yǎng)化程度越嚴(yán)重。

    1.2.3 潛在性富營(yíng)養(yǎng)化評(píng)價(jià)模式

    REDFIELD[14]研究表明,浮游植物吸收氮磷摩爾比通常為16∶1(即Redfield值),浮游植物一般都按Redfield值來攝取營(yíng)養(yǎng)鹽,必然導(dǎo)致有一部分氮(對(duì)磷限制水體而言)或磷(對(duì)氮限制水體而言)相對(duì)過剩。郭衛(wèi)東等[15]認(rèn)為水體只有得到適量的氮(對(duì)磷限制水體而言)或磷(對(duì)氮限制水體而言)的補(bǔ)充,讓氮磷比(N/P)的值接近Redfield值,才能使這部分氮或磷對(duì)富營(yíng)養(yǎng)化的貢獻(xiàn)體現(xiàn)出來,這種現(xiàn)象被稱為潛在性富營(yíng)養(yǎng)化。郭衛(wèi)東等以此概念為基礎(chǔ),提出了潛在性富營(yíng)養(yǎng)化評(píng)價(jià)模式,具體劃分原則和富營(yíng)養(yǎng)化分級(jí)見表1。

    表1 潛在性富營(yíng)養(yǎng)化評(píng)價(jià)模式營(yíng)養(yǎng)級(jí)的劃分原則

    注:“—”表示該營(yíng)養(yǎng)級(jí)的判定不受該指標(biāo)限制。

    1.2.4 灰色聚類分析

    灰色聚類分析應(yīng)用灰色系統(tǒng)理論的思想和方法,根據(jù)關(guān)聯(lián)矩陣或灰數(shù)的白化權(quán)函數(shù)將一些調(diào)查指標(biāo)或監(jiān)測(cè)對(duì)象聚集成若干個(gè)可定義的類別,既能反映系統(tǒng)的灰色性,又能體現(xiàn)各指標(biāo)的綜合作用[16],可較客觀呈現(xiàn)海水富營(yíng)養(yǎng)化狀況,已被運(yùn)用到水體富營(yíng)養(yǎng)化的評(píng)價(jià)中[17-22]。本文選取DO、COD、Chl-a、DIN、PO4-P 5個(gè)指標(biāo)作為評(píng)價(jià)海水富營(yíng)養(yǎng)化的聚類指標(biāo),采用均值法對(duì)這5個(gè)指標(biāo)實(shí)測(cè)值和評(píng)價(jià)標(biāo)準(zhǔn)的灰類值進(jìn)行數(shù)據(jù)的無(wú)量綱處理,白化函數(shù)、聚類權(quán)和聚類系數(shù)的確定及構(gòu)造聚類向量判定各聚類對(duì)象等級(jí)的方法和公式,參見文獻(xiàn)[22]。5個(gè)指標(biāo)的評(píng)價(jià)標(biāo)準(zhǔn)和富營(yíng)養(yǎng)化的等級(jí)劃分主要參照《海水水質(zhì)標(biāo)準(zhǔn)》(GB 3097—1997)[23]和已有研究中對(duì)膠州灣富營(yíng)養(yǎng)化程度的分級(jí)標(biāo)準(zhǔn)及方法[17,22],灰色聚類分析的評(píng)價(jià)標(biāo)準(zhǔn)和等級(jí)如表2所示。最后依據(jù)最大隸屬度原則,即調(diào)查站位的富營(yíng)養(yǎng)化等級(jí)為其最大聚類系數(shù)所對(duì)應(yīng)的灰類。

    表2 灰色聚類分析的評(píng)價(jià)標(biāo)準(zhǔn)和等級(jí)

    2 結(jié)果和討論

    2.1 富營(yíng)養(yǎng)化指數(shù)計(jì)算結(jié)果與分析

    富營(yíng)養(yǎng)化指數(shù)計(jì)算結(jié)果和空間分布狀況(圖2)分析表明,該海域E值區(qū)間為0.18~3.10,平均值為1.13。春季E值范圍為0.43~3.00,有13個(gè)站位(占65%)E>1,主要分布于萊州灣西側(cè)和南側(cè)海域,其中最大值出現(xiàn)在萊州灣西南側(cè)的P18站位;夏季E值區(qū)間為0.18~3.10,4個(gè)站位(占20%)的E>1,主要位于萊州灣西北側(cè)海域,其中最大值出現(xiàn)在黃河口外側(cè)的P2站位;春夏季對(duì)比結(jié)果表明,春季富營(yíng)養(yǎng)化的范圍大于夏季。

    圖2 富營(yíng)養(yǎng)化指數(shù)的空間分布

    根據(jù)式(1)可知,E值的大小受COD、DIN、PO4-P濃度高低的影響。調(diào)查結(jié)果表明,COD春季的濃度范圍為1.59~1.98 mg/L,平均值為1.72 mg/L;夏季的濃度范圍為0.770~1.77 mg/L,平均值為1.30 mg/L;DIN春季的濃度范圍為0.220~0.612 mg/L,平均值為0.404 mg/L;夏季的濃度范圍為0.207~0.624 mg/L,平均值為0.360 mg/L。PO4-P春季濃度范圍為0.005 02~0.012 4 mg/L,平均值為0.008 96 mg/L;夏季的濃度范圍為0.003 77~0.014 6 mg/L,平均值為0.007 35 mg/L。與標(biāo)準(zhǔn)GB 3097—1997[23]對(duì)比結(jié)果表明,PO4-P、COD春夏季各站位均符合一類海水水質(zhì)標(biāo)準(zhǔn)。DIN春季劣于二類的站位比例達(dá)四分之三, 一半站位為四類和劣四類(符合二、三、四類及劣四類的站位比例分別為25%、25%、20%、30%);夏季劣于二類的站位比例也達(dá)七成, 三成站位達(dá)四類、劣四類(符合二、三、四類及劣四類標(biāo)準(zhǔn)的站位比例分別為30%、40%、20%、10%)。可見,DIN是調(diào)查海域的主要污染要素。春季劣四類區(qū)域分布在萊州灣西南部和黃河口鄰近海域,夏季分布在黃河口鄰近海域。60%的站位PO4-P濃度春季高于夏季,75%的站位DIN濃度春季略高于夏季,95%的站位COD濃度春季高于夏季。配對(duì)t檢驗(yàn)結(jié)果表明,PO4-P、COD濃度春季和夏季存在顯著差異,而DIN春季和夏季不存在顯著差異。究其原因,一方面是夏季水溫較春季高,浮游植物生長(zhǎng)旺盛,氮磷營(yíng)養(yǎng)物質(zhì)被大量消耗[24]。調(diào)查結(jié)果顯示,春季和夏季Chl-a的平均濃度分別為1.31、4.37 mg/m3;另外,可能與降雨有關(guān),2013年山東省平均降雨量5月(104.3 mm)較歷年同期偏多122.9%[25],而8月(80.1 mm)則較歷年同期偏少47.6%[26]。此外,大氣的氮沉降也是海洋中氮的重要來源[27]。

    2.2 潛在性富營(yíng)養(yǎng)化評(píng)價(jià)模式計(jì)算結(jié)果與分析

    潛在性富營(yíng)養(yǎng)化評(píng)價(jià)模式計(jì)算結(jié)果和富營(yíng)養(yǎng)化程度分布狀況(圖3)分析表明,春季和夏季各站位的營(yíng)養(yǎng)級(jí)均只包括磷限制中度營(yíng)養(yǎng)(ⅣP)、磷限制潛在性富營(yíng)養(yǎng)(ⅥP)兩類,未出現(xiàn)其他7類營(yíng)養(yǎng)級(jí)別,這說明磷限制為萊州灣營(yíng)養(yǎng)鹽結(jié)構(gòu)的主要特征。其中,夏季磷限制潛在性富營(yíng)養(yǎng)站位的比例略高于春季。春季的ⅣP、ⅥP站位比例分別為30%、70%,夏季的ⅣP、ⅥP站位比例分別為25%、75%。從分布上看,春夏季磷限制區(qū)域主要位于灣中部和東部。

    圖3 潛在性富營(yíng)養(yǎng)化評(píng)價(jià)模式結(jié)果站位分布圖

    潛在性富營(yíng)養(yǎng)化評(píng)價(jià)模式主要包括DIN、PO4-P濃度和N/P。調(diào)查結(jié)果表明,該海域春季的N/P的范圍為48.36~150.39,平均值為100.76,夏季N/P的范圍為48.05~205.26,平均值為117.84,兩個(gè)季節(jié)N/P的最小值均遠(yuǎn)大于Redfield值,其中N/P>100的站位比例分別為50%(春)和65%(夏),因此該海域的營(yíng)養(yǎng)鹽結(jié)構(gòu)為磷限制,夏季氮磷比失衡狀況比春季嚴(yán)重,各站位的N/P均大于30,依據(jù)潛在性富營(yíng)養(yǎng)化評(píng)價(jià)模式的營(yíng)養(yǎng)等級(jí)劃分標(biāo)準(zhǔn),均呈現(xiàn)為磷限制,這時(shí)其等級(jí)只與DIN濃度和N/P有關(guān)。出現(xiàn)磷限制現(xiàn)象的主要原因,可能是化肥用量劇增但比例不當(dāng),地表水把未被利用的過量氮肥匯入河水,而磷酸鹽從河流徑流的補(bǔ)充有限,致使入海徑流氮磷比很高[15],高氮磷比河水的注入必然會(huì)使河口及其鄰近海域水體中N/P升高。另外,也與大氣氮沉降[27]和水交換周期長(zhǎng)有關(guān)。氮濃度過高和氮磷比失衡會(huì)引起浮游植物種群結(jié)構(gòu)的改變或富營(yíng)養(yǎng)化,甚至引發(fā)赤潮[28]。萊州灣2001年的浮游植物多樣性物種組成和多樣性指數(shù)較1989年均發(fā)生了很大變化[4],1990—2010年期間,該海域共發(fā)生赤潮11次[3]。因此,應(yīng)加強(qiáng)富營(yíng)養(yǎng)化海域監(jiān)測(cè),特別關(guān)注營(yíng)養(yǎng)鹽失衡嚴(yán)重區(qū)域,采取措施減少氮的排放并加大氮移除力度,改善營(yíng)養(yǎng)鹽結(jié)構(gòu),避免富營(yíng)養(yǎng)化和營(yíng)養(yǎng)鹽結(jié)構(gòu)失衡程度的加劇。

    2.3 灰色聚類分析計(jì)算結(jié)果與分析

    基于灰色聚類的評(píng)價(jià)結(jié)果(表3)和富營(yíng)養(yǎng)化程度分布狀況(圖4)表明,春季有1個(gè)站位(P17,位于小清河口鄰近海域)Ⅲ級(jí)灰類的聚類系數(shù)最大,聚類結(jié)果為Ⅲ級(jí),表明這個(gè)區(qū)域水體污染較重,為富營(yíng)養(yǎng)水平。其余19個(gè)站位Ⅱ級(jí)灰類的聚類系數(shù)最大,聚類結(jié)果為Ⅱ級(jí),表明這些區(qū)域水體受到一定的污染,為中度富營(yíng)養(yǎng)水平,其中P2、P12、P18、P19 4個(gè)站位雖處于中度富營(yíng)養(yǎng)水平,但這些站位的Ⅲ級(jí)灰類的聚類系數(shù)卻較大,均大于0.400,與Ⅱ級(jí)的灰類的聚類系數(shù)相差很小,表明這些區(qū)域具有潛在的富營(yíng)養(yǎng)化風(fēng)險(xiǎn),應(yīng)加以關(guān)注。夏季有5個(gè)站位(P1、P2、P12、P16、P20,主要位于黃河口和小清河口鄰近海域及萊州灣東南部)Ⅲ級(jí)灰類的聚類系數(shù)最大,聚類結(jié)果為Ⅲ級(jí),表明這個(gè)區(qū)域水體污染較重,為富營(yíng)養(yǎng)水平;1個(gè)站位(P9)Ⅰ級(jí)灰類的聚類系數(shù)最大,聚類結(jié)果為Ⅰ級(jí),表明這個(gè)區(qū)域海水質(zhì)量較好,為貧營(yíng)養(yǎng)水平;其余14個(gè)站位Ⅱ級(jí)灰類的聚類系數(shù)最大,聚類結(jié)果為Ⅱ級(jí),表明這些區(qū)域水體受到一定的污染,為中度富營(yíng)養(yǎng)水平,其中P4、P13、P18 3個(gè)站位雖處于中度富營(yíng)養(yǎng)水平,但這些站位的Ⅲ級(jí)灰類的聚類系數(shù)也較大,均大于0.300,與Ⅱ級(jí)的灰類的聚類系數(shù)相差較小,表明這些區(qū)域具有潛在的富營(yíng)養(yǎng)化風(fēng)險(xiǎn),應(yīng)加以關(guān)注。

    表3 灰色聚類和富營(yíng)養(yǎng)化指數(shù)評(píng)價(jià)結(jié)果

    圖4 灰色聚類分析結(jié)果站位分布圖

    灰色聚類分析選取了DO、COD、Chl-a、DIN、PO4-P作為評(píng)價(jià)富營(yíng)養(yǎng)化的聚類指標(biāo),聚類結(jié)果受這5個(gè)指標(biāo)影響。調(diào)查結(jié)果表明,DO春季的濃度范圍為8.88~10.16 mg/L,平均值為9.16 mg/L;夏季的濃度范圍為5.20~7.41 mg/L,平均值為6.54 mg/L。Chl-a春季的濃度范圍為0.764~2.13 mg/m3,平均值為1.31 mg/m3;夏季的濃度范圍為0.460~9.11 mg/m3,平均值為4.37 mg/m3。所有站位的DO春季均高于夏季,85%的站位Chl-a春季低于夏季,配對(duì)t檢驗(yàn)結(jié)果表明,DO、Chl-a濃度春季和夏季存在顯著差異。這主要是由于夏季水溫較春季高,浮游植物生長(zhǎng)旺盛有關(guān)。5個(gè)指標(biāo)與灰色聚類分析的評(píng)價(jià)標(biāo)準(zhǔn)對(duì)比結(jié)果表明,各站位PO4-P濃度均低于Ⅰ級(jí)標(biāo)準(zhǔn),因此其聚類權(quán)重較低;DO、COD的監(jiān)測(cè)值低于Ⅰ級(jí)標(biāo)準(zhǔn)或介于Ⅰ~Ⅱ標(biāo)準(zhǔn),其聚類權(quán)重次之;Chl-a、DIN的監(jiān)測(cè)數(shù)值較大且變動(dòng)范圍大,故其聚類權(quán)重高,是影響聚類結(jié)果的主要指標(biāo)。由于DIN春夏季濃度不存在顯著差異,大部分站位的Chl-a濃度夏季高于春季,因此夏季富營(yíng)養(yǎng)級(jí)站位數(shù)多于春季。

    2.4 計(jì)算結(jié)果比較與討論

    富營(yíng)養(yǎng)化指數(shù)和灰色聚類分析結(jié)果相同之處是均篩選出春季萊州灣西南側(cè)呈富營(yíng)養(yǎng)化、夏季西北側(cè)的黃河口富營(yíng)養(yǎng)化。但依據(jù)海洋富營(yíng)養(yǎng)化是海水中營(yíng)養(yǎng)物質(zhì)過度增加并導(dǎo)致生態(tài)系統(tǒng)有機(jī)質(zhì)增多、低氧區(qū)形成、藻華暴發(fā)一些異常改變的過程[29]這一定義,富營(yíng)養(yǎng)化應(yīng)包括營(yíng)養(yǎng)鹽的增加且還應(yīng)包括生態(tài)系統(tǒng)的異常改變,因此其評(píng)價(jià)應(yīng)包括相關(guān)的生態(tài)指標(biāo)。但富營(yíng)養(yǎng)化指數(shù)公式不含有生態(tài)指標(biāo),且當(dāng)COD、PO4-P、DIN中某一指標(biāo)很高或很低時(shí),會(huì)影響計(jì)算結(jié)果的偏高或偏低,會(huì)掩蓋其他指標(biāo)的影響,由于萊州灣具有DIN濃度高的特點(diǎn),故萊州灣的富營(yíng)養(yǎng)化指數(shù)計(jì)算結(jié)果會(huì)夸大氮營(yíng)養(yǎng)鹽的作用[18,23]?;翌惥垲惙治鐾ㄟ^白化數(shù)據(jù)的灰化處理,克服了評(píng)價(jià)標(biāo)準(zhǔn)邊界過于明確對(duì)評(píng)價(jià)結(jié)果的影響,既可較好體現(xiàn)富營(yíng)養(yǎng)化程度的灰色性和調(diào)查數(shù)據(jù)的不確定性,又可反映單個(gè)指標(biāo)或多個(gè)指標(biāo)對(duì)于污染的貢獻(xiàn),并可通過各個(gè)指標(biāo)的聚類系數(shù)反映局部差異,能夠較客觀合理的對(duì)調(diào)查結(jié)果進(jìn)行綜合性評(píng)價(jià)[18-23],本研究中基于灰色聚類計(jì)算的富營(yíng)養(yǎng)化結(jié)果不僅確定了富營(yíng)養(yǎng)化站位,還篩選出存在較大潛在富營(yíng)養(yǎng)風(fēng)險(xiǎn)的站位,并且灰色聚類分析較富營(yíng)養(yǎng)化指數(shù)增加了可體現(xiàn)富營(yíng)養(yǎng)化的水中氧濃度水平的DO和生態(tài)指標(biāo)Chl-a等參與評(píng)價(jià),因此灰色聚類分析較富營(yíng)養(yǎng)化指數(shù)得出的結(jié)論更客觀。但這兩種方法均未考慮氮磷比,可能會(huì)忽略氮或磷限制表現(xiàn)出富營(yíng)養(yǎng)化的潛在性,應(yīng)結(jié)合營(yíng)養(yǎng)鹽的氮磷比結(jié)構(gòu)來確定富營(yíng)養(yǎng)化水平更合理。

    潛在性富營(yíng)養(yǎng)化評(píng)價(jià)模式考慮了氮磷比,但也和富營(yíng)養(yǎng)化指數(shù)一樣缺少體現(xiàn)富營(yíng)養(yǎng)出現(xiàn)時(shí)底層缺氧、浮游植物大量繁殖等現(xiàn)象的指標(biāo),并且對(duì)萊州灣海域的評(píng)價(jià)結(jié)果為磷限制中度營(yíng)養(yǎng)(ⅣP)和磷限制潛在性富營(yíng)養(yǎng)(ⅥP),結(jié)果均屬磷限制類,結(jié)論較粗,不符合海洋管理精細(xì)化的要求,由于萊州灣海域氮高磷低的特性,其評(píng)價(jià)結(jié)果可作為灰色聚類分析的補(bǔ)充。

    以上說明3種方法各有優(yōu)劣,當(dāng)前富營(yíng)養(yǎng)化的評(píng)估應(yīng)綜合運(yùn)用這些方法,從不同需求獲取有價(jià)值的信息應(yīng)用到海洋管理中。下一步需要研究提出一種計(jì)算結(jié)果更科學(xué)、結(jié)論更便于理解和應(yīng)用的富營(yíng)養(yǎng)化評(píng)估模型,以滿足日常海洋調(diào)查和監(jiān)測(cè)的需要。

    3 結(jié)論

    1) 春夏季PO4-P、COD各站位均符合一類海水水質(zhì)標(biāo)準(zhǔn),DIN春夏季劣四類的站位比例分別為30%、10%,DIN是調(diào)查海域的主要污染要素,劣四類區(qū)域春季分布在萊州灣西南部和黃河口鄰近海域、夏季分布在黃河口鄰近海域。春夏兩季的N/P平均值分別為100.76、117.84,N/P>100的站位比例分別為50%、65%,因此萊州灣的營(yíng)養(yǎng)鹽結(jié)構(gòu)為磷限制,夏季氮磷比失衡狀況較春季更為嚴(yán)重。

    2) 春夏季富營(yíng)養(yǎng)化指數(shù)的平均值為1.13,春夏季E>1站位比例分別為65%、20%,E>1區(qū)域春季主要分布在萊州灣西側(cè)和南側(cè)、夏季主要分布在萊州灣西北側(cè)海域,春季富營(yíng)養(yǎng)化的范圍大于夏季?;疑垲惙治鼋Y(jié)果表明,春季Ⅱ、Ⅲ級(jí)的站位比例分別為95%、5%,夏季Ⅱ、Ⅲ級(jí)的站位比例分別為70%、25%,Ⅱ級(jí)中的部分站位Ⅲ級(jí)灰類的聚類系數(shù)較大,具有較大潛在富營(yíng)養(yǎng)化風(fēng)險(xiǎn)。富營(yíng)養(yǎng)級(jí)(Ⅲ級(jí))和潛在富營(yíng)養(yǎng)化風(fēng)險(xiǎn)的區(qū)域春季主要位于萊州灣西部、夏季主要位于萊州灣西北部和萊州灣南部,夏季富營(yíng)養(yǎng)化范圍大于春季。潛在性富營(yíng)養(yǎng)化評(píng)價(jià)模式結(jié)果表明,春夏季各站位的營(yíng)養(yǎng)級(jí)均只包括ⅣP、ⅥP兩類,春夏季磷限制區(qū)域主要位于萊州灣中部和東部。

    3) 灰色聚類分析較富營(yíng)養(yǎng)化指數(shù)得出的結(jié)論更客觀,潛在性富營(yíng)養(yǎng)化評(píng)價(jià)模式的結(jié)果可作為灰色聚類分析的補(bǔ)充,需綜合運(yùn)用這些方法評(píng)估海水富營(yíng)養(yǎng)化,從不同需求獲取有價(jià)值的信息應(yīng)用到海洋管理中。應(yīng)加強(qiáng)富營(yíng)養(yǎng)化海域監(jiān)測(cè),特別關(guān)注營(yíng)養(yǎng)鹽失衡嚴(yán)重區(qū)域。

    [1] 中國(guó)海灣志編纂委員會(huì).中國(guó)海灣志:第三分冊(cè)[M].北京:海洋出版社,1991:1-77.

    [2] 鄧景耀,金顯仕.萊州灣及黃河口水域漁業(yè)生物多樣性及其保護(hù)研究[J].動(dòng)物學(xué)研究,2000,21(1):76-82.

    [3] 于寧,于建生,呂振波,等.山東海域赤潮災(zāi)害特征及預(yù)警報(bào)管理[J].生態(tài)學(xué)雜志,2012,31(5):1 272-1 281.

    [4] 郝彥菊,王宗靈,朱明遠(yuǎn),等.萊州灣營(yíng)養(yǎng)鹽與浮游植物多樣性調(diào)查與評(píng)價(jià)研究[J].海洋科學(xué)進(jìn)展,2005,23(2):197-204.

    [5] 孫丕喜,王波,張朝暉,等.萊州灣海水中營(yíng)養(yǎng)鹽分布與富營(yíng)養(yǎng)化的關(guān)系[J].海洋科學(xué)進(jìn)展,2006,24(3):329-335.

    [6] 夏斌,張曉理,崔毅,等.夏季萊州灣及附近水域理化環(huán)境及營(yíng)養(yǎng)現(xiàn)狀評(píng)價(jià)[J].漁業(yè)科學(xué)進(jìn)展,2009,30(3):103-111.

    [7] 劉義豪,楊秀蘭,靳洋,等.萊州灣海域營(yíng)養(yǎng)鹽現(xiàn)狀及年際變化規(guī)律[J].漁業(yè)科學(xué)進(jìn)展,2011,32(4):1-5.

    [8] 張雪,張龍軍,侯中里,等.1980—2008年萊州灣主要污染物的時(shí)空變化[J].中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版),2012,42(11):91-98.

    [9] 張錦峰,高學(xué)魯,莊文,等.萊州灣漁業(yè)資源與環(huán)境變化趨勢(shì)分析[J].海洋湖沼通報(bào),2014,3:82-90.

    [10] 國(guó)家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局,中國(guó)國(guó)家標(biāo)準(zhǔn)化管理委員會(huì).海洋監(jiān)測(cè)規(guī)范 第3部分:樣品采集、貯存與運(yùn)輸:GB 17378.3—2007[S].北京:中國(guó)標(biāo)準(zhǔn)出版社,2008:2-7.

    [11] 國(guó)家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局,中國(guó)國(guó)家標(biāo)準(zhǔn)化管理委員會(huì).海洋監(jiān)測(cè)規(guī)范 第4部分:海水分析:GB 17378.4—2007[S].北京:中國(guó)標(biāo)準(zhǔn)出版社,2008:1-162.

    [12] 張文彤,鄺春偉.SPSS統(tǒng)計(jì)分析基礎(chǔ)教程[M].2版.北京:高等教育出版社,2011:213-235.

    [13] 鄒景忠,董麗萍,秦保平.渤海灣富營(yíng)養(yǎng)化和赤潮問題的初步探討[J].海洋環(huán)境科學(xué),1983,2(2):41-54.

    [14] REDFIELD A C.The biological control of chemical factors in the environments[J].American Scientist,1958,46:205-221.

    [15] 郭衛(wèi)東,章小明,楊逸萍,等.中國(guó)近岸海域潛在性富營(yíng)養(yǎng)化程度的評(píng)價(jià)[J].臺(tái)灣海峽,1998,17(1):64-70.

    [16] 孫維萍,于培松,潘建明.灰色聚類法評(píng)價(jià)長(zhǎng)江口、杭州灣海域表層海水中的重金屬污染程度[J].海洋學(xué)報(bào),2009,31(1):79-84.

    [17] 姚云,沈志良.膠州灣海水富營(yíng)養(yǎng)化水平評(píng)價(jià)[J].海洋科學(xué),2004,28(6):14-17.

    [18] 陳鳴淵,俞志明,宋秀賢,等.利用模糊綜合方法評(píng)價(jià)長(zhǎng)江口海水富營(yíng)養(yǎng)化水平[J].海洋科學(xué),2007,31(11):47-54

    [19] 胡麗慧,潘安,李鐵松,等.灰色聚類法在升鐘水庫(kù)水體富營(yíng)養(yǎng)化評(píng)價(jià)中的應(yīng)用[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2008,27(6):2 407-2 412.

    [20] 吳祥慶,黎小正,秦振發(fā),等.灰色聚類法在西津水庫(kù)水體富營(yíng)養(yǎng)化評(píng)價(jià)中的應(yīng)用[J].安徽農(nóng)業(yè)科學(xué),2010,38(22):11 926-11 928.

    [21] 張學(xué)超,劉營(yíng),胡紅智,等.基于灰色聚類的威海雙島灣海域富營(yíng)養(yǎng)化水平評(píng)價(jià)[J].水產(chǎn)科學(xué),2013,32(10):605-608.

    [22] 夏斌,陳碧鵑,辛福言,等.灰色聚類法在膠州灣富營(yíng)養(yǎng)化評(píng)價(jià)中的應(yīng)用[J].漁業(yè)科學(xué)進(jìn)展,2011,32(5):114-120.

    [23] 國(guó)家環(huán)境保護(hù)局.海水水質(zhì)標(biāo)準(zhǔn):GB 3097—1997[S/OL].[2016-07-05].http://kjs.mep.gov.cn/hjbhbz/bzwb/shjbh/shjzlbz/199807/t19980701_66499. shtml.

    [24] 馮士筰,李鳳岐,李少菁.海洋科學(xué)導(dǎo)論[M].北京:高等教育出版社,1999:138-140.

    [25] 山東省水文局.山東省五月水情月報(bào)[R].濟(jì)南:山東省水文局,2013:1-2.

    [26] 山東省水文局.山東省八月水情月報(bào)[R].濟(jì)南:山東省水文局,2013:1-2.

    [27] 高會(huì)旺,姚小紅,郭志剛,等.大氣沉降對(duì)海洋初級(jí)生產(chǎn)過程與氮循環(huán)的影響研究進(jìn)展[J].地球科學(xué)進(jìn)展,2014,29(12):1 325-1 332.

    [28] 張輝,石曉勇,張傳松,等.北黃海營(yíng)養(yǎng)鹽結(jié)構(gòu)及限制作用時(shí)空分布特征分析[J].中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版),2009,39(4):773-780.

    [29] 俞志明,沈志良,陳亞瞿,等.長(zhǎng)江口水域富營(yíng)養(yǎng)化[M].北京:科學(xué)出版社,2011:1-12.

    Characteristics of Environment Parameters and Eutrophication Assessment of Seawater in Laizhou Bay in Spring and Summer, 2013

    XU Yandong1, 2,WEI Xiao1,LI Jiahui1,WU Xingwei1,MA Yuanqing1,SUN Wei1

    1.Shandong Marine Resource and Environment Research Institute, Shandong Provincial Key Laboratory of Restoration for Marine Ecology, Yantai 264006, China 2.College of Environment Science and Engineering, Ocean University of China, Qingdao 266100, China

    Based on field data of environmental parameters of seawater in Laizhou Bay in May (on behalf of spring) and August (on behalf of summer) of 2013, their characteristics were investigated and eutrophication degree of seawater was successively assessed by eutrophication index, potential eutrophication assessment model (PEAM) and grey clustering analysis (GCA). The results revealed that dissolved inorganic nitrogen (DIN) is the principal contamination parameter for water quality of Laizhou Bay. Mean values ofN∶Pratio in spring and summer are 100.76 and 117.84, respectively. By means of potential eutrophication assessment model, we found that there are merely two levels of eutrophication, ⅣPand ⅥP, to which all the sampling sites in spring and summer have attained. Thus, the structure of nutrients is characterized by phosphorus-limitation in Laizhou Bay. As for the eutrophication index,Evalues exceed 1, the guideline, for 65% of all the sampling sites in spring and for 20% of them in summer. The results of grey clustering analysis showed that the proportions of mid-eutrophic (Ⅱ) and eutrophic (Ⅲ) stations are 95% and 5% respectively in spring while are 70% and 25% respectively in summer, and some mid-eutrophic (Ⅱ) ones have comparatively high potential risk of eutrophication.

    eutrophication;environment parameters;grey clustering analysis;N/Pratio;Laizhou Bay

    2015-05-11;

    2015-07-22

    國(guó)家海洋公益性行業(yè)專項(xiàng)經(jīng)費(fèi)項(xiàng)目(201105006, 201205001);山東省科技發(fā)展計(jì)劃(2014GSF117030)

    徐艷東(1980-),男,山東蒼山人,在讀博士,副研究員。

    X824

    A

    1002-6002(2016)06- 0063- 07

    10.19316/j.issn.1002-6002.2016.06.10

    猜你喜歡
    評(píng)價(jià)
    SBR改性瀝青的穩(wěn)定性評(píng)價(jià)
    石油瀝青(2021年4期)2021-10-14 08:50:44
    中藥治療室性早搏系統(tǒng)評(píng)價(jià)再評(píng)價(jià)
    自制C肽質(zhì)控品及其性能評(píng)價(jià)
    寫作交流與評(píng)價(jià):詞的欣賞
    基于Moodle的學(xué)習(xí)評(píng)價(jià)
    關(guān)于項(xiàng)目后評(píng)價(jià)中“專項(xiàng)”后評(píng)價(jià)的探討
    HBV-DNA提取液I的配制和應(yīng)用評(píng)價(jià)
    有效評(píng)價(jià)讓每朵花兒都綻放
    模糊數(shù)學(xué)評(píng)價(jià)法在水質(zhì)評(píng)價(jià)中的應(yīng)用
    治淮(2013年1期)2013-03-11 20:05:18
    保加利亞轉(zhuǎn)軌20年評(píng)價(jià)
    亚洲无线在线观看| 久久伊人香网站| 精品欧美国产一区二区三| 亚洲aⅴ乱码一区二区在线播放 | 国产高清videossex| 丰满人妻熟妇乱又伦精品不卡| 国产精品一区二区免费欧美| 99精品在免费线老司机午夜| 91九色精品人成在线观看| 国产精品电影一区二区三区| 一个人免费在线观看的高清视频| 国产精品二区激情视频| 国产精品免费一区二区三区在线| www.999成人在线观看| 男女床上黄色一级片免费看| av在线天堂中文字幕| 欧美在线一区亚洲| 国产成人精品无人区| 免费看a级黄色片| 亚洲熟女毛片儿| 国产精品1区2区在线观看.| 好看av亚洲va欧美ⅴa在| 中文字幕精品免费在线观看视频| 欧美乱码精品一区二区三区| 国产精品国产高清国产av| 在线观看66精品国产| 亚洲 国产 在线| 欧美午夜高清在线| 99精品欧美一区二区三区四区| 免费在线观看完整版高清| 亚洲av成人av| а√天堂www在线а√下载| 久久精品成人免费网站| 久久午夜亚洲精品久久| 女人被狂操c到高潮| 成人18禁在线播放| 亚洲第一电影网av| 国产精品野战在线观看| 50天的宝宝边吃奶边哭怎么回事| 精品久久久久久成人av| 高清在线国产一区| 亚洲国产中文字幕在线视频| 午夜福利成人在线免费观看| 可以免费在线观看a视频的电影网站| svipshipincom国产片| 波多野结衣一区麻豆| 免费在线观看亚洲国产| 亚洲成a人片在线一区二区| 欧美成人免费av一区二区三区| 亚洲天堂国产精品一区在线| 中文字幕高清在线视频| 色婷婷久久久亚洲欧美| 国产午夜福利久久久久久| 男女午夜视频在线观看| 在线观看一区二区三区| 99热只有精品国产| 丰满的人妻完整版| 最新在线观看一区二区三区| 精品午夜福利视频在线观看一区| 操美女的视频在线观看| 亚洲天堂国产精品一区在线| 丝袜美足系列| 国产免费男女视频| 精品一区二区三区视频在线观看免费| 伊人久久大香线蕉亚洲五| 国产xxxxx性猛交| 午夜福利,免费看| 免费观看人在逋| 色精品久久人妻99蜜桃| 精品国产一区二区久久| 日韩欧美在线二视频| 亚洲精品av麻豆狂野| 一边摸一边抽搐一进一出视频| 精品国产亚洲在线| 美女午夜性视频免费| 一a级毛片在线观看| 天堂动漫精品| av视频在线观看入口| 精品欧美一区二区三区在线| 亚洲午夜精品一区,二区,三区| 久久国产精品男人的天堂亚洲| 亚洲欧美一区二区三区黑人| 正在播放国产对白刺激| 老熟妇仑乱视频hdxx| 国产黄a三级三级三级人| 久久香蕉激情| 不卡av一区二区三区| 欧美在线黄色| 嫁个100分男人电影在线观看| 久久人人爽av亚洲精品天堂| 伊人久久大香线蕉亚洲五| or卡值多少钱| 欧美日韩亚洲国产一区二区在线观看| 在线观看免费视频日本深夜| 亚洲精品粉嫩美女一区| 性欧美人与动物交配| 757午夜福利合集在线观看| 久热这里只有精品99| 亚洲人成电影观看| 亚洲国产欧美网| 欧美日本亚洲视频在线播放| 日本 欧美在线| 久久国产精品人妻蜜桃| www.自偷自拍.com| 90打野战视频偷拍视频| 精品欧美国产一区二区三| 99香蕉大伊视频| 神马国产精品三级电影在线观看 | av片东京热男人的天堂| av天堂久久9| 女性被躁到高潮视频| 成人av一区二区三区在线看| 黄网站色视频无遮挡免费观看| 成人国语在线视频| 成人18禁在线播放| 女性生殖器流出的白浆| 日日干狠狠操夜夜爽| 黄片大片在线免费观看| bbb黄色大片| 亚洲男人的天堂狠狠| 91精品国产国语对白视频| 777久久人妻少妇嫩草av网站| 国产av在哪里看| 天堂动漫精品| a级毛片在线看网站| 99久久精品国产亚洲精品| 亚洲欧美日韩高清在线视频| 精品福利观看| 伊人久久大香线蕉亚洲五| 色综合亚洲欧美另类图片| 一夜夜www| 国产97色在线日韩免费| 日韩中文字幕欧美一区二区| 黄色视频,在线免费观看| 大码成人一级视频| 色哟哟哟哟哟哟| 亚洲av电影在线进入| 老司机午夜福利在线观看视频| 老司机在亚洲福利影院| 久久久久久国产a免费观看| e午夜精品久久久久久久| 亚洲一区中文字幕在线| 国产黄a三级三级三级人| 十分钟在线观看高清视频www| 日韩精品中文字幕看吧| 变态另类丝袜制服| 成人三级做爰电影| 国产高清视频在线播放一区| cao死你这个sao货| 国产1区2区3区精品| 亚洲成人精品中文字幕电影| 国产av又大| 亚洲中文日韩欧美视频| 精品人妻在线不人妻| 午夜福利视频1000在线观看 | 91麻豆精品激情在线观看国产| 男人操女人黄网站| 手机成人av网站| 日韩欧美三级三区| 久久久久国产一级毛片高清牌| 亚洲欧美日韩另类电影网站| 亚洲性夜色夜夜综合| 亚洲男人的天堂狠狠| 亚洲精品中文字幕在线视频| 亚洲av熟女| 亚洲情色 制服丝袜| 婷婷精品国产亚洲av在线| 人成视频在线观看免费观看| 久久天堂一区二区三区四区| 国产国语露脸激情在线看| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲免费av在线视频| 激情在线观看视频在线高清| 成人免费观看视频高清| 国产aⅴ精品一区二区三区波| 亚洲专区中文字幕在线| 男女下面插进去视频免费观看| 9色porny在线观看| 成人三级黄色视频| 欧美色视频一区免费| 老司机深夜福利视频在线观看| 精品乱码久久久久久99久播| 精品日产1卡2卡| 欧美亚洲日本最大视频资源| 一级片免费观看大全| 精品国内亚洲2022精品成人| 久久婷婷人人爽人人干人人爱 | 日本五十路高清| 狠狠狠狠99中文字幕| 国产在线精品亚洲第一网站| 国产一区在线观看成人免费| 99久久国产精品久久久| 亚洲中文字幕日韩| 18禁美女被吸乳视频| 久久午夜综合久久蜜桃| 国产又爽黄色视频| 国产精品美女特级片免费视频播放器 | 亚洲国产欧美网| 日韩欧美一区视频在线观看| 亚洲一区高清亚洲精品| 亚洲专区中文字幕在线| 久久性视频一级片| 免费无遮挡裸体视频| 亚洲精品国产精品久久久不卡| 婷婷丁香在线五月| 亚洲伊人色综图| 这个男人来自地球电影免费观看| 欧美精品啪啪一区二区三区| 黑丝袜美女国产一区| 亚洲国产中文字幕在线视频| 精品久久久久久,| 国产亚洲av嫩草精品影院| 亚洲一区二区三区不卡视频| 99精品欧美一区二区三区四区| 亚洲欧美激情综合另类| 在线观看免费视频日本深夜| 精品免费久久久久久久清纯| 国内精品久久久久精免费| 亚洲第一av免费看| x7x7x7水蜜桃| 又紧又爽又黄一区二区| 久久精品国产亚洲av高清一级| 午夜免费鲁丝| 欧美黑人精品巨大| 热99re8久久精品国产| 一级黄色大片毛片| 久久久久久久精品吃奶| 国产精品电影一区二区三区| 欧美激情 高清一区二区三区| 国产精品秋霞免费鲁丝片| 欧美精品亚洲一区二区| 性色av乱码一区二区三区2| 黄色a级毛片大全视频| 中文字幕人成人乱码亚洲影| 日本vs欧美在线观看视频| 亚洲五月色婷婷综合| 成在线人永久免费视频| 久久午夜综合久久蜜桃| 激情视频va一区二区三区| 长腿黑丝高跟| 久久精品国产亚洲av高清一级| 天堂√8在线中文| 高清毛片免费观看视频网站| 91大片在线观看| 亚洲电影在线观看av| 国产91精品成人一区二区三区| 日韩av在线大香蕉| 色综合欧美亚洲国产小说| 成人国产综合亚洲| 久久亚洲精品不卡| 一个人观看的视频www高清免费观看 | 国产欧美日韩一区二区三| 丝袜人妻中文字幕| 久久热在线av| 亚洲精品中文字幕在线视频| 亚洲精品国产精品久久久不卡| 极品人妻少妇av视频| 咕卡用的链子| 中文字幕人成人乱码亚洲影| 桃红色精品国产亚洲av| 国产精品电影一区二区三区| 最好的美女福利视频网| 亚洲国产精品成人综合色| 桃色一区二区三区在线观看| 国产极品粉嫩免费观看在线| 91大片在线观看| 亚洲国产高清在线一区二区三 | 啦啦啦韩国在线观看视频| 精品国产一区二区三区四区第35| 嫁个100分男人电影在线观看| 国产成人av激情在线播放| 伊人久久大香线蕉亚洲五| 黑丝袜美女国产一区| 国产高清videossex| 精品国产超薄肉色丝袜足j| 中文字幕人妻丝袜一区二区| 国产色视频综合| 黑人欧美特级aaaaaa片| 亚洲精品国产区一区二| 美女高潮喷水抽搐中文字幕| 午夜福利视频1000在线观看 | 久久天堂一区二区三区四区| 成年人黄色毛片网站| 成人欧美大片| 男女午夜视频在线观看| 99国产精品一区二区蜜桃av| 国产精品永久免费网站| 999精品在线视频| 国产成人精品在线电影| 欧美一级毛片孕妇| 国产亚洲av高清不卡| 变态另类丝袜制服| 欧美色欧美亚洲另类二区 | 亚洲性夜色夜夜综合| 人妻丰满熟妇av一区二区三区| 日韩欧美一区视频在线观看| 午夜久久久在线观看| 91成人精品电影| 国产麻豆成人av免费视频| 不卡一级毛片| 国产av又大| 亚洲五月色婷婷综合| 久久久久久久精品吃奶| 欧美不卡视频在线免费观看 | 伊人久久大香线蕉亚洲五| 免费看美女性在线毛片视频| 欧美日韩福利视频一区二区| 丝袜在线中文字幕| 后天国语完整版免费观看| 午夜激情av网站| 51午夜福利影视在线观看| 色综合婷婷激情| 19禁男女啪啪无遮挡网站| 黄片播放在线免费| 在线免费观看的www视频| www国产在线视频色| 婷婷精品国产亚洲av在线| 后天国语完整版免费观看| 亚洲熟女毛片儿| 色精品久久人妻99蜜桃| 变态另类成人亚洲欧美熟女 | 妹子高潮喷水视频| 久久九九热精品免费| 久热爱精品视频在线9| 人成视频在线观看免费观看| 午夜激情av网站| 大型av网站在线播放| 中文字幕精品免费在线观看视频| 搡老妇女老女人老熟妇| 岛国在线观看网站| 国产精品久久久久久人妻精品电影| 成人永久免费在线观看视频| 91av网站免费观看| 欧美乱色亚洲激情| 搡老妇女老女人老熟妇| 欧美激情久久久久久爽电影 | 久久精品91蜜桃| 色精品久久人妻99蜜桃| 9热在线视频观看99| av视频在线观看入口| 免费在线观看视频国产中文字幕亚洲| 成在线人永久免费视频| 一本综合久久免费| 亚洲七黄色美女视频| 国产成+人综合+亚洲专区| 日韩视频一区二区在线观看| 精品久久久精品久久久| 精品免费久久久久久久清纯| 一边摸一边抽搐一进一出视频| 欧美成人性av电影在线观看| 国产乱人伦免费视频| 色播亚洲综合网| 亚洲av五月六月丁香网| 亚洲av第一区精品v没综合| 国产伦一二天堂av在线观看| 性色av乱码一区二区三区2| netflix在线观看网站| or卡值多少钱| 无限看片的www在线观看| 欧美日韩亚洲综合一区二区三区_| 窝窝影院91人妻| 欧美日本中文国产一区发布| 精品高清国产在线一区| 亚洲人成77777在线视频| 波多野结衣巨乳人妻| 成人免费观看视频高清| 韩国av一区二区三区四区| 亚洲av电影在线进入| 亚洲av成人一区二区三| 欧美激情高清一区二区三区| 国产成+人综合+亚洲专区| 国产伦人伦偷精品视频| 9热在线视频观看99| 久久久久久久久久久久大奶| 欧美日韩精品网址| 又大又爽又粗| 久久婷婷成人综合色麻豆| 一区在线观看完整版| 亚洲一区二区三区色噜噜| 久久久国产成人免费| 久久久久久人人人人人| 日韩一卡2卡3卡4卡2021年| 国产av在哪里看| 国产区一区二久久| 国产又色又爽无遮挡免费看| 黄色丝袜av网址大全| 欧美午夜高清在线| 啦啦啦观看免费观看视频高清 | 国产精品乱码一区二三区的特点 | 亚洲精品在线美女| 日韩有码中文字幕| 一进一出抽搐gif免费好疼| 一边摸一边抽搐一进一小说| 一本大道久久a久久精品| 久久中文看片网| 高清毛片免费观看视频网站| 在线观看午夜福利视频| 精品一区二区三区av网在线观看| 大型黄色视频在线免费观看| 欧美成人一区二区免费高清观看 | 在线观看免费视频日本深夜| 日韩欧美三级三区| 日韩欧美在线二视频| 欧美精品亚洲一区二区| 久久中文字幕人妻熟女| 色尼玛亚洲综合影院| 好男人电影高清在线观看| 午夜免费激情av| 亚洲自偷自拍图片 自拍| 久久天堂一区二区三区四区| 老司机午夜福利在线观看视频| 亚洲国产欧美日韩在线播放| 中亚洲国语对白在线视频| 9191精品国产免费久久| 国产精品久久久人人做人人爽| 中文字幕色久视频| 91字幕亚洲| 麻豆av在线久日| 9191精品国产免费久久| 国产91精品成人一区二区三区| 最新在线观看一区二区三区| 黄色 视频免费看| 欧美成人性av电影在线观看| 热re99久久国产66热| 在线免费观看的www视频| 色综合欧美亚洲国产小说| 波多野结衣高清无吗| 制服诱惑二区| 少妇 在线观看| 18禁美女被吸乳视频| 亚洲av片天天在线观看| 国产一区二区三区综合在线观看| 国产激情久久老熟女| 国产激情欧美一区二区| 亚洲精品一区av在线观看| 亚洲人成电影观看| 天天一区二区日本电影三级 | 性欧美人与动物交配| 一区二区三区激情视频| 日韩av在线大香蕉| 无限看片的www在线观看| 琪琪午夜伦伦电影理论片6080| 欧美久久黑人一区二区| 精品人妻在线不人妻| 99久久综合精品五月天人人| av中文乱码字幕在线| 午夜精品久久久久久毛片777| 每晚都被弄得嗷嗷叫到高潮| 亚洲狠狠婷婷综合久久图片| 欧美精品啪啪一区二区三区| 午夜免费激情av| 纯流量卡能插随身wifi吗| 亚洲欧美激情综合另类| 一区福利在线观看| 午夜成年电影在线免费观看| 国产野战对白在线观看| 极品人妻少妇av视频| 午夜精品国产一区二区电影| 午夜福利视频1000在线观看 | 色精品久久人妻99蜜桃| 麻豆久久精品国产亚洲av| 亚洲人成伊人成综合网2020| 久久国产精品男人的天堂亚洲| 正在播放国产对白刺激| 女性生殖器流出的白浆| 亚洲最大成人中文| 国产亚洲av嫩草精品影院| 91字幕亚洲| 欧美一级毛片孕妇| 在线视频色国产色| 在线观看午夜福利视频| aaaaa片日本免费| 999精品在线视频| 制服诱惑二区| 国产一级毛片七仙女欲春2 | 国产三级黄色录像| 黄片播放在线免费| 男人舔女人下体高潮全视频| 欧美日韩一级在线毛片| 国产视频一区二区在线看| 91av网站免费观看| 午夜日韩欧美国产| 午夜福利,免费看| 精品久久蜜臀av无| 此物有八面人人有两片| 一区在线观看完整版| 老司机在亚洲福利影院| 99国产综合亚洲精品| 亚洲 欧美一区二区三区| 欧美黑人精品巨大| 亚洲一区高清亚洲精品| 午夜福利高清视频| 欧美日韩福利视频一区二区| 老汉色∧v一级毛片| 男人操女人黄网站| 亚洲一区二区三区色噜噜| 大香蕉久久成人网| 欧美日韩黄片免| 真人一进一出gif抽搐免费| 国产欧美日韩一区二区三区在线| 真人一进一出gif抽搐免费| 亚洲国产高清在线一区二区三 | 久久久国产精品麻豆| 91九色精品人成在线观看| 亚洲欧美一区二区三区黑人| 女人爽到高潮嗷嗷叫在线视频| 久久精品影院6| 99国产精品99久久久久| 一个人免费在线观看的高清视频| 在线国产一区二区在线| 天堂√8在线中文| 窝窝影院91人妻| 成人免费观看视频高清| 99精品在免费线老司机午夜| 国产一区二区在线av高清观看| 人妻丰满熟妇av一区二区三区| 美女高潮喷水抽搐中文字幕| 精品久久蜜臀av无| 啦啦啦免费观看视频1| 国产精品 国内视频| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精华一区二区三区| 成年版毛片免费区| 国产高清有码在线观看视频 | 变态另类丝袜制服| 久久国产精品人妻蜜桃| 美女高潮喷水抽搐中文字幕| 91老司机精品| 亚洲成av片中文字幕在线观看| 淫妇啪啪啪对白视频| 99国产精品免费福利视频| 精品国内亚洲2022精品成人| 欧美日韩中文字幕国产精品一区二区三区 | 级片在线观看| 亚洲av美国av| 亚洲欧美日韩高清在线视频| 变态另类丝袜制服| 波多野结衣巨乳人妻| 国产精品精品国产色婷婷| 精品久久久精品久久久| 亚洲精品中文字幕一二三四区| 成年版毛片免费区| 又大又爽又粗| 亚洲精品国产区一区二| 中文字幕人妻熟女乱码| 日本五十路高清| 中出人妻视频一区二区| 亚洲精品久久成人aⅴ小说| 色尼玛亚洲综合影院| 久热这里只有精品99| 级片在线观看| 久久中文看片网| 精品久久久久久久久久免费视频| 亚洲熟妇中文字幕五十中出| 日韩欧美国产在线观看| 精品人妻在线不人妻| 久久伊人香网站| 免费av毛片视频| 我的亚洲天堂| videosex国产| 久久精品亚洲精品国产色婷小说| 两个人视频免费观看高清| 老司机在亚洲福利影院| 国产成人影院久久av| 亚洲精品中文字幕在线视频| 久久性视频一级片| 亚洲一区二区三区色噜噜| 久久香蕉激情| 久久久久亚洲av毛片大全| 国产区一区二久久| 日本欧美视频一区| svipshipincom国产片| 免费在线观看完整版高清| 亚洲色图综合在线观看| 国产成人精品在线电影| 精品久久久久久久人妻蜜臀av | 香蕉丝袜av| 精品国产国语对白av| 日韩三级视频一区二区三区| 国产99久久九九免费精品| 国产一区二区三区在线臀色熟女| 日韩国内少妇激情av| 18美女黄网站色大片免费观看| 亚洲一区二区三区色噜噜| 非洲黑人性xxxx精品又粗又长| 亚洲美女黄片视频| 老司机福利观看| 亚洲第一青青草原| 国内精品久久久久精免费| 19禁男女啪啪无遮挡网站| 国产亚洲精品综合一区在线观看 | 午夜福利影视在线免费观看| 精品免费久久久久久久清纯| 黄色成人免费大全| 国产欧美日韩一区二区精品| 黑人巨大精品欧美一区二区蜜桃| 首页视频小说图片口味搜索| 久久狼人影院| 国产精品一区二区免费欧美| 久久久久久大精品| 757午夜福利合集在线观看| 青草久久国产| av中文乱码字幕在线| www.精华液| 中文字幕人妻丝袜一区二区| 香蕉久久夜色| 一级片免费观看大全| 国产午夜福利久久久久久| 国产国语露脸激情在线看| 中文字幕人成人乱码亚洲影| 黑人欧美特级aaaaaa片| 丝袜美腿诱惑在线| 夜夜看夜夜爽夜夜摸| 免费在线观看亚洲国产| 精品人妻在线不人妻|