• <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ì)近年來(lái)萊州灣海水環(huán)境要素的系統(tǒng)調(diào)查和基于潛在性富營(yíng)養(yǎng)化評(píng)價(jià)模式和灰色聚類分析富營(yíng)養(yǎng)化的研究未見(jiàn)報(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ì)采樣站位見(jiàn)圖1。

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

    2013年5月(春季)、8月(夏季)開(kāi)展了兩個(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)行。采樣過(guò)程中同時(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值來(lái)攝取營(yíng)養(yǎng)鹽,必然導(dǎo)致有一部分氮(對(duì)磷限制水體而言)或磷(對(duì)氮限制水體而言)相對(duì)過(guò)剩。郭衛(wèi)東等[15]認(rèn)為水體只有得到適量的氮(對(duì)磷限制水體而言)或磷(對(duì)氮限制水體而言)的補(bǔ)充,讓氮磷比(N/P)的值接近Redfield值,才能使這部分氮或磷對(duì)富營(yíng)養(yǎng)化的貢獻(xiàn)體現(xiàn)出來(lái),這種現(xiàn)象被稱為潛在性富營(yíng)養(yǎng)化。郭衛(wèi)東等以此概念為基礎(chǔ),提出了潛在性富營(yíng)養(yǎng)化評(píng)價(jià)模式,具體劃分原則和富營(yíng)養(yǎng)化分級(jí)見(jiàn)表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í)的方法和公式,參見(jiàn)文獻(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%)??梢?jiàn),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]。此外,大氣的氮沉降也是海洋中氮的重要來(lái)源[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í)別,這說(shuō)明磷限制為萊州灣營(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),地表水把未被利用的過(guò)量氮肥匯入河水,而磷酸鹽從河流徑流的補(bǔ)充有限,致使入海徑流氮磷比很高[15],高氮磷比河水的注入必然會(huì)使河口及其鄰近海域水體中N/P升高。另外,也與大氣氮沉降[27]和水交換周期長(zhǎng)有關(guān)。氮濃度過(guò)高和氮磷比失衡會(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ì)過(guò)度增加并導(dǎo)致生態(tài)系統(tǒng)有機(jī)質(zhì)增多、低氧區(qū)形成、藻華暴發(fā)一些異常改變的過(guò)程[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]?;翌惥垲惙治鐾ㄟ^(guò)白化數(shù)據(jù)的灰化處理,克服了評(píng)價(jià)標(biāo)準(zhǔn)邊界過(guò)于明確對(duì)評(píng)價(jià)結(jié)果的影響,既可較好體現(xiàn)富營(yíng)養(yǎng)化程度的灰色性和調(diào)查數(shù)據(jù)的不確定性,又可反映單個(gè)指標(biāo)或多個(gè)指標(biāo)對(duì)于污染的貢獻(xiàn),并可通過(guò)各個(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)來(lái)確定富營(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ǔ)充。

    以上說(shuō)明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)化和赤潮問(wèn)題的初步探討[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)過(guò)程與氮循環(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à)
    日本黄色日本黄色录像| 国产亚洲欧美精品永久| 一边亲一边摸免费视频| 国产亚洲最大av| 午夜免费鲁丝| 成人无遮挡网站| 好男人视频免费观看在线| 精品久久久精品久久久| 麻豆成人午夜福利视频| 国产乱来视频区| 欧美日韩亚洲高清精品| 老熟女久久久| 在线观看一区二区三区激情| 亚洲三级黄色毛片| 国产精品人妻久久久久久| 国产精品99久久久久久久久| 欧美激情极品国产一区二区三区 | 在线观看免费日韩欧美大片 | 欧美高清性xxxxhd video| 久久毛片免费看一区二区三区| 欧美区成人在线视频| 一区二区av电影网| 91在线精品国自产拍蜜月| 色视频www国产| 赤兔流量卡办理| av在线蜜桃| 爱豆传媒免费全集在线观看| 欧美日韩国产mv在线观看视频 | 国产欧美另类精品又又久久亚洲欧美| 能在线免费看毛片的网站| 久久久久久久国产电影| 岛国毛片在线播放| 亚洲美女黄色视频免费看| 我要看日韩黄色一级片| 51国产日韩欧美| 80岁老熟妇乱子伦牲交| 亚洲av欧美aⅴ国产| 又粗又硬又长又爽又黄的视频| 热99国产精品久久久久久7| 老女人水多毛片| 亚洲无线观看免费| 九九在线视频观看精品| 女性生殖器流出的白浆| 精品国产一区二区三区久久久樱花 | 精品一区二区三区视频在线| 欧美成人a在线观看| 高清欧美精品videossex| 亚洲国产最新在线播放| 久久久色成人| 久久99热6这里只有精品| 中文字幕亚洲精品专区| 老师上课跳d突然被开到最大视频| 国产色婷婷99| 80岁老熟妇乱子伦牲交| 高清视频免费观看一区二区| 成年美女黄网站色视频大全免费 | 美女主播在线视频| 啦啦啦中文免费视频观看日本| 亚洲精品中文字幕在线视频 | 久久人妻熟女aⅴ| 最近最新中文字幕免费大全7| 王馨瑶露胸无遮挡在线观看| 久久精品国产a三级三级三级| 免费av不卡在线播放| 亚洲怡红院男人天堂| av卡一久久| 久久精品国产鲁丝片午夜精品| 国产男女超爽视频在线观看| 美女国产视频在线观看| 最近手机中文字幕大全| 亚洲第一av免费看| 亚洲av日韩在线播放| 国产亚洲欧美精品永久| 最近2019中文字幕mv第一页| 80岁老熟妇乱子伦牲交| 国产极品天堂在线| 美女xxoo啪啪120秒动态图| 麻豆精品久久久久久蜜桃| 日日摸夜夜添夜夜添av毛片| 国产精品偷伦视频观看了| 青青草视频在线视频观看| 干丝袜人妻中文字幕| 国产在线视频一区二区| 男女下面进入的视频免费午夜| 午夜福利视频精品| 国产乱人偷精品视频| 99视频精品全部免费 在线| 免费不卡的大黄色大毛片视频在线观看| av不卡在线播放| 成人美女网站在线观看视频| 制服丝袜香蕉在线| 精品久久国产蜜桃| 国产精品国产三级国产专区5o| 色5月婷婷丁香| 一区在线观看完整版| 久久 成人 亚洲| 亚洲不卡免费看| 国产av码专区亚洲av| 一本色道久久久久久精品综合| 大香蕉97超碰在线| 男女无遮挡免费网站观看| 国产乱人视频| 一级二级三级毛片免费看| 国产精品久久久久久精品古装| 久久久久国产网址| 久久99热这里只有精品18| 春色校园在线视频观看| 男人添女人高潮全过程视频| 亚洲精品乱码久久久久久按摩| 激情 狠狠 欧美| 婷婷色av中文字幕| 网址你懂的国产日韩在线| 高清日韩中文字幕在线| 久久久色成人| 久久久欧美国产精品| 91aial.com中文字幕在线观看| 日韩欧美精品免费久久| 国产日韩欧美亚洲二区| 欧美xxⅹ黑人| 麻豆国产97在线/欧美| 欧美xxxx性猛交bbbb| 亚洲性久久影院| 少妇丰满av| 亚洲精品国产av成人精品| 中文字幕久久专区| 亚洲一区二区三区欧美精品| 校园人妻丝袜中文字幕| 国产在线一区二区三区精| 色网站视频免费| 国产成人精品婷婷| 国产中年淑女户外野战色| 人人妻人人爽人人添夜夜欢视频 | 欧美日韩在线观看h| 亚洲欧美精品专区久久| 国产欧美亚洲国产| 伦精品一区二区三区| 久久久久久九九精品二区国产| 欧美日韩视频精品一区| 国产精品精品国产色婷婷| 免费观看的影片在线观看| 国产精品伦人一区二区| 中文字幕人妻熟人妻熟丝袜美| 赤兔流量卡办理| 亚洲成人av在线免费| 国产精品国产三级专区第一集| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲av.av天堂| 久久久久久人妻| 国产亚洲av片在线观看秒播厂| 夜夜爽夜夜爽视频| 国产免费一级a男人的天堂| 亚洲国产精品专区欧美| 在线观看美女被高潮喷水网站| 春色校园在线视频观看| 免费少妇av软件| 人体艺术视频欧美日本| 97超视频在线观看视频| 精品久久久精品久久久| 在线亚洲精品国产二区图片欧美 | 18禁在线播放成人免费| 在现免费观看毛片| 久久久久久人妻| 日韩国内少妇激情av| 国产精品国产三级专区第一集| 日韩av免费高清视频| 大香蕉97超碰在线| 乱码一卡2卡4卡精品| 成人亚洲精品一区在线观看 | 亚洲va在线va天堂va国产| 国产久久久一区二区三区| www.av在线官网国产| 亚洲美女黄色视频免费看| 久久影院123| 在线免费十八禁| 免费少妇av软件| 欧美三级亚洲精品| 一本—道久久a久久精品蜜桃钙片| 最近手机中文字幕大全| 国产成人a区在线观看| 国产av码专区亚洲av| 日韩制服骚丝袜av| 日日摸夜夜添夜夜爱| 中文字幕制服av| 亚洲久久久国产精品| 天天躁夜夜躁狠狠久久av| 久久久久久久久大av| 校园人妻丝袜中文字幕| 国产视频内射| 80岁老熟妇乱子伦牲交| 91久久精品电影网| 丰满人妻一区二区三区视频av| 毛片一级片免费看久久久久| 国产精品蜜桃在线观看| 国产精品不卡视频一区二区| 亚洲第一区二区三区不卡| 亚洲怡红院男人天堂| 在线观看免费高清a一片| 亚洲成人一二三区av| 精品少妇黑人巨大在线播放| 在线观看av片永久免费下载| 久久ye,这里只有精品| 亚洲人与动物交配视频| 国产淫语在线视频| 午夜福利网站1000一区二区三区| 国产亚洲午夜精品一区二区久久| 一区二区三区乱码不卡18| 涩涩av久久男人的天堂| 一个人看的www免费观看视频| 亚洲一级一片aⅴ在线观看| 五月天丁香电影| 久久精品国产a三级三级三级| 色5月婷婷丁香| 黄片无遮挡物在线观看| 亚洲国产精品一区三区| 久久精品国产亚洲网站| 欧美日本视频| 国产黄频视频在线观看| 日韩成人av中文字幕在线观看| 欧美3d第一页| 伦精品一区二区三区| 欧美日韩视频精品一区| 久久99热这里只频精品6学生| 亚洲精品国产av蜜桃| 岛国毛片在线播放| 免费播放大片免费观看视频在线观看| 免费大片18禁| 超碰97精品在线观看| 日日摸夜夜添夜夜添av毛片| 日本-黄色视频高清免费观看| 亚洲欧洲国产日韩| videossex国产| 日韩av在线免费看完整版不卡| 最近中文字幕2019免费版| 在线观看一区二区三区| 色视频在线一区二区三区| 如何舔出高潮| 亚洲色图av天堂| 日日撸夜夜添| 精品视频人人做人人爽| 亚洲精品乱码久久久久久按摩| 舔av片在线| 蜜桃亚洲精品一区二区三区| 夫妻性生交免费视频一级片| 国产精品国产三级国产av玫瑰| 亚洲av免费高清在线观看| 在线观看免费视频网站a站| 男人舔奶头视频| 久久99热这里只频精品6学生| 亚洲成人一二三区av| 国产成人一区二区在线| 日韩制服骚丝袜av| 超碰97精品在线观看| 日韩欧美 国产精品| 美女内射精品一级片tv| 国产精品99久久久久久久久| 九色成人免费人妻av| 久久精品国产a三级三级三级| av在线观看视频网站免费| 日本欧美国产在线视频| 女性生殖器流出的白浆| 国产精品精品国产色婷婷| 国产日韩欧美在线精品| 久久精品久久久久久噜噜老黄| 国产淫片久久久久久久久| 免费看不卡的av| 人妻制服诱惑在线中文字幕| 精品一品国产午夜福利视频| 男女免费视频国产| 黑人猛操日本美女一级片| 在线免费观看不下载黄p国产| 久久鲁丝午夜福利片| 国产免费福利视频在线观看| av免费在线看不卡| 日韩亚洲欧美综合| 另类亚洲欧美激情| 亚洲欧美一区二区三区黑人 | 亚洲丝袜综合中文字幕| 亚洲av国产av综合av卡| 黄色怎么调成土黄色| 欧美高清性xxxxhd video| 国产有黄有色有爽视频| 国产亚洲欧美精品永久| 人妻系列 视频| 岛国毛片在线播放| 最近最新中文字幕免费大全7| 在线免费十八禁| 国产伦在线观看视频一区| 亚洲av成人精品一二三区| 热99国产精品久久久久久7| 国产午夜精品一二区理论片| 中文欧美无线码| 插阴视频在线观看视频| 亚洲色图av天堂| 超碰av人人做人人爽久久| 在线观看免费视频网站a站| 十分钟在线观看高清视频www | 国产日韩欧美在线精品| 国产精品秋霞免费鲁丝片| 一级av片app| 久久久午夜欧美精品| 女性被躁到高潮视频| 国产av精品麻豆| 免费观看无遮挡的男女| 国产精品99久久久久久久久| 国产精品一区二区在线不卡| av专区在线播放| 国产一区亚洲一区在线观看| 极品少妇高潮喷水抽搐| 日韩精品有码人妻一区| 91精品国产九色| 自拍偷自拍亚洲精品老妇| 亚洲欧美一区二区三区国产| 两个人的视频大全免费| 2021少妇久久久久久久久久久| 人人妻人人澡人人爽人人夜夜| 亚洲精品久久午夜乱码| 国产精品不卡视频一区二区| 男女国产视频网站| h视频一区二区三区| 日韩三级伦理在线观看| 观看免费一级毛片| 日日撸夜夜添| 丝瓜视频免费看黄片| 一区在线观看完整版| 国产亚洲一区二区精品| 日韩av在线免费看完整版不卡| 久久久久网色| 三级国产精品欧美在线观看| 国产在视频线精品| 欧美日韩综合久久久久久| 最近中文字幕高清免费大全6| 亚洲国产精品999| 国产视频首页在线观看| 亚洲av中文字字幕乱码综合| 又粗又硬又长又爽又黄的视频| 最近最新中文字幕大全电影3| 哪个播放器可以免费观看大片| 免费在线观看成人毛片| 亚洲av中文av极速乱| 国产 精品1| 人妻夜夜爽99麻豆av| 免费观看的影片在线观看| 99精国产麻豆久久婷婷| 日韩欧美精品免费久久| 欧美区成人在线视频| 国产免费一区二区三区四区乱码| 人人妻人人澡人人爽人人夜夜| 91久久精品国产一区二区成人| 国产成人freesex在线| 偷拍熟女少妇极品色| 中国美白少妇内射xxxbb| 久久ye,这里只有精品| 小蜜桃在线观看免费完整版高清| 久久久精品94久久精品| 国产精品国产av在线观看| 国国产精品蜜臀av免费| 伊人久久国产一区二区| 2022亚洲国产成人精品| 噜噜噜噜噜久久久久久91| 丝袜脚勾引网站| 久久亚洲国产成人精品v| 国产成人精品婷婷| 国产在线一区二区三区精| 91午夜精品亚洲一区二区三区| 欧美xxxx黑人xx丫x性爽| 99re6热这里在线精品视频| 免费av中文字幕在线| 青春草国产在线视频| h日本视频在线播放| 中文精品一卡2卡3卡4更新| 成人二区视频| av在线蜜桃| 一二三四中文在线观看免费高清| 亚洲第一av免费看| 欧美日韩精品成人综合77777| 男人爽女人下面视频在线观看| 五月开心婷婷网| 夫妻午夜视频| 涩涩av久久男人的天堂| 日本与韩国留学比较| 国产亚洲精品久久久com| 久久久久国产网址| a级毛片免费高清观看在线播放| 伊人久久国产一区二区| 欧美亚洲 丝袜 人妻 在线| 国产精品久久久久成人av| 十八禁网站网址无遮挡 | 中国国产av一级| 亚洲三级黄色毛片| 九九久久精品国产亚洲av麻豆| 亚洲精品第二区| 欧美+日韩+精品| 美女中出高潮动态图| 亚洲激情五月婷婷啪啪| 亚洲精品一二三| 国产亚洲91精品色在线| 中文字幕亚洲精品专区| 交换朋友夫妻互换小说| 亚洲成人手机| 国产精品不卡视频一区二区| 18禁裸乳无遮挡动漫免费视频| 欧美日韩亚洲高清精品| 亚洲av国产av综合av卡| av不卡在线播放| 女人久久www免费人成看片| 亚洲精品国产成人久久av| 涩涩av久久男人的天堂| 少妇人妻一区二区三区视频| 乱系列少妇在线播放| 全区人妻精品视频| 嫩草影院入口| 日本免费在线观看一区| 亚洲精品日本国产第一区| 中文天堂在线官网| 午夜激情福利司机影院| av播播在线观看一区| 亚洲成人一二三区av| 看非洲黑人一级黄片| 亚洲精品乱码久久久v下载方式| 美女国产视频在线观看| 亚洲三级黄色毛片| a级毛片免费高清观看在线播放| 中文资源天堂在线| 日日摸夜夜添夜夜爱| av天堂中文字幕网| 热99国产精品久久久久久7| 久久久久久九九精品二区国产| 少妇熟女欧美另类| 日韩中字成人| videos熟女内射| 又粗又硬又长又爽又黄的视频| 80岁老熟妇乱子伦牲交| 日本av手机在线免费观看| h日本视频在线播放| 夜夜骑夜夜射夜夜干| 1000部很黄的大片| 国产视频内射| 少妇被粗大猛烈的视频| 老女人水多毛片| 亚洲内射少妇av| 中文资源天堂在线| 美女福利国产在线 | 国产精品人妻久久久影院| 午夜免费观看性视频| 美女xxoo啪啪120秒动态图| 伊人久久精品亚洲午夜| 一级毛片我不卡| 午夜福利网站1000一区二区三区| 久久99热6这里只有精品| 在线观看一区二区三区| 又黄又爽又刺激的免费视频.| 久久精品国产亚洲av涩爱| 国产精品一区二区在线不卡| 伦理电影大哥的女人| 一级毛片我不卡| av网站免费在线观看视频| 国产亚洲一区二区精品| 精品久久久久久电影网| 久久精品熟女亚洲av麻豆精品| 舔av片在线| 大香蕉久久网| 大香蕉97超碰在线| av在线观看视频网站免费| 中文字幕免费在线视频6| 国产精品国产三级国产av玫瑰| a级一级毛片免费在线观看| 色视频在线一区二区三区| 久久国产精品大桥未久av | 联通29元200g的流量卡| 少妇裸体淫交视频免费看高清| 波野结衣二区三区在线| 亚洲国产色片| 亚洲天堂av无毛| 高清在线视频一区二区三区| 97在线视频观看| 亚洲最大成人中文| 亚洲综合精品二区| freevideosex欧美| 亚洲图色成人| 97在线视频观看| 久久鲁丝午夜福利片| 日韩免费高清中文字幕av| 天天躁夜夜躁狠狠久久av| 国产v大片淫在线免费观看| av免费在线看不卡| 亚洲精华国产精华液的使用体验| 亚洲aⅴ乱码一区二区在线播放| 日本免费在线观看一区| 国产精品国产三级国产专区5o| 日韩不卡一区二区三区视频在线| 一级毛片aaaaaa免费看小| 综合色丁香网| 欧美三级亚洲精品| 日本黄色片子视频| 色吧在线观看| 欧美老熟妇乱子伦牲交| videossex国产| 在线观看av片永久免费下载| 久久久久久久国产电影| 国产亚洲欧美精品永久| 国产在线视频一区二区| 在现免费观看毛片| 欧美日韩视频精品一区| 亚洲av成人精品一区久久| 丰满迷人的少妇在线观看| 欧美最新免费一区二区三区| 久久精品夜色国产| 国产伦精品一区二区三区视频9| 亚洲丝袜综合中文字幕| 一二三四中文在线观看免费高清| 边亲边吃奶的免费视频| 在线亚洲精品国产二区图片欧美 | 亚洲综合色惰| 蜜桃在线观看..| 久久午夜福利片| 亚洲精品乱码久久久久久按摩| 精品99又大又爽又粗少妇毛片| 一级片'在线观看视频| 日本黄色日本黄色录像| 亚洲精品日本国产第一区| 亚洲天堂av无毛| 亚洲第一av免费看| 久久99热这里只频精品6学生| 午夜福利在线观看免费完整高清在| 日本黄色片子视频| 各种免费的搞黄视频| 亚洲无线观看免费| 久久国内精品自在自线图片| 国产真实伦视频高清在线观看| 又黄又爽又刺激的免费视频.| 国产精品女同一区二区软件| 欧美三级亚洲精品| 亚洲精品456在线播放app| 爱豆传媒免费全集在线观看| 久久久久久久久久成人| 99久国产av精品国产电影| 日本av免费视频播放| 亚洲人成网站高清观看| 综合色丁香网| 人妻夜夜爽99麻豆av| 久久久久久久大尺度免费视频| 久久久国产一区二区| www.av在线官网国产| 亚洲国产精品成人久久小说| 又大又黄又爽视频免费| 少妇的逼水好多| 欧美成人a在线观看| 亚洲一区二区三区欧美精品| 97在线人人人人妻| 一级二级三级毛片免费看| 亚洲电影在线观看av| 成年免费大片在线观看| 国产精品人妻久久久影院| 国产在线视频一区二区| 下体分泌物呈黄色| 国产伦精品一区二区三区视频9| 免费看不卡的av| 午夜福利影视在线免费观看| 亚洲欧美日韩无卡精品| 成人美女网站在线观看视频| av在线播放精品| 22中文网久久字幕| 国产久久久一区二区三区| 久久ye,这里只有精品| 插逼视频在线观看| 国产一级毛片在线| 有码 亚洲区| 欧美最新免费一区二区三区| 国产精品一二三区在线看| 日韩三级伦理在线观看| 免费大片黄手机在线观看| 国产精品99久久99久久久不卡 | 免费播放大片免费观看视频在线观看| 成人午夜精彩视频在线观看| 能在线免费看毛片的网站| 黄片无遮挡物在线观看| 亚洲三级黄色毛片| 欧美成人a在线观看| 纵有疾风起免费观看全集完整版| 99热全是精品| 少妇被粗大猛烈的视频| 老师上课跳d突然被开到最大视频| 亚洲成人av在线免费| 日韩中文字幕视频在线看片 | 久久热精品热| 妹子高潮喷水视频| 最新中文字幕久久久久| 九九爱精品视频在线观看| 丰满人妻一区二区三区视频av| 一级毛片aaaaaa免费看小| 亚洲国产精品一区三区| 成人免费观看视频高清| 国产精品久久久久成人av| 日韩av在线免费看完整版不卡| 男人狂女人下面高潮的视频| 少妇人妻久久综合中文| av专区在线播放| 视频中文字幕在线观看| 黄色日韩在线| 亚洲国产精品成人久久小说| 免费人妻精品一区二区三区视频| 欧美丝袜亚洲另类| 精品久久久久久久久亚洲| 又爽又黄a免费视频| 亚洲自偷自拍三级| 日韩在线高清观看一区二区三区| 久久99热这里只频精品6学生| 人体艺术视频欧美日本| 黄色视频在线播放观看不卡| 深夜a级毛片| 久久久久久久久久久丰满| 一边亲一边摸免费视频| 高清视频免费观看一区二区| 精品久久久精品久久久| 欧美精品人与动牲交sv欧美|